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ABSTRACT 

We have measured the differential energy spectra of cosmic-ray 
positrons and negatrons with energies between ~11 and 1500 MeV during 
the period 1968-1971 using a balloon-borne magnetic spectrometer. These 

measurements fill a gap in the previously existing data and permit us 
to determine, within quantitative limits, the interstellar spectra of 
cosmic-ray positrons and electrons (e + + e ). Knowledge of these 
spectra provides a crucial tool for studies of the distribution and 
density of matter and magnetic fields in the interstellar medium and 
the origin and dynamics of energetic particles contained in the fields. 

From a study of the near -Earth electron spectra and their 
relationship to the interstellar spectrum derived from the galactic 
non- thermal-radio-background emission, and from a study of the near- 
Earth positron spectra and their relationship to the interstellar 
positron spectrum calculated from collisions of cosmic-ray nuclei with 
the interstellar matter, we have found that the differential energy 
spectrum of interstellar electrons may be represented as a power-law, 
j a T for 100 MeV ~ T ~ 2 GeV, but must flatten considerably at 
lower energies. From the measured electron charge composition, which 
we find to be little affected by solar modulation, we have concluded 
that the majority of cosmic-ray electrons with energies above ~10 MeV 
are not the result of nuclear collisions in the galaxy but presumably 
originate in "primary" sources. 

In the energy range of our measurements the near-Earth 
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intensities of cosmic-ray positrons and electrons, as well as the 
intensity of cosmic-ray nuclei, are significantly lower than their 
interstellar intensities because the particles are scattered by magnetic 
irregularities imbedded in the outward- flowing plasma of the solar 
wind. Long-term changes in the scattering properties of the inter- 
planetary medium, i.e. in the cosmic-ray diffusion coefficient, k, are 
responsible for the observed long-term variations in the near-Earth 
cosmic -ray intensities which are as large as a factor of 10 from 
"solar minimum" to "solar maximum". We have used the cosmic-ray 
positron and electron spectra as tools to study the solar modulation 
mechanism. By using numerical solutions of the cosmic-ray transport 
equation to relate the near-Earth electron spectra to the interstellar 
electron spectrum, we have found that the magnetic rigidity dependence 
of the interplanetary cosmic-ray diffusion coefficient at rigidities 
from ~100 MV to ~10 GV may be represented as k a with b increasing 
from 0 to ~l-2 with increasing rigidity. However, from a comparison 
of the near-Earth and interstellar positron spectra we find that below 
~60 MV the diffusion coefficient must increase with decreasing rigidity. 

The magnitude of the diffusion coefficient at 1 AU derived 
from the electron and positron modulation studies depends on the 
assumed radial dependence of k. In order to place limits on this 
radial dependence and to make estimates of the size of the solar 
modulation region, we have also evaluated diffusion coefficients from 
measurements of the power spectrum of the interplanetary magnetic field 
near 1 AU. Assuming K( r ) Q r 11 , we have found that n ~ 1.1 in order 
that the calculated modulation beyond 1 AU agretfwith the observed 
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modulation. For k independent of radius, we obtained consistency 
between the diffusion coefficients derived by the two methods for 
boundary distances of the solar modulation region in the range of 
6-25 AU. . , 

These diffusion coefficients derived from the electron 
modulation study must also apply to the cosmic-ray nuclei. As a 
consistency check, we have used the electron diffusion coefficients 
to calculate solutions of the transport equation for cosmic-ray 
protons and He nuclei for four different time periods from 1965 to 
1970. Assuming a particular, time-independent form for the interstellar 
spectra of these particles, we have derived spectra at 1 AU which are 
consistent with the observations over the full range of intensity 
variations observed during this solar half cycle. 
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I. INTRODUCTION 

Cosmic-ray electrons were known to exist long before their 
discovery near Earth in 1960 by Earl (1961) and Meyer and Vogt (1961). 
Radio astronomers have observed the synchrotron radiation from 
relativistic electrons in such places as the Sun, Jupiter, the inter- 
stellar medium, supernovae envelopes, and other galaxies. Thus they 
are almost universal in nature. 

Because of their universality and their energy losses due 
to synchrotron radiation and inverse Compton collisions with photons, 
cosmic-ray electrons represent unique probes for determining physical 
conditions in the universe. For example, from an analysis of the 
observed galactic synchrotron background radiation, information on the 
galactic magnetic field, the structure of the interstellar medium, and 
the average interstellar electron spectrum can be obtained. On a 
larger scale, an argument for galactic confinement of the bulk of 
cosmic rays is that if cosmic-ray electrons were present in the same 
numbers throughout the universe as they are in our galaxy, then inverse 
Compton collisions with the universal black-body photons would give 
rise to an isotropic flux of x-rays far in excess of what is observed. 
Hence the observation and interpretation of the cosmic-ray electron 
flux has important implications on the distribution of matter and 
fields in both interstellar and intergalactic space. 

The origin of the electron component of cosmic rays has 
long been debated. Their existence in expanding supernovae shells 
suggests that they are directly accelerated in such sources. On the 
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other hand, collisions of the cosmic-ray nuclei with the interstellar 
matter give rise to "secondary" electrons through pion decays. Above 
~10 MeV the calculated fraction of positrons in the collision-source 
model is much higher than that observed near Earth, giving strong 
evidence for the existence of sources of directly accelerated negatrons 
(Beuermann et al., 1970, Fanselow et ali, 1969). Below ~10 MeV 
secondary knock-on negatrons outnumber those produced in nuclear inter- 
actions (Abraham et al., 1966), and the calculated intensity is 
consistent with the observed average flux (Cline and Porreca, 1970). 
However, the flux of low-energy electrons is highly variable, even 
during solar quiet times (McDonald et al. , 1972). The origin of these 
variations is uncertain, and hence the origin of the low-energy 
particles themselves remains in doubt. 

Many solutions to problems in cosmic-ray astrophysics 
depend on a knowledge of the energy spectra of the particles at their 
source. However, near Earth we observe the spectra which are modulated 
by the outward- flowing solar plasma. The study of this long-term 
modulation, which is anti-correlated with the 11 year sunspot cycle, 
has two immediate alms: 1) the determination of the local interstellar 
spectra of cosmic rays, and 2) information on the state of the inter- 
planetary medium through which the particles diffuse and lose energy. 
Although electrons comprise only a small fraction of the total cosmic- 
ray flux, the study of their modulation provides us with important 

*In this thesis the designations "positron" and "negatron" will be used 
whenever the charge sign is relevant to the discussion. The term 
"electron" will refer either to the sum e + + e“ or to the electron 
component of cosmic rays without regard to sign. 
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advantages over nuclei studies in realizing these aims. Among these 
advantages are: 

1) The possibility of independent knowledge of the inter- 
stellar spectra of electrons and positrons from the non-thermal-radio- 
background data and calculations of galactic nuclear collisions, 
respectively. For the nuclei studies we can only estimate the inter- 
stellar spectrum by extrapolating the high-energy, near-Earth data to 
low energies using an arbitrary power law. 

2) The relatively high sensitivity of the near-Earth 
electron spectrum to the interplanetary cosmic-ray diffusion coefficient 
and the interstellar spectrum. It has now been realized that the low- 
energy O^lGeV/nucleon) spectra of nuclei, where most of the available 
data fall, are shaped primarily by convection and adiabatic deceler- 
ation and are relatively insensitive to the low-energy values of both 
the interstellar spectrum and diffusion coefficient (Goldstein et al., 
1970a; Rygg and Earl, 1971; Urch and Gleeson, 1972; Garrard, 1973). 
Because electrons are relativistic in the energy region of 10 - 1000 
MeV, they lose energy through adiabatic deceleration at a slower rate 
than the nuclei. The total effect of adiabatic deceleration is 
diminished further for electrons because the diffusion coefficient is 
proportional to velocity. Hence, electrons diffuse much faster than 
nuclei of the same energy and, therefore, lose less energy in pene- 
trating to the Earth from the boundary. 

We emphasize that the study of cosmic-ray electron spectra 
involves a great number of interrelated topics in astrophysics. For 
example, the study of the modulation of electrons and positrons has 
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direct bearing on the state of the interplanetary medium and the 
interstellar intensity of both electrons and nuclei. Many galactic 
parameters, e.g. the magnetic field strength and the temperature of 
the interstellar medium, are involved in relating the interstellar 
electron intensity to synchrotron radiation in the galaxy. Similarly, 
the calculation of the negatron and positron spectra from galactic 
nuclear collisions depends on physical conditions in the interstellar 
medium. The propagation and confinement of cosmic rays also depend on 
these galactic parameters. Hence, the interpretation of the cosmic- 
ray electron flux observed near Earth has bearing on the condition of 
both local and interstellar space. 

Information on the galactic parameters used in these studies 
involves a wide variety of experimental and theoretical physics. The 
magnitude of the average magnetic field in the galaxy has been esti- 
mated to be in the range 3-5 p-gauss on the basis of the dynamical 
balance of the cosmic-ray pressure with the pressure of the galactic 
magnetic field (Parker, 1969a) . This range is in rough agreement with 
the observations of Faraday rotation and dispersion of pulsar signals. 
The dispersion measurements also yield information on the number 
density of thermal electrons in interstellar space. The temperature 
and number density of the thermal electrons in interstellar clouds 
are obtained from observations of 21 - cm absorption. These para- 
meters of the interstellar medium are important in determining the 
absorption of synchrotron radiation. 

Similarly, physical conditions in interplanetary space are 
inferred from a variety of sources. For example, observations of the 
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power spectrum of the interplanetary magnetic field yield important 
information on the cosmic-ray diffusion coefficient. Further infor- 
mation on this diffusion coefficient comes from studies of the modu- 
lation of cosmic-ray nuclei and from studies of solar-flare particle 
propagation (Lupton, 1972) . 

In this thesis we shall discuss these interrelated 
phenomena and attempt to form a consistent picture of our knowledge 
of cosmic-ray electrons. We shall make use of cosmic-ray positron 
and negatron data derived from observations with Caltech instruments 
in the range from ~11-1500 MeV over the period 1968-1971. We shall 
supplement our data with those of other experimenters to cover the 
solar half-cycle beginning in 1965. We shall discuss the use of 
numerical solutions to the equation describing particle propagation in 
the interplanetary medium in determining the parameters governing the 
modulation of electrons. To calculate the diffusion coefficient in 
the interplanetary medium necessary to explain the observations, we 
shall need a knowledge of the local interstellar spectrum of electrons 
For this purpose we reanalyze the non-thermal-radio-background data 
and derive a band of possible electron spectra above ~100 MeV. With 
these spectra we arrive at the rigidity dependence of the diffusion 
coefficient. An independent method of calculating this dependence 
is to use the power spectrum of the interplanetary magnetic field. 

We shall compare the diffusion coefficients calculated in the two ways 
to gain information on both the rigidity and radial dependences of the 
diffusion coefficient. Below ~100 MeV the radio data yield little 
information on the interstellar electron spectrum. At these energies 
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we use the near-Earth positron data and the calculated interstellar 
positron spectrum to obtain information on solar modulation and the 
interstellar spectrum. We shall also check the electron-derived solar 
modulation parameters for their applicability to cosmic-ray nuclei. 

Investigations similar to portions of the study presented 
here have been carried out by Beuermann et al. (1969, 1970) and Urch 
and Gleeson (1972). Beuermann et al. used their 1968 cosmic-ray 
positron data to discuss the absolute modulation of positrons below 
~200 MeV. The present study significantly extends their work by 
including more recent positron data and by including a detailed dis- 
cussion of solar modulation of both positrons and electrons from 
~10 MeV to 10 GeV over the solar half-cycle beginning in 1965. Urch 
and Gleeson (1972) derived the rigidity dependence of the diffusion 
coefficient above a few hundred MeV from the near-Earth electron data 
and an interstellar electron spectrum calculated from the non-thermal- 
radio-background data. These diffusion coefficients were then used in 
fitting the cosmic-ray proton and He-nuclei data. The present study 
extends their work in several ways, e.g. by 1) including a detailed 
study of the transport equation for electrons using numerical solutions, 

2) including the positron data in the study to provide information 

on the interstellar spectrum and diffusion coefficient at low energies 

3) quantitatively correlating positron and electron modulation results 
to determine consistent interstellar spectra of these particles, and 

4) quantitatively correlating power-spectra data and electron modulation 
results to derive information on both the radial and rigidity depen- 
dence of the diffusion coefficient. 
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II. DETECTOR SYSTEM 

A. Overall Description 

The positron and negatron data presented in this thesis 
were derived from observations with a balloon-borne magnetic spectro- 
meter near the top of the atmosphere. The instrument determined the 
charge sign and magnetic rigidity (momentum divided by charge) of 
particles by measuring their deflection in a magnetic field. 

Observations have been performed with the detector in two forms, 
hereafter referred to as MOD-1 and MOD- 2. MOD-1 was used in 1968 
and 1969 and has been described in detail by Rice (1970). It has a 
1000-gauss permanent magnet and an effective rigidity range of 6-200 MV. 
MOD -2 was flown in 1970 and 1971 and is identical to MOD-1 except that 

v 

it employs a 2300-gauss magnet and an additional gas Cerenkov counter. 
Its rigidity range is 15-1500 MV. A brief overall description of 
MOD-2 will be given for completeness, but emphasis will be placed only 
on the modifications to the original instrument. A schematic cross- 
section of MOD-2 is shown in Figure II-l. 

An "event" (observation of a charged particle) is defined 
by a triple coincidence between Telescope Counter $1 (Tl), Telescope 

V >✓ 

Counter $2 (T2) , and the Lucite Cerenkov Counter (LC), and the 
absence of a pulse from any of the guard counters. This coincidence 
produces the fast-gate pulse (FG) which triggers the high voltage 
to the spark chambers and initiates the data read-out cycle. The two 
4-gap spark chambers are used to define the particle's trajectory 
before and after passing through the gap' of the permanent magnet. An 
exploded view of a spark-gap module is shown in Figure II-2. The 
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wires are .0022" diameter silver-coated beryllium copper and are 
evenly spaced at 48 per inch. The active area of each module is 
5" x 9". The magnetos trie tive technique is used in determining the 
spark locations. We measure the time delay between a fiducial 
pulse (fiducial wires are located outside the active chamber area at 
each end of the module) and a subsequent spark pulse. The spatial 
resolution is approximately gaussian with a standard deviation of • 
~.008". If more than one spark is present in a chamber, the location 
of the spark nearest the pickup coil is recorded and a multiple- 
spark-indicator (MSI) bit is set. The modules are continuously 
flushed during flight by standard "spark-chamber neon" (90$ neon and 
10$ helium). An ethanol admixture acts as a quenching agent. 

The lucite Cerenkov counter was retained from the MOD-1 
version and serves two functions: 

1) It eliminates a large portion of the cosmic-ray nuclei flux, 
thereby increasing the effective live time for electron events. 

2) It eliminates approximately 96$ of the upward -moving splash 
albedo particles. 

V 

The velocity threshold for Cerenkov radiation in lucite is 0.67 c 

which corresponds to rigidity thresholds of 0.46, 845, and 1690 MV 

for electrons, protons, and alpha particles, respectively. Electronic 

data handling effectively increases these thresholds by ~15$. 

v v 

The gas Cerenkov counter (GC) was added to the MOD-1 detector 
system in order: 

1) to eliminate contamination due to cosmic-ray nuclei above L& ' 
threshold; because of the larger MOD- 2 magnet, these particles 
would be indistinguishable from high-energy electrons, and 
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2) to further discriminate against upward -moving particles. 

Its velocity threshold is 0.9984 c which corresponds to rigidity 
thresholds of 0 .0091, 16.8, and 33.6 GV for electrons, protons, and 
alpha particles, respectively. Each of the two phototubes of the 
counter (see Figure II-l) acts independently; a coincidence between 
the fast-gate pulse and the output of one phototube generates a data 
bit which is recorded 1 as part of the event's data word. 

The specially designed magnet guard counter (MA) is shown 
with the magnet in an exploded isometric projection in Figure II-3. 
The pole faces and the upper surface of the magnet are covered with a 
plastic scintillator which prevents the analysis of particles that 
might interact or scatter in the magnet. The 3-cm x 12-cm open 
passage, together with Tl and T2, determines the acceptance cone of 
the detector. 

The sides and top of the instrument, except for the 
telescope aperture, are surrounded by guard counters. These counters 
are in active anti-coincidence and eliminate particles which enter 
the detector from outside the acceptance cone and which might inter- 
act, providing particles which trigger the telescope counters. 

A general block diagram of the electronics system is shown 
in Figure II-4. An Accutron clock is used as a timing device. It 
drives a 4-bit time scaler (16 minute cycle) whose output is used 
to control the data-collection cycle. During the first 15 minutes 
(referred to as Phase A) of each cycle, particles which satisfy the 
coincidence requirements initiate a readout of the data, which 
requires 350 msec for completion. Each data word, consisting of 8 
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spark locations, the MSI bits, the two g£ bits, the time, and 
temperature, is recorded on 16-channel magnetic tape. The remaining 
minute of the 16-minute cycle is referred to as Phase B; during this 
time the normal coincidence trigger input is blocked and the following 
rates are scaled: g£, T1 a T2, T1 ~ T2 * l£, MA, and TA (sum of all 

guard counters except MA) . These rates are monitored to check counter 
performance, to detect variations in background radiation, and to 
determine detector dead time due to the guard counters. At the 
beginning of Phase B an internal trigger is generated which results 
in the application of high voltage to the chambers and the initiation 
of the readout cycle. Since no particle is normally present in the 
chambers, the spacing between the fiducials is thus recorded to 
provide a check of the digitizing circuitry. 

The atmospheric pressure during flight is recorded by a 
photobarograph, a device which photographs a Wallace-Tiernan aneroid 
barometer (FA 160), a clock, and a thermometer at 5-minute intervals. 
The barometer is calibrated before and after each flight and is 
accurate to ± 0.1 mb at 2.4 mb, our typical float altitude. Usually, 
two redundant photobarographs were flown on each flight. 

In the following we discuss in more detail the gas Cerenkov 
counter and the 2300-gauss magnet, the two major additions to the 
MOD-1 detector. 
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B. The Gas Cerenkov Counter 

The gas Cerenkov counter (see Figure II-l) contains sulfur 
hexafluoride at 2.2 atmospheres absolute pressure. This conf iguration 
has a velocity threshold of v = 0.9984c. The absolute kinetic energy 
threshold for various ! particles are: 

electrons 8.62 MeV 

muons 1 . 78 GeV 

pions 2.50 GeV 

protons 15.7 GeV 

alpha particles 62.9 GeV 

V 

The two flat mirrors serve to reflect the Cerenkov light 
into the phototube faces. The mirrors are constructed of 1/8-inch 
lucite and are aluminized on their upper surfaces. The conical 
mirrors are made of spun aluminum with their interior surfaces 
aluminized, A coating of magnesium fluoride covers all mirror 
surfaces to retard oxidation which would: otherwise cause poor 
reflectivity at ultraviolet wavelengths. The conical and flat mirrors 
are mounted on a thin aluminum basket which is not shown in the figure. 
Of these pieces, only the flat mirrors are within the acceptance cone 
of the detector. 

The -phototubes are EMI 9531 QB (Whittaker Corp., Plainview, 

N. Y.) which have quartz faces 3 1/2 inches in diameter. A 1/2-inch- 
thick fused-silica window (Corning Glass Works, Orange, California) is 
mounted in front of each phototube to protect them from the gas 
pressure. High voltage for the tubes is supplied by DC -DC converters 
(Crestronics , Crestline, California) whidh are mounted inside the 
phototube housings. In addition, electronic pulse-discr imination and 


i 
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coincidence circuitry is mounted inside the housings which also serve 
as a shield from the spark noise. A coincidence between the fast- 
gate pulse and the discriminator output of either phototube generates 
a data bit for the phototube involved. 

The counter was calibrated at the Caltech Synchrotron using 
300-MeV positrons. Forty-five incident beam directions were chosen 
to cover the acceptance cone of the detector. For each incident 
direction a pulse-height analysis was made, and from the resulting 
distribution the mean number of photoelectrons emitted from the 
cathode was determined. For each phototube this number varied from 
approximately 4-10 over the range of incident directions. The average 
over direction was about 6 photoelectrons. The electronic 
discrimination level was set just above the one photoelectron level, 
which results in an average efficiency of approximately 98$. However, 
the efficiency of the counter slowly decreases with time due to 
oxidation of the mirror surfaces. Therefore, we determined the 

_ V 

efficiency of the gas Cerenkov counter directly from the flight 
data by a procedure described in Appendix A.3.c. We found that 
the efficiency was approximately 93$ and 84$ for 1970 and 1971, 
respectively. 1 

The effective energy "threshold" is not precisely determined 
since the actual number of photoelectrons emitted from the photo- 
cathode is Poisson distributed about the mean; hence, any particle 
above the absolute velocity threshold has a finite probability of 
producing enough photoelectrons to trigger the discriminator. In 
Figure II-5 we show a plot of the Cerenkov light output for singly- 
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charged particles versus y = — ^ where W is the total energy of the 

me 

particle and m is the rest mass. Particles with y greater than about 3 
times the absolute threshold value produce more than 90$ of the output 
level of completely relativistic particles. We roughly estimate the 
"effective" threshold as 50$ of the full output. The corresponding 
effective kinetic -energy thresholds for the various particles are: 

electrons 12 MeV 

muons 2.5 GeV 

pions 3.4 GeV 

protons 22 GeV 

alpha particles 67.5 GeV 

The noise rate of the phototubes was monitored during the 
Phase B period. At float altitude typical values of the combined 
noise rate were 600/sec and 1200/sec in 1970 and 1971, respectively. 

The probability of an accidental coincidence is given approximately 
by the product of the noise rate and the sum of the widths of the 
fast-gate and the discriminator output pulses, which was about 10 ^ 
seconds. Thus these probabilities were roughly .0006 and .001 in 
1970 and 1971, respectively. 
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C . The Magnet 

The 2300-gauss magnet (Indiana General Corp . , Valparaiso, 
Ind.) used in the MOD-2 configuration is 1 shown in Figure II-3. The 
construction is similar to the 1000-gauss magnet employed in the MOD-1 
configuration as described by Rice (1970). Alnico-8 permanent magnets 
are used, and a magnetic circuit of steel reduces external fields as 
much as possible. The three orthogonal components of the field were 
measured at 1-cm intervals throughout the volume accessible to the 
particles out to a distance of 8 cm above and below the magnet. In 
Figure II-6 we show a plot of these three components along three 
representative paths through the magnet gap. The locations refer to 
a right-handed coordinate system with the origin at the center of the 
gap. The z-axis is vertical and positive upward; the x-axis is 
perpendicular to the ipole faces and positive toward the south pole 
(see Figure II-3) . The magnetic field was monitored before and after 
each flight by a permanently mounted Hall effect device (F . W. Bell, 
Inc., Columbus, Ohio); no change in the field strength greater than 
~10 gauss was noted on any of the flights. 

The geometrical factor of the MOD-2 detector was determined 
at 7 different rigidities between 12 and 400 MV by the Monte-Carlo 
method described by Rice (1970). In Figure II-7 we show the geo- 
metrical factor as a function of rigidity. The error bars represent 
the fluctuation due to the finite number (1000) of valid trajectories 
used in the calculation. 

The deflection of a particle of rigidity R (MV) in a 
magnetic field B is given by 
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0(rad) = 


3 x 10 
R 



(II-l) 


where B (gauss) is the component of the field normal to the trajectory 
x 

and d^(cm) is an increment of distance along the path. The line 
integral in equation II-l is referred to as the magnetic path, M. 

The Monte-Carlo program that calculates the geometrical factor also 
computes the value of M along each particle trajectory. In Figure II-8 
we show the mean values of R times 9 in MV-radians at 7 rigidities 
between 12 and 400 MV. The solid error bars refer to the rms 
deviation and the dashed error bars represent the extreme values. At 
each rigidity the value of R0 is within 2$ of 8.85 MV-radians. The 
rms deviation is typically less than 2$ of the mean and the extreme 
values are within about 7 <f of the mean. Since the resolution of the 
detector, FWHM (see Appendix A. 2) is ^ 16$ we can use the mean value 
with negligible error. Thus we use the following approximate 
relationship between deflection angle and rigidity for all 
particles : 


R = MV (MOP - 2 ) (II- 2a) 

The corresponding relationship for the MOD-1 detector, using a 1000- 
gauss magnet, is 


R = 


3.55 


MV 


(MOD-1) 


(II- 2b) 


(Rice, 1970). 

The complete sheathing of the magnet by the magnet guard 
counter eliminates particles which interact in the magnet pole faces 

i 
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and thus eliminates the necessity of detailed trajectory reconstruction. 
It is thus sufficient to read out spark locations in the y-z pro- 
jection only, which saves considerable data storage and detector 


live time. 
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III. BALLOON FLIGHTS 

The Caltech data presented in this thesis were derived from 
10 high-altitude balloon flights launched from Ft. Churchill, Manitoba 
during the summers of 1968, 1969, 1970, and 1971. (The data from 1968 
were previously published. (Beuermann et al., 1969, 1970)) We summarize 
in Table III-l the relevant information on these flights. 

In Figure III-l we show trajectories of three typical flights. 
We also show in the figure the invariant latitude contours, calculated 
from the internal field only (Cain et al;, 1967), in order to indicate 
the trajectories in the geomagnetic field. 

Figure III- 2 shows two typical altitude profiles. The solid 

curve is from flight '71C2 and is representative of the altitude profiles 

of eight of the flights. In each of these eight flights the launch was 

2 

timed so that the instrument passed through 70 g/cm altitude after the 
evening transition to low geomagnetic cutoff. This timing ensured that 
the ascent data used in the separation of atmospheric secondaries were 
not contaminated by return albedo electrons. (See Chapter IV.) The 
dashed curve in Figure III-2 is from flight 69C1. This step profile 
and a similar one from flight 69C2 were used to more accurately define 
the atmospheric depth dependence of the electron flux in 1969. 

The relationship of our flights to the 11-year solar modula- 
tion cycle is shown in Figure III-3. We have plotted the daily average 
of the hourly count rate of the Deep River neutron monitor (Steljes, 
1965-1971) for the period 1962-1971. Ground-based neutron monitors 
record the near-Earth flux of high-energy cosmic-ray nuclei 
(^ 1 GeV/nucleon) and serve as a convenient continuing reference of the 
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high-energy cosmic-ray intensity. The neutron monitor intensities reach 
a maximum during the period of minimum solar activity (1965-1966) and a 
minimum during solar maximum (1969-1970). The dates of our flights are 
marked with vertical lines in Figure III-3. The flights cover the 
period near solar maximum. We shall supplement our data with those of 
others to cover the solar half-cycle beginning in 1965. 

Since we are interested in the long-term modulation of galactic 
cosmic -ray electron spectra, it is important to identify short-term 
variations that might affect our measurements. Short-term fluctuations 
are generally associated with solar activity. The energy spectra of 
particles emitted from the sun are usually quite steep and, hence, the 
effect of solar emission is most significant at low energies. In 
addition, Forbush decreases usually follow large solar flares and produce 
a general depression of the galactic cosmic-ray flux below several GeV. 

We have examined the following sources of data relating to 
solar activity during the period of our balloon flights: 

1) ESSA bulletins (ESSA Solar Geophysical Data, 1968-1971), which 

contain, for example, data from solar proton monitors on Explorer 34 

and 41 satellites (E > 10 MeV) and on the ATS-1 satellite 

P ~ 

(Ep > 5 MeV), geomagnetic indices, and daily average of neutron 
monitor rates. 

2) Caltech cosmic-ray experiment on the 0G0-6 satellite, which provided 
information on low-energy protons and electrons from June 1969 - 
July 19 70. 

3) Goddard Space Flight Center cosmic-ray experiment on the IMP series 
of satellites, which provided almost continuous data on 3 - 12 MeV 
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electrons from November 1963 - September 1969 . (A summary of the 
electron counting rates from these experiments is given by 
McDonald et al., 1972.) 

The observations of 1968, 1969, and 1970 were made during or just after 
the recovery phases of Forbush decreases, as indicated by the neutron 
monitor counting rate (see Figure III-3). However, from examination of 
the data from the other sources, we have concluded that short-term solar 
activity, e.g. solar flares, did not contaminate our electron fluxes. 

We have also examined the Goddard Space Flight Center 3-12 MeV 
electron data for evidence of the large quiet-time increases observed 
below ~25 MeV (L'Heureux et al . , 1971; McDonald et al., 1972). We have 
concluded that our 19'68 balloon flights corresponded to quiescent flux 
levels. No comparisons could be made for the summer of 1969 since the 
published 3 - 12 MeV data extend only to mid-March 1969. However, the 
raw fluxes from different flights during the 1969 summer are not 
significantly different. Hence, we may assume that our data are 
typical of undisturbed times. 

We have made similar comparisons of the raw fluxes from 
different flights for the summers of 1968, 1970, and 1971. In the 
absence of significant differences we have combined the data from the 
flights of the same summer for greater statistical accuracy. 
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IV. DATA ANALYSIS 

The data analysis procedure for the observations made with the 

detector in the MOD-1 configuration has been previously described (Rice, 

1970). The procedure for the MOD-2 data is similar; however, the 

v 

addition of the gas Cerenkov counter and the larger magnet does require 
some new considerations. We shall give a general outline of the proce- 
dure; details of the .analysis technique are described in Appendix A. 

The basic information provided in the data word for an event 
consists of: the spark location in each of the eight spark-gap modules 

(4 above and 4 below the magnet), the multiple-spark-indicator bits, 
v 

the gas Cerenkov bits, and the time and temperature (see Chapter II) . 

In the initial phase of the analysis we sort the events according to 
spark chamber performance. The selection criteria for this sorting are 
reviewed in Appendix A. 1. a. Roughly 15 fo of the data are rejected from 
analysis in applying these criteria. 

In the process of determining chamber performance, the 
particle trajectory through each spark chamber (4 spark modules) is 
determined for the analyzable events by making a least-squares fit to 
the measured spark locations to a straight line. The bending angle 
through the magnet is then computed. A trajectory-consistency check is 
made to determine whether the calculated trajectories in the two spark 
chambers are consistent with the bending expected in the magnetic field 
for the computed deflection angle. The selection criterion established 
for this test is described in Appendix A.l.b. The criterion depends on 
the resolution of the instrument (see Appendix A. 2) and is such that 
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there is very little probability that an event with a misfit trajectory, 
etc., is accepted. Because of resolution effects approximately 11$ of 
valid MOD-1 events and about 7$ of valid MOD-2 events are also rejected. 

The deflection angle computed from the trajectories is 
inversely proportional to particle rigidity (equations II-2a and b). 

The relationship between this computed rigidity and the true particle 
rigidity involves a study of the resolution of the detector. The 
ability of the detector to measure the rigidity of a particle is affected 
primarily by 1) multiple scattering within the chambers or magnet gap, 
and 2) intrinsic angular resolution resulting from the spatial resolution 
(^.008") of each spark location measurement. Both effects, as well as 
the results of calibrations at the Caltech Synchrotron, are discussed 
in Appendix A. 2. The result is that the angular probability distribution, 
P(0, 0 ' )d0 1 = probability that a particle with rigidity R corresponding 
to deflection angle 0 will actually be observed to have deflection angle 
between 0' and 0' + d0', is approximately Gaussian, i.e., 


P(0,0') = 


a Jin 
u 


„-( 0 - 0 ' ) 2 „ 
exp [— 2~ 1 


(IV- 1) 


2a 


0 


where a is the standard deviation. From Appendix A. 2 we have: 

U 


= V( .170) 2 + ( .0025 ) 2 


0 


= •J ( . 


( .0680) 2 + ( .0025) 2 


MOD-1 


MOD- 2 


(IV- 2a) 
(IV-2b) 


The deflection resolution P, FWHM, is given by: 
2.36a, 


P s 


0 


0 


= X* o? + <^P ) 2 


MOD - 1 


(IV- 3a) 
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- vW + <^P> 2 


MOD- 2 


(IV-3b) 


A plot of the resolution is shown in Figure A- 7. Note that the effect 

of multiple scattering (constant term in equations IV-3a and b) is 

dominant below ~100 MV and is insignificant above ~500 MV. For MOD-2 

the angular resolution has a minimum value of ~16$ and increases to 

100$ at ~1500 MV (0 % .006 radians). 

Once the bending angles are determined for the analyzable 

events, the data are sorted into deflection-angle (energy) bins over 

appropriate time intervals. In determining the flux of electrons, only 

v 

data tagged with a gas Cerenkov bit are used. The raw flux in units of 
2 

particles/ (m sec sr MeV) for a given time interval is defined by: 

N. 

F a (IV-4) 


1 c L G iV^f£^i 


where : 


i = energy interval index 

= number of events in ith energy interval during 
specified time interval 
t = total live time during the time interval 

i_i 

= average geometrical factor for the ith energy 
interval 

D^ = spark chamber efficiency factor 

AT^ = width of ith energy interval 

^ .. = gas Cerenkov efficiency factor, 
eff 
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For convenience, we list in Table IV-1 the deflection-angle intervals, 
the corresponding energy intervals, and the average geometrical factors 
for both MOD-1 and MOD-2 data. The parameters t^, D^/ (typically 
0.75 - 0.85), and if (0.93 in 1970 and 0.84 in 1971) are discussed 
in Appendix A . 3 . 

In this thesis we are interested in discussing the implications 
of the flux of primary electrons, i.e. galactic particles which have 
penetrated through the interplanetary medium to 1 AU. Therefore, in 
analyzing electron fluxes observed near Ft. Churchill, Manitoba it is 
important to distinguish between fluxes at rigidities above the geo- 
magnetic cutoff rigidity, which consist of primary electrons and 
atmospheric secondaries, and those fluxes below cutoff, which consist 
of re-entrant albedo electrons and atmospheric secondaries. Recent 
calculations (Smart, 1971; Smart and Shea, 1972), based on a magneto- 
spheric model with magnetic fields of both internal and external origin, 
have shown that the geomagnetic cutoff rigidity for Ft. Churchill 
(/\ 70°) is approximately 150 MV during local daytime (~0600 to ~1800 

local magnetic time) and has an abrupt transition (due to the asymmetry 
of the magnetosphere) to a value below 20 MV for local nighttime. 

Since the intensity of re-entrant albedo electrons is found to be 

larger than that of primary electrons, this cutoff rigidity transition 

i 

is observed in balloon-borne electron detectors as a change in the 
counting rate of low-energy electrons (~ 150 MeV). Such transitions 
have been observed near Ft. Churchill in the data from the Caltech 
instrument (Rice, 1970), as well as in the data from other experiments 
(Jokipii et al., 1967; Hovestadt and Meyer, 1970; Israel and Vogt, 1969). 
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TABLE 

VI 

-1 


A. MOD-1 

Parameters 





Index 

(i) 

Deflection 
Angle Interval 
(Radians) 

Energy Interval 
at Detector 
(MeV) 

Average 

Geometrical 

Factor 






(cm^ sr) 

1 

.6 - .3 

5.4 


11.3 

2.52(e + ) 2.17(e" 

2 

.3 - .144 

11.3 

- 

24.1 

3.50 

3 

.144 - .072 

24.1 


48.8 

3.70 

4 

.072 - .036 

48.8 

- 

98.1 

3.70 

5 

.036 - .018 

98.1 

- 

197 

3.70 

* 

The asymmetry in the geometrical factor in the lowest energy range is 
due to a slight asymmetry in the geometry of the detector. 

B. MOD-2 

Parameters 





Index 

(i) 

Deflection 
Angle Interval 
(Radians) 

Energy Interval 
at Detector 
(MeV) 

Average 

Geometrical 

Factor 






(cm^ sr) 

1 

.6 - .3 

14.3 

- 

29 .0 

.2.14 

2 

.3 - .15 

29.0 

- 

58.5 

3.14 

3 

.15 - .072 

58.5 

- 

122 

3.62 

4 

.072 - .036 

122 

- 

245 

3.80 

5 

.036 - .018 

245 

- 

491 

3.80 

6 

.018 - .009 

491 

- 

983 

3.80 

7 

.009 - .006 

983 

- 

1475 

3.80 
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An example of a typical transition is shown in Figure IV-1. We show 
the hourly count rate plotted versus local time for the five lowest 
energy intervals of flight 71C2. Positrons (dotted histogram) and 
negatrons (solid histogram) are shown separately. In the low-energy 
intervals (~ 245 MeV) ; we use only the indicated nighttime (low-cutoff) 
period in deriving electron intensities. : Since no night-day transition 
was observed above 245 MeV on any of the flights of 1970 or 1971, the 
total float period is used in computing the fluxes for the three highest 
energy intervals for these years. The raw flux measurements at float 
altitude for 1968, 1969, 1970, and 1971 will be presented in Chapter V. 

These properly selected raw fluxes at float altitude consist 
of primary cosmic-ray electrons and secondary electrons generated in 
the atmosphere above the detector. We have also considered the following 
possible sources of contamination: upward-moving particles (splash 

albedo and those due to y-ray interactions in the lucite Kerenkov 
counter), atmospheric muons and pions, secondaries produced in the gas 
Cerenkov counter, high-energy cosmic-ray nuclei above gas Cerenkov 
threshold (which could be incorrectly identified as electrons), accidental 
gas Cerenkov coincidences, and spark chamber misalignment. For the 
highest energy interval of the MOD- 2 data it was necessary to make small 
corrections for contamination from high-energy nuclei, accidental gas 
Cerenkov coincidences, and spark chamber misalignment. In all other 
MOD-2 energy intervals the corrections were negligible. In the case of 
MOD-1 data, small corrections in the lower energy intervals for upward- 
moving particles were necessary. The investigation of all the sources 

I 

of contamination mentioned above is described in Appendix A. 4. The 
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values of the significant corrections are listed in the data tables 
of Chapter V . 

2 

At our typical float altitude (2.4 g/cm ) atmospheric 
secondary electrons represent a large contribution to the flux below 
a few hundred MeV. The procedure for Separating these particles from 
the primary electrons has been described in detail (Rice, 1970). We 
shall only briefly describe the method. 

The atmospheric depth dependences of the positron and 

* 

negatron rates in a given energy interval are determined from the 
data collected during ascent and descent ^ We represent these observed 
rates, r^*, by: 

r i ± (d) = a i ± s i ± (d) + b^ P i ± (d) (IV-5) 

where d is the atmospheric depth, s^ ± (d) and p^*(d) represent the 
calculated depth dependence of the rates of secondary and primary 
positrons or negatrons, respectively, and a^ and b^ are parameters 
giving the relative contribution of each component. We use the 
calculations of Beuermann (19 71) to evaluate the functions s^(d) and 
p.' t: (d). The s^(d), the secondaries generated by the nuclei component 
of the cosmic rays, are calculated using an incident nuclei spectrum 
adjusted from year to year according to changes in the Mt. Washington 


Local rates (number observed in a given energy interval per second) 
rather than fluxes (N/m2 . sec • sr *MeV) are used for convenience because 
the average geometrical factor depends on the energy dependence of the 
spectrum which changes with depth. 
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neutron monitor rate (Lockwood, private communication) and in the 
integral flux of nuclei above 400 MeV/nucleon which is measured 
directly by our detector. To calculate the p^(d) it is necessary 
to assume specific forms for the primary positron and negatron spectra 
incident at the top of the atmosphere. 

A least-squares fit is made to determine the values of a.* 

± a 4 . a + 

and b. , as well as their standard deviations a. ‘ and b. . The 
1 I 1 1 

local rate at float altitude of primary positrons and negatrons is 
then given by 


e . * = b. p.^Cd = float altitude) 
1 1 1 


(IV- 6 ) 


The rates are converted to fluxes and corrected to the top of the 
atmosphere by a procedure described in Appendix A. 5. The corrected 
fluxes are then used to estimate a new primary input spectrum in an 
iterative calculation i In practice the derived spectrum at the top 
of the atmopshere is not very sensitive to the assumed input spectrum, 
and the process converges quickly. 

The procedure described above for the separation of the 
primary and secondary components by a fitting technique is used for 
the lowest energy intervals (<; 245 MeV) where the growth curves can 
be measured with reasonable statistical accuracy. For the highest- 
energy intervals (MOD-2) these growth curves are statistically not as 
well defined. At these energies, however, the atmospheric secondaries 
correction is sufficiently small so that the secondary component can 
be calculated and simply subtracted from the observed flux to give the 
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primary contribution. 

Two examples of the fitting technique for determining the 
primary and secondary : components of the measured flux are shown in 
Figure IV-2. Figure IV-2a illustrates a case in which a relatively 
large contribution of residual primaries is obtained; Figure IV-2b 
shows a case in which zero primary flux is determined. The upturn 
at large atmospheric depths in the residual primary curve of 
Figure IV- 2a is due to the energy dependence of the incident primary 
spectrum and its changes due to energy loss in the atmosphere. 



31 


V. RESULTS 

In this capter we present the results of the analysis 
described in Chapter IV and Appendix A. Since the energy range and the 
data analysis procedure were somewhat different for the two detector 
configurations of our instrument, we discuss the results from MOD-1 
and MOD-2 observations separately. ! 

The 1968 and 1969 observations were made with the MOD-1 
detector configuration!. The absence of the gas Cerenkov counter and 
the smaller 1000-gauss bending magnet restricted the rigidity range to 
6-200 MV and required corrections at the low energies for gamma-ray- 
induced background and splash albedo. The method of correcting for 
gamma-ray interactions in the lucite Cerenkov counter has been described 
by Rice (1969). His results have been changed slightly and the esti- 
mated errors reduced as a result of further calibrations at the Caltech 
Synchrotron. In addition, the splash-albedo corrections of Rice (1970) 
have been reduced by one-third. This change in the correction was made 
after calibrations showed that electrons entering the detector from the 
backward direction had a higher probability of being rejected from 
analysis than forward -moving particles (see Appendix A. 4. a). The 
results for 1968 are shown in Table V-l and Figure V-l, and the 1969 
results are given in Table V-2 and Figure V-2. In some cases a small 
negative primary flux was obtained from the fitting procedure, indi- 
cating that the data were dominated by atmospheric secondaries. In 
these cases a l-cr upper limit above zero flux has been listed. In 
1969 the atmospheric contamination was more severe than in 1968 because 
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1) only lower float altitudes could be reached during the 1969 obser- 
vations and 2) the primary fluxes were lower due to increased solar 
modulation. For low-energy positrons the atmospheric secondaries 
contamination is smaller than for low-energy negatrons since the knock- 
on component is not present. As a result, measurable low-energy 
positron fluxes were obtained for both years. 

The 1970 and 1971 data were collected with the MOD-2 detector. 
The addition of the gas Cerenkov counter and the larger 2300-gauss 
magnet eliminated the necessity to correct for upward -moving particles 
and also allowed an extension of the rigidity range to 15-1500 MV. The 
observed fluxes were corrected to the top of the atmosphere using the 
matrix-inversion procedure described in Appendix A. 5. The 1970 results 
are shown in Tables V-3 and V-4 and Figure V-3, and the 1971 data are 
given in Tables V-5 and V-6 and Figure V-4. (The subdivision of the 
tables, corresponding to the low-energy 245 MeV) and high-energy 
data, is due to the different analysis procedures used at low and high 
energies (see Chapter IV).) 

The relatively large error limits of the data reflect the 
difficulty in measuring the charge composition of the electron spectrum 
within the atmosphere. In our energy range of 6-1500 MeV, which is of 
major interest to solar modulation studies, a magnet spectrometer, such 
as the Caltech detector, is the only instrument which can effectively 
determine this composition. Such instruments unambiguously determine 
the charge sign and offer excellent energy resolution since the electrons 
pass through little mass in traversing the magnet spectrometer. The 
upper limits and error bars of the data primarily reflect the contamina- 
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tion from atmospheric secondaries at balloon altitudes. The size and 
weight of the instrument necessary to record the relatively low flux of 
cosmic-ray electrons have made it unsuitable as a payload for 
satellites so far. Thus balloon-borne spectrometers are presently 
our only source of information on positron data in this important 
energy range. Indeed, during the period 1968-1971 the Caltech positron 
observations represent the only published data in this energy range. 

Despite the uncertainties, these positron data allow important 
definitive conclusions regarding the low-energy interstellar electron 
spectrum and the low-energy cosmic-ray diffusion coefficient. Indeed, 
below ~100 MeV the positrons represent the only direct tool for studying 
solar modulation. 
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VI. DISCUSSION OF SOLAR MODULATION 

A. Introduction and Statement of the Problems 

The study of the solar modulation of electrons provides 
information on physical conditions of the interplanetary medium, e.g. 
the cosmic-ray diffusion coefficient, as well as information on the 
local interstellar spectra of electrons and positrons. These spectra 
carry the signature of their origin, i.e. of their sources and of the 
interstellar medium in which they were stored. The electron modulation 
study also contributes to the understanding of the modulation of cosmic 
ray nuclei, providing important parameters for the deduction of their 
interstellar spectra. 

The Caltech electron data shown in Chapter V were acquired 
over the period 1968-1971 and extend over the energy range from 
~ll - 1500 MeV. We show in Figure VI-1 our 1968 and 1971 electron 
fluxes together with selected data from other authors which extend over 
the energy range 10 MeV to 10 GeV and cover the period since the last 
solar minimum in 1965-66. The effects of the long-term solar-cycle 
variation are readily apparent in the hundred MeV range in Figure Vl-la 
For example, at ~300 MeV there is about a factor of 10 difference in 
the 1965-66 and 1970 fluxes. Above ~10 GeV no distinct long-term 
variations have been observed and hence we shall ignore this region of 
the spectrum in our solar modulation study. Below ~25 MeV observations 
from detectors on the IMP and 0G0-5 satellites have shown short-term 
variations by factors of ~3-5 over time intervals of a few days 
(L'Heureux et al . , 1972; McDonald et al., 1972). These increases occur 
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during solar quiet times, are essentially energy independent from 
3-25 MeV, and are often anti-correlated with low-energy solar-proton 
events. It has been suggested that the low-energy electrons observed 
during quiet-time increases, as well as during quiescent times, are of 
galactic origin. Because of these short-term variations, the long-term 
modulation at these energies is not well-determined. McDonald et al. 
have put an upper limit of a factor of 2.3 on the intensity variation 
from solar minimum to solar maximum. It is important to note that the 
short-term variations are not observed above 25 MeV (L'Heureux et al., 
1972). I ■ 

Our observations using balloon flights, each of which lasts 
~20 hours, and which are separated by a few days, are not well suited 
to the study of the short-term variations. Ihe period of our 1968 data 
corresponds to quiet-time conditions as observed in the 3-12 MeV 
electron fluxes from the IMP-4 satellite (McDonald et al., 1972). The 
data from Simnett and McDonald (1969) for 1967 and L'Heureux et al . 

(1972) for 1968 in Figure Vl-la represent the average flux level 
during solar quiet times. This thesis, therefore, addresses itself 
to the long-term effects of solar modulation only. 

We shall present a quantitative analysis of the solar- 
modulation process which uses the electron and positron data for its 
basis but which provides a consistent picture of the modulation for all 
cosmic-ray particles. Some of the outstanding problems of solar 
modulation studies are: 

1) Interstellar Electron and Positron Spectra 

In order to make deductions on the absolute solar modulation of 




48 


electrons it is necessary to have an estimate of their interstellar 
spectrum. Previous investigations have used either a) a power-law 
extrapolation to low energies of the observed high-energy spectrum 
(>10 GeV), which is expected to be little affected by solar 
modulation (Meyer et al . , 1971; Schmidt, 1972) or b) a spectrum 
above a few hundred MeV calculated from the non- thermal -rad io- 
background data with a power-law extrapolation to lower energies 
(Burger and Swanenburg, 1971; Lezniak and Webber, 1971; Urch and 
Gleeson, 1972). These extrapolated interstellar spectra differ 
considerably, e.g. about a factor of 10 at 100 MeV. Nonetheless, 
by using different approximations to the transport equation and 
different diffusion coefficients, the authors have made the 
different interstellar spectra consistent with the data observed 
near Earth. In addition, the interstellar positron spectrum has 
been calculated by several authors, e.g. Ginzburg and Syrovatskii 
(1964), Hayakawa et al . (1964), and more recently by Ramaty and 

Lingenfelter (1966, 1968), Perola et al. (1968) , Beedle (1970), 
and Aral (1971). These calculated intensities also differ by 
factors of ~10. 

2) Diffusion Tensor 

Another problem in solar modulation studies is the evaluation 
of the interplanetary cosmic-ray diffusion coefficient, £. The 
diffusion of cosmic-ray particles results in part from pitch-angle 
scattering due to the irregular fluctuations of the interplanetary 
magnetic field. Several authors have derived equations relating 
the diffusion coefficient to the power spectrum of the interplanetary 
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magnetic field (Jokipii, 1966, 1967, 1971; Hasselmann and Wibberentz, 
1968; Roelof, 1968; Earl, 1972b). However, observations of the power 
spectrum are available for only a few limited time periods and 
generally do not cover a large enough frequency range to 
establish the rigidity dependence of the diffusion coefficient 
below a few hundred MV. Moreover, these measurements have all been 
made relatively near Earth and thus there are few observational 
data on the radial dependence of the diffusion coefficient . A 
further point of controversy is the question of the separability 
of k . Some authors (Burger and Tanaka, 1970; Burger, 1971; Burger 
and Swanenburg, 1971; L'Heureux et al . , 1972) argue that the 
diffusion coefficient must be a non-separable function of radius 
and rigidity in order to fit the cosmic-ray nuclei and electron 
data, whereas others have assumed a separable function in interpreting 
the data (Gleeson and Axford, 1968; Goldstein, Fisk and Ramaty, 1970; 
Fisk, 1971; Gleeson and Urch, 1971; Lezniak and Webber, 1971; 

Meyer et al., 1971; Urch and Gleeson, 1972; Garrard, 1973). 

3) Modulation Region 

Considerable speculation exists concerning the heliocentric 
radial distance to the "boundary" of the modulation region. Solar- 
flare studies have generally indicated a rather nearby boundary 
in the vicinity of 3-6 AU, whereas some studies of solar modulation 
have used much larger boundary distances, e.g. ~25 AU (Burger and 
Swanenburg, 1971; Fisk, 1971). 

4) Analytic Approximations to the Transport Equation 

Several investigators have used sufficiently simple approxi- 
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mations to the cosmic-ray transport equation such that analytical 
solutions are obtained (e.g. Gleeson and Axford, 1967, 1968; Meyer 
et al., 1971; Earl, 1972a jSchmidt, 1972). These analytic solutions 
diverge at low energies and hence the interstellar spectra and 
diffusion coefficients derived from these approximations 
necessarily differ. 

In this chapter we shall attempt to resolve some of these 
discrepancies. We shall first briefly review the basic physics of 
solar modulation and discuss our results of a numerical analysis of the 
cosmic-ray transport equation for electrons. We shall then discuss the 
results of a self-consistent study of solar modulation. The major 
elements of this study are; 

1) A new calculation of the possible range of interstellar electron 
spectra from the galactic non- thermal -radio -background data . From 
this range we shall discuss the absolute modulation of electrons 
above ~100 MeV. 

2) An interstellar positron spectrum from nuclear collisions 

in the interstellar medium . This spectrum is chosen by requiring 
agreement between electron and positron modulation above ~100 MeV. 

3) The rigidity dependence of the diffusion coefficient from ~10 MeV to 

» 

~10 GeV derived from comparisons of numerical solutions of the 
transport equation, using the interstellar spectra derived in 1) 
and 2), with the best available near-Earth data . At low energies 
(< 100 MeV) the results allow us to comment on the so-far unknown 


interstellar electron spectrum. 
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4) A comparison of these diffusion coefficients with those derived 
from the available power spectra of the interplanetary magnetic 
field . From this comparison we comment on the possible radial 
variation of k and make estimates of the size of the modulation 
region. 

5) Numerical solutions of the transport equation for the cosmic-ra 


nuclei us ine the electron-and positron-derived diffusion coefficients. 
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B . Background Physics 

The physics of the interplanetary medium responsible for solar 

modulation is well-established. The solar corona is dynamically unstable 

and expands outward from the Sun at supersonic velocity (Parker, 1963). 

The magnetic field of the Sun is frozen into the hot, fully ionized, 

highly conducting plasma and is swept outward into interplanetary space. 

The rotation of the Sun (with angular velocity Q) causes the magnetic 

lines of force to have, on the average, the shape of an Archimedes' 

spiral. Superimposed ' on this average shape are irregular fluctuations. 

Charged particles penetrating into the interplanetary medium are 

spiraling about this solar magnetic field and those with gyroradii 

comparable to the wavelength of the fluctuation undergo resonant 

pitch -angle scattering. This effect gives rise to a random walk of 

particles along the average field-line direction. In addition, 

particles are transported perpendicular to the field lines because the 

lines of force also execute a random walk. In general, this diffusive 

process is described by a diffusion tensor (Jokipii, 1971) which 

includes other effects such as curvature drift and the gradient of the 

average magnetic field. These curvature and gradient drifts are not 

expected to be significant for galactic energetic particles 

(Jokipii, 1970). In this thesis we shall ignore these terms of the 

diffusion tensor and treat only k and k , the diffusion coefficients 

II x 

parallel and perpendicular to the interplanetary magnetic field, 
respectively. 

The fluctuations which scatter the charged particles are being 
convected outward with the solar wind; hence, cosmic-ray particles are 
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convected away from the Sun producing a radial gradient which leads to 
diffusion in the opposite direction. In ! addition, the particles are 
scattering from field-line fluctuations which are, on the average, 
moving away from each other; hence, in these collisions the particles 
lose energy (adiabatic deceleration) . 

In a recent review, Jokipii (1971) describes how the inclusion 
of these effects leads to the following equation describing the 
propagation of galactic cosmic-rays in the interplanetary medium: 

1 i 

V-(VU) - (oru) _ v.(k;vU) =0 (VI-1) 


where U(r,T) is the number of particles per unit volume per unit energy 
at radial distance r with kinetic energy T (U = 4jrj/fJc where j is the 
intensity and pc is the particle velocity), V is the solar-wind velocity, 
a(T) is a parameter given by 


a(T) = 


5 ?,n T _ T + 2m 
3 f, n p T+m 


where p is the particle momentum, m is the particle rest energy, and 
K is the particle diffusion tensor. For electrons with energy above a 
few MeV, a is essentially unity. The three terms in equation VI-1 
represent, respectively, convection, adiabatic deceleration, and 
diffusion of charged particles in the interplanetary medium. 

Some direct information is available on the principal 
parameters, and k, which enter into the solution of equation VI-1. 
The solar-wind velocity is relatively constant from year to year and 
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has an average value near 400 km/sec (Gosling et al., 1971). Existing 
measurements of the power spectrum of the interplanetary magnetic field 
near 1 AU yield information on the rigidity dependence of the diffusion 
coefficient. In addition, further information is provided by solar- 
flare studies, since propagation of energetic flare particles is 
governed by a time -dependent equation of the same general form as 
equation VI-1. It is the long-term variation of these parameters, 
particularly k, which produces the time variation of the cosmic-ray 
intensities . 


/ 



55 


C . Review of Analytic Approximations tothe Transport Equation 

In order to treat equation VI -1 analytically or numerically 
it is necessary to construct a simplified model of physical conditions 
in the interplanetary medium. Several reviews have been published on 
the physics of the solar wind (Parker, 1965a ; Dessler, 1967; Parker, 
1969a) . In general, the presence of the solar wind is attributed to the 
radial expansion of the solar corona. The wind velocity becomes 
supersonic beyond a few tens of solar radii and roughly maintains a 

constant magnitude until the shock-termination boundary is reached. 

2 

At this point the stream-flow pressure, which falls as 1/r , is no 
longer able to sweep back the interstellar medium. This termination 
boundary is not necessarily the same as the cosmic -ray modulation 
boundary since the fluctuations of the interplanetary magnetic field, 
which scatter the particles, may be damped out in a shorter distance. 

On this basis we assume the solar-wind velocity, ^f, to be independent 
of heliocentric radius, r, and furthermore, for simplicity, we assume 
both V and the cosmic-ray density, U, are independent of angle about the 
Sun . 

Under these assumptions, the transport equation VI-1 becomes 


\ £ (A, - f & <<*0 - *2 h £> - 0 


(VI-2) 


where % is a scalar quantity, the radial diffusion coefficient, defined 
by (Jokipii, 1971) 

2.2 

K = K*. = K cos 9 + ic sin 6 (VI-3) 

rr II .L 
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where 0 is the angle between the radius vector from the Sun and the 

outward direction along the interplanetary magnetic field (~48° at 

1 AU), and «; and ^ are the parallel and perpendicular components of 

II x 

the diffusion tensor, respectively. 

Even with the assumptions made so far, no general, analytic 
solutions to equation VI -2 have been found. Several further approxi- 
mations to the equation have been made which lead to analytic solutions . 
We shall discuss the most important of these approximations as they 
apply to electrons. (Garrard (1973) has given a detailed discussion 
relevant to the cosmic-ray nuclei). 

1. The Diffusion-Convection (DC) Approximation 

If we neglect the adiabatic deceleration term in equation VI-2 

we obtain: 


a- (r 2 VU - r 2 K S2) . o 

ar * dr' 

In the absence of sources or sinks at the origin this equation may be 
written: 


VU = 



(VI -4) 


which is a statement of the balance between the outward current of 
particles due to convection and an inward current due to diffusion. The 
solution to the DC equation is: 


U(r,T) = U (co , T ) exp 



(VI -5a) 


If we assume there is a boundary at distance D beyond which V/^ is zero. 
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then 


U(r,T) = U(D,T) exp 


u 

■/ 


K ( r ' T) 


dr' 


= U(D,.T) e"t^ r,T ^ 


where the quantity ^ , defined by: 

D 


V 


♦■ (r,T) I tC ( r ' T) 


7^7 dr' 


(VI -5b) 


(VI-5c) 


(VI-6) 


is called the "modulation parameter" or simply the "modulation". It 
will be shown to be the determining parameter in the study of the solar 
modulation of electrons. It has also been found to be the important 
parameter in the discussion of the high-energy (> few GeV) solution 

t 

of the transport equation for nuclei (Garrard, 1973) . 

It is interesting to note that if the near-Earth and inter- 
stellar electron spectra are known, then the modulation parameter at 
1 AU, (1, T) , is determined in the DC approximation from equation VI-5c, 

i .e . , 

♦ d.T) = in (VI-7) 


If the radial and energy dependences of the diffusion coefficient are 
separable, i.e. |<(r,T) = K^( r ) then the energy dependence of k 

is determined from ^(1,T) (see equation VI-6). 

2. The Force-Field (FF) Approximation 

Gleeson and Axford (1967, 1968) have derived an approximate 
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solution in the case of small modulation by making use of the radial 
differential current density (or streaming), S, defined as: 

s = vu ‘ K §7 ‘ 3 |t (aTU) (VI ' 8) 

The second term on the right represents the contribution from diffusion. 
The remaining two terms represent the effective radial current due to 
the transformation between a frame of reference at rest with respect to 
the solar wind and the observer's reference frame (the Comp ton- Get ting 
effect). The first term represents the contribution due to convection. 
The origin of the third term can be visualized by imagining that we 
observe particles at a single energy T with a "directional" detector. 
Then if we point the detector toward the Sun, the velocity of the solar- 
wind frame effectively adds energy, +&T, to the particles we are 
observing. Thus in our frame we observe the rate of particles 
corresponding to the intensity at energy T-^T in the solar -wind frame. 

If we point the detector radially away from the Sun we observe the rate 
of particles corresponding to a different part of the solar-wind-frame 
spectrum, i.e. corresponding to the intensity at T+^T. Since these 
intensities are usually different, there is an effective radial current. 

Gleeson and Axford present arguments to show that S is 
negligible whenever VL/k « 1, where L is a length characteristic of 
the radial variation of the diffusion coefficient. If one assumes 
S = 0 and that «; is a separable function of radius and energy, then one 
obtains the so-called force-field solution (Gleeson and Axford, 1968) 

j( r »W) = i(D, W + $) 

2 2 .2 2 
W - m (W+$) 


- m 


(VI -9) 



59 


where j is the particle intensity, W is the total energy of the particle, 
m is its rest energy, and § is a spectral shift parameter which is 
determined from the diffusion coefficient. 

This simple equation (VI-9) has been used by several 
investigators in interpreting electron data. We shall show that the 
zero -streaming assumption breaks down at low energies 200 MeV) for 
certain forms of the diffusion coefficient and the unmodulated spectrum. 
In general, it is difficult to predict the range of applicability of 
the FF solution. In the following section we shall investigate 
some cases of interest in the light of results of a numerical 
solution of the full cosmic-ray transport equation. 

3. The Convection -Adiabatic Deceleration (CAD) Approximation 

If one assumes that the diffusion term (containing k ^) in 
equation VI -2 is small compared to the other terms, we have 


*2 I? <r2 U) - S tf (OTU) “ 0 (VI - 10) 

r 

Rygg and Earl (1971) solved this equation, assuming a to be constant, 
by 

„/ n^/O- - 1 <1 , „3/2a. 

U(r,T) = T J (r T ) 

For a = 1 (relativistic electrons), 

U(r,T) = T 2 ^ (r T 3/2 ) (VI-11) 

The function 3 is an arbitrary function to be determined by the 
boundary condition. If the boundary condition is U(r,Tg) = constant = 
U(D,Tq) at some boundary energy Tq (i.e. no modulation for T > T^), the 
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solution reduces to: 

U(r,T) = (T/T q ) 2 U(D,T Q ) (VX-12) 

or 

j(r,T) = AT 2 

where A is a constant. The corresponding solution for non-relativistic 
cosmic -ray nuclei is j = AT. 

It is interesting to note that in this model particles at the 
boundary with energy T < T^ do not propagate into 1 AU since k (T < T Q ) 
is assumed zero. However, the intensity of particles with energy 
T > Tq is the same at all radial distances since «; (T > T^) is assumed 
infinite. Therefore, in this model particles arrive near Earth with 
energy T < Tq only by being decelerated from higher energies. 
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D . General Results from Numerical Solutions to the Transport Equation 

A numerical solution to the transport equation VI-2 has 
advantages over the analytic approximations. For example, numerical 
solutions can be readily obtained for any specified radial and energy 
dependence of y , whereas the analytic approximations are often 
restricted to certain functional forms of the diffusion coefficient. 

In addition, the numerical solution can be obtained for all values of 
radius and energy of interest; these solutions can then be used to test 
the validity of the analytic approximations. 

We have constructed a numerical solution to equation VI-2 

based on the Crank-Nicholson technique outlined by Fisk (1971). In this 

technique the continuous radius-energy plane is replaced by a grid with 

maximum radial distance D and an energy range from T . to T . The 

° mxn max 

transport equation is expressed as a finite difference equation in terms 
of grid location. By specifying three boundary conditions, the resulting 
set of simultaneous equations can be solved for all the radius-energy 
grid points. 

The boundary conditions which we must specify are: 

1) r = D: We assume that beyond the boundary the modulation is 

negligible at all energies. Thus U(D,T) = galactic spectrum. For 
electrons and positrons we have information on these spectra from 
the non-thermal-radio-background data and galactic nuclear 
collisions, respectively. For nuclei one generally has to rely on 
power-law extrapolations of the high-energy data. 

2) r = 0: In order to eliminate source-like solutions at r = 0, we 

transform the equations so that the dependent variable is X =<frU 
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and require X -* 0 as r 0. This transformation implies that the 

solutions are not valid for small r (^0.2 AU typically). 

3) T = T : We assume that at sufficiently high energy, T , the 

max max 

modulation is negligible. Thus U(r,T ) = U(D,T ) for all r. 

max max 

In our work we assume T =10 GeV for electrons. 

max 

1 . Analysis of Analytic Approximations 
a. The CAD Approximation 

The solution of the CAD (convection-adiabatic deceleration) 

2 

approximation is j = AT . If we examine the spectra of Figure VI-1 we 

find that, in view of the error limits, the data for 1969-1971 might 

2 

be consistent with j a T over a limited energy range (~100 - 500 MeV). 

It has been suggested (Luhmann, 1971; Earl, 1972a) that such a segment 

is due to the validity of the CAD approximation (the low-energy turn-up 

is attributed to the dominance of diffusion below ~100 MeV). How- 

2 

ever, we shall now show that such a turn-over to j a T would not be due 

to the dominance of adiabatic deceleration and convection over diffusion. 

In fact, it will be noted later in the discussion of the DC solution 

(Section VI.D.l.b) that such a turn-over can also exist in that solution 

which totally ignores adiabatic deceleration. 

In discussing the CAD solution it is useful to define the 

3 3 

phase-space density, F (the number of particles per unit volume (d r d p) 
in the six-dimensional r, p phase space), by: 

F = j/p 2 

2 

Therefore, j o T is equivalent to F = constant for relativistic 
electrons. From the numerical solution of the transport equation over 
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the entire radius-energy plane it is easy to determine if F is constant 

over some portion of the plane and therefore to determine over what 

region the CAD approximation is valid. As an illustration, we show in 

Figures VI-2a and b the numerical solution of equation VI-3 using the 

interstellar spectrum and diffusion coefficient assumed by Meyer et al . 

(1971) in interpreting the 1970 data. This diffusion coefficient in 
2 

cm /sec is represented by: 


k(R) 





4.94 x 10 
4.94 x 10 


17 PR R > R 

17 pR c R <; R 


c 


c 


= 440 MV 


(VI -13) 


where pc is the particle velocity, R is the particle rigidity, and a 
constant radial dependence for k is assumed. Figure VI-2a shows 
electron intensity versus kinetic energy at 1 AU. As we shall point 
out shortly the numerical solution and the force-field approximation, 
used by Meyer et al., disagree at low energies and hence the assumed 
parameters dq* not lead to a good fit to the low-energy data. However, 
we use the solution at this point for purposes of illustration. In 
Figure VI-2b we show a plot of the contours of equal phase-space 
density in r,T space for the numerical solution. A region containing 
few contour lines would imply F constant and would possibly indicate 
that the CAD approximation is valid. In the figure adjacent contours 
are separated by a factor of 2 in F. (Note that the contour lines are 
horizontal in the large r region only because of the assumed boundary 
distance of 3 AU.) We observe only a slight spreading in the lines near 
500 MeV for r & 1.5 AU. This spreading reflects the turn-over in the 
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spectrum and is not necessarily an indication that diffusive effects 

2 

are small (it is an indication that j oc T over a limited energy range, 
however). Since no large area with F = constant is indicated, we 
conclude that the CAD approximation is inappropriate for the case 
considered. 

We expect the CAD approximation to be valid when diffusion can 
be neglected, i.e. when the diffusion coefficient is small or 
equivalently when the modulation parameter, * (see equation VI-6), is 
large. We shall show below (Section VI.D.2.a) that the numerical 
solutions of the transport equation at radius r using a given inter- 
stellar electron spectrum are primarily determined by ^(r,T) (defined by 
equation VI-6). Therefore, to determine the region of r,T space in which 
the CAD approximation is valid, we determine the minimum ^ which yields 
F constant. As an illustration, we show in Figures VI-2c and d (solid 
lines) the numerical solution for the same parameters as in Figures VI-2a 
and b, except that R c , the rigidity at which the diffusion coefficient 
changes form, is lowered from 440 to 100 MV. The effect of lowering 
is to lower the diffusion coefficient by roughly a factor of 4 at 

rigidities below 100 MV. It is evident from Figure VI-2c that the 

2 

intensity at 1 AU is roughly proportional to T at low energies 

2 

(compare with dashed j a T line), and in Figure VI-2d a large region 
where F is nearly constant does exist. Thus by sufficiently lowering 
the diffusion coefficient, we find, as expected, that the CAD approxima- 
tion (which ignores diffusion), is approximately valid, although the 
spectrum so obtained does not resemble the observed spectrum. The 
dashed lines in Figure VI-2d are curves along which ^(r,T) is constant. 
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We find, roughly, that % must be small enough such that ^ ~ 15 for the 

CAD approximation to apply. Using the nominal interstellar spectrum 

discussed in Section VI.E.l, we shall find that solutions at 1 AU 

consistent with the observations are obtained with ^ ~ 6 . In order 

for the CAD approximation to be valid near 1 AU the absolute 

interstellar electron intensity would be required to be a factor of 
-6 

— — ** 8000 (see equation VI-5c) larger than our calculated intensity, 
e 

2 

We therefore conclude that the interpretation of a possible j = AT 
segment in the 1969-1971 electron spectra is almost certainly not the 
manifestation of the validity of the CAD approximation. 

\ 

In the above analysis we have assumed that the F = constant 

region in Figure VI-2d is a result of the validity of the CAD 

approximation (for ^ ^ 15). We might ask the question: Is the 

F = constant region necessarily due to the dominance of adiabatic 

deceleration over diffusion? To clarify this question we show in 

Figure VI-2e the phase-space density contours for the diffusion - 

convection solution, i.e. the solution ignoring adiabatic deceleration, 

for the identical parameters as in Figure VI- 2d. Here we find no large 

F = constant region, which indicates that, indeed, adiabatic deceleration 

is responsible for the large blank area in Figure VI-2d. Note, however, 

the region of spreading in the DC phase- space density contours 

2 

(Figure VI-2e), which is a result of the fact that j = AT segments over 
a limited energy range can also be produced by approximations which 
ignore adiabatic deceleration. 

b. The FF and DC Approximations 

Meyer et al. (1971) and Schmidt (1972) have related modulated 
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spectra observed near Earth to the interstellar spectrum by using the 
force-field approximation. Some of the subtleties of this approximation 
have been pointed out by Garrard (1973). For example, the FF solution 
ignores adiabatic deceleration, although it does include the Compton- 
Getting effect. Furthermore, the conditions which produce zero 
streaming, required by the approximation, are not clear. 

In order to determine the region of applicability of the FF 
and the DC approximations, we show in Figure VI-3a a comparison of these 
solutions with the numerical solution (FN) of the full transport 
equation. For all three models we use the same diffusion coefficient, 
assumed independent of radius with a boundary at 3 AU and with the 

R > R = 300 MV 
c 

R R 

c 

This diffusion coefficient and the interstellar spectrum shown in the 
figure were used by Meyer et al. in interpreting their 1968 data using 
the force- field approximation. Below ~100 MeV the FF result diverges 
significantly (factor of ~10 too small) from the full numerical 
solution. Schmidt has fit the same data (shown in Figure Vl-la) using 
the FF approximation with a steeper interstellar spectrum and a some- 
what different diffusion coefficient. Since the force-field approxi- 
mation is inconsistent with the full numerical solution at low energies, 
the parameters used in deriving these approximate solutions are 
necessarily inconsistent with ones we derive from fits to the data 
based on the numerical solution. Note that the DC solution is a better 


rigidity dependence (cm /sec): 


K(R) = 


7.15 x 10 17 pR 


7.15 x 10 17 pR 
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approximation (within a factor of ~2) of the numerical solution over 

the whole energy range depicted than is the FF solution. 

It is difficult to predict under what circumstances the force- 

field solution is a reasonable approximation. We find that the 

numerical and force- field solutions are more consistent if the inter- 

-2.5 

stellar electron spectrum is flatter than T * at low energies, 

(below a few hundred MeV) . To illustrate this improvement we show in 
Figure VI-3b a similar comparison of solutions as in Figure VI-3a except 
that we have used the interstellar positron spectrum calculated by 
Ramaty and Lingenfelter (1968). This galactic spectrum flattens out 
gradually below ~1 GeV and eventually turns over below ~50 MeV. Both 
the DC and the FF solutions are within a factor of ~2 of the full 
numerical solution over most of the energy range from 10 MeV to 10 GeV. 
The region of validity of the FF solution probably also depends on the 
rigidity dependence of the diffusion coefficient (Urch and Gleeson, 
1972). 

We find that the DC solution is a fairly good approximation 
for a wide range of interstellar spectra and diffusion coefficients. 
Thus, if k is assumed to be a separable function of r and T, then, for a 
given interstellar spectrum and the spectra observed near 1 AU, the 
energy dependence of the appropriate diffusion coefficient can be 
estimated reasonably accurately by computing ty(l,T) from equation VI-7. 
The effect of adiabatic deceleration is to shift the solution in energy. 
For example, compare in Figure VI-3a the peak position in the numerical 
solution near 500 MeV with the peak in the DC solution. The observed 
shift indicates that for the assumed parameters the fractional energy 
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loss near 500 MeV is roughly 0.25. A more detailed discussion of the 
energy loss of the galactic electrons which penetrate to 1 AU is given 
in the next section. These energy losses mean that small adjustments 
in the diffusion coefficient estimated on the basis of the DC approxi- 
mation are necessary in order to yield a good fit to the data using 
the numerical solution of the full transport equation. 

2. General Remarks Concerning the Solution of the Transport Equation 
In this section we present some general results from a study 
of both the numerical and analytic solutions of the transport equation 
(VI-2). In what follows the diffusion coefficient is assumed 
to be a separable function of radius and rigidity, k ( r,R)=PK-^ (r)^^) • 
Several authors have argued on the basis of their modulation studies 
that k must be a non-separable function of radius and rigidity. 

(Burger and Tanaka, 1970; Burger, 1971; Burger and Swanenburg, 1971; 
L'Heureux et al., 1972). However, Gleeson and Urch (1972) have pointed 
out that these arguments are based on the assumption that the rigidity 
dependence of an assumed separable diffusion coefficient does not 
change form from year to year. No necessity for such a restriction has 
been suggested. We shall find that adequate fits to the observed 
cosmic-ray data can be made using separable diffusion coefficients with 
different rigidity dependences for different years. In this regard, we 
note that recent hysteresis studies of neutron monitor data have shown 
that the rigidity dependence of the modulation parameter changed 
abruptly several times during the last solar cycle. (Carmichael and 
Stoker, 1970; Carmichael and Katzman, 1971; Stoker and Carmichael, 1971; 
Kane, 1972). Thus we feel that there is, as yet, no compelling 
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observational evidence for either separability or non- separability . 

The assumption of separability provides us with a convenient framework 
in which to study the radial and rigidity dependences of the diffusion 
coefficient . 

a. Effects of the Radial Dependence of the Diffusion 
Coefficient 

There are few observational data on the radial dependence 
of k • Jokipii and Coleman (1968), from analysis of the inter- 
planetary magnetic field data of Mariner IV, find no drastic 
changes in the parallel diffusion coefficient, % , between 1 and 1.5 AU. 

'I 

The solar-flare studies of Lupton (1972) imply that a diffusion 
coefficient independent of radius, r, inside 1 AU is more consistent 
with the data than one which varies linearly with r. Sari (1972a), 

using power spectra from Pioneer 6 magnetic -field data, finds k (50 MeV) 

-2.7 ^ 

varies approximately as r ' between 0.82 and 1 AU. Observational 

data on the radial dependence of k are non-existent. We recall that k 

l 

depends on both and , whose radial dependences are uncertain, and 
11 J. 

on 6 (see equation VI-3), which is given by 

0(r) = tan' 1 (^) (VI-14) 

where Q is the angular velocity of the Sun. Thus the radial dependence 
of k is highly uncertain. 

We do not expect large variations in the calculated electron 
spectra if we change the radial dependence of since the DC solution, 
which depends only on the integrated effect of k from the Earth to the 
boundary, is a reasonable approximation to the transport equation. How- 
ever, it is useful to investigate just how much variation we do obtain 
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by varying the radial dependence of k, in order to see how well the 

parameter <j, determines the solution. For this purpose we use an 

illustrative set of radial functions shown in Figure VI-4a. In 

addition to the e r * and r° dependences used by many authors, we include 
3 

a 2 + r dependence modeled after the diffusion coefficient suggested 
by Ng and Gleeson (1971) to explain solar-flare observations, a 
[(r-3/2)^ + 1/4] function constructed by Garrard (1973) as an analytic 
representation of a possible scatter-free region near the Sun, a r ^ 
function suggested recently by Jokipii (1972) as a possible radial 
dependence beyond 1 AU, and a simple r^ dependence. For convenience the 
diffusion coefficient is assumed to be infinite beyond 3 AU. For the 
cases illustrated we use the single rigidity dependence of shown in 
Figure VI -4b. 

In Figure VI -5a we show the numerical solutions of the 
transport equation using an interstellar spectrum derived from the non- 
thermal- radio data (see Section VI.E.l) and each of the 6 radial 
dependences of k shown in Figure VI -4a. Also shown are the Caltech and 
Chicago data (Schmidt, 1972) for 1968. The absolute magnitudes of the 
diffusion coefficients are normalized such that 


Since 



= 1950 MV (appropriate for 1968) 


(VI -15) 


*(1,R) 


_J 

Pk 2 (R) 
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and (S^CR) i- s the same f° r each solution, we are assuming that each 
solution has the same ^-value at 1 AU. As we see from the figure, all 6 
curves are nearly identical, indicating that the solution for electrons 
is essentially independent of the radial dependence of «; . 

We have also investigated the effect of varying the boundary 
distance D while maintaining ^(1,T) constant (by adjusting the magnitude 
of k). In Figure VI-5b we show the solution for ic^( r ) = constant for 
boundary distances of 3, 5, 10, 15 and 30 AU. The slight increase in 
energy loss with boundary distance is evident from the small shift of 
the curves; however, the differences are small compared to the 
experimental uncertainties in the data. 

We conclude from these studies that ^ is indeed the determining 

I 

parameter in the study of the modulation of electrons, 
b. Effects of Adiabatic Deceleration 

In Section VI.D.l we demonstrated that the approximation 
ignoring energy loss allows a reasonable first estimate of the diffusion 
coefficient. In order to make refinements to ^ we need to include the 
effect of adiabatic deceleration. 

The rate of energy loss through adiabatic decleration is 

given by 

If - - | a(T) T (v-?) 

(Parker, 1965b). In the case of relativistic electrons (a = 1) and a 

— * 

constant radial solar-wind velocity, V, this equation reduces to 


1 dT 
T dt 


2V 

3r 


(VI-16) 
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Thus the particles lose more energy near the Sun than near the boundary. 
In addition, if the diffusion coefficient is large, particles diffuse 
in from the boundary very quickly and thus lose less energy than they 
would if the diffusion coefficient were small. Thus we expect the 
energy loss to depend on both the magnitude of the diffusion 
coefficient and the boundary distance, D: 

In Figure VI-6 we illustrate an energy-loss calculation. The 
unmodulated spectra, shown by solid lines, are of the form 

U(D,T) = A exp £ -50 (^n T/T Q ) 2 

where T^ takes on the values in MeV of 25, 50, 100, 200, 400, 800, 1600 
and 3200, We use the diffusion coefficient of equation VI-13 which is 
assumed independent of radius within a boundary distance of 3 AU. The 
envelope formed by the peaks of the unmodulated spectra corresponds to 
the galactic secondary positron spectrum calculated by Ramaty and 
Lingenfelter (1968). The corresponding spectra at 1 AU, representing 
numerical solutions of the transport equation, are shown as dotted 
lines in the figure. As a consequence of adiabatic deceleration, the 
peaks in the spectra at 1 AU are shifted in energy from the corresponding 
peaks in the unmodulated spectra. (Note that the shape of the near- 
Earth spectrum, as indicated by an envelope of the peaks of the dotted 
curves, is determined by the particular choice of the rigidity dependence 
of the diffusion coefficient.) 

For k independent of radius and energy (^440 MeV in this 
example) the fractional energy loss is independent of energy (see 
equation VI-16) . As we mentioned before, however, the fractional loss 
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does depend on the magnitude of k and the boundary distance D. Since 
the electron spectra calculated from the transport equation depend 
primarily on the modulation parameter , we show in Table VI-1 the 
values of the fractional energy loss (for Tq = 50 MeV) for several 
representative combinations of ^(1 ,Tq) and D, assuming the functional 
dependence of the diffusion coefficient given by equation VI-13. The 
increase in the fractional energy loss with increasing ^ for a given 
value of D is due to the inverse relationship between ^ and «:> i.e. 

^ (Larger implies smaller k implies larger energy loss.) 

For a given D is approximately proportional to k and the effect of 
a larger boundary (larger energy loss) is roughly offset by the effect 
of the correspondingly larger k (smaller energy loss) . From the table the 
fractional energy losses for electrons are less than roughly 0.5 for 
typical values of ^ and D. Above 440 MeV, where the diffusion 
coefficient is proportional to rigidity (in this example), we found 
that the fractional energy loss (for a given ^ and D) is smaller than 
that shown in the table. 

c. Discussion of the "Flat" Portion of Electron Spectra and 
Rough Estimates of the Modulation (^) 

We now discuss the characteristic "flat" region from 
approximately 100 MeV to 1 GeV of the observed electron spectra shown 
in Figure VI-1. We have shown that the simple diffusion-convection 
solution is a reasonable first approximation to the numerical solution 
of the full transport equation (VI-2) over the energy range 10 MeV - 
10 GeV. The force-field solution is a somewhat better approximation 
above a few hundred MeV but breaks down rather badly in some cases below 
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TABLE VI -1 

Fractional Energy Loss* for Representative Values of the 
Modulation Parameter, ^(ljT), and the Boundary Distance, D. 


D \ 

(AU) \ 

(1,T) _ 

1.25 

2.5 

5.0 

3 


.19 

.27 

.35 

7 


.25 

.35 

.47 

10 


.25 

.36 

.50 

15 


.25 

.39 

.52 


+ 

The values apply to the energy range (<440 MeV in this example) where 
the diffusion coefficient is independent of energy. 

* 1 • j ' 

V(D-l) 1 

For k = constant, = — 


K 
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~100 MeV. We shall use both approximations in the discussion. The 
main points of this discussion are: 

1) to show that the energy dependence of the observed spectrum is 
consistent with our understanding of 

a) the interstellar spectrum 

b) the diffusion 1 coefficient and 

c) the DC and FF approximations above ~100 MeV and 

2) to demonstrate that an estimate of the absolute modulation can be 
made from a knowledge of only the energy dependences (and not the 

. absolute magnitudes) of the diffusion coefficient and the inter- 
stellar spectrumy together with an observation of a relative peak 
in the observed spectrum near Earth. 

The latter point is interesting because the energy dependence of the 
interstellar electron . spectrum, for example, can be deduced more 
accurately from the non-thermal-radio-background data than the absolute 
interstellar electron intensity. (See Section VI.E.l.) 

We first examine the simple diffusion-convection solution given 
by equation VI-5c. Since the differential intensity is given by 
j = we can rewrite equation VI -5c (for r = 1 AU) as: 

j(l,T) = j (D, T) e"* (1 » T) (VI-17) 

where ^ is related to the diffusion coefficient y through equation VI-6. 
Measurements of the power spectrum of the interplanetary magnetic field 
are consistent with « a PR^ where R is the particle rigidity, pc the 
particle velocity (p = 1 for relativistic electrons), and b is rigidity 
dependent ranging from ~0.5 at low rigidities (~1 GV) to ~1.5 at high 
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rigidities (^10 GV) . (See Section VI. E. 4.) Using this form for the 
diffusion coefficient, equation VI-17 reads 

j(l, T) - j(D,T) exp [- 73 / (f (b)T b )] 


where Tj is defined by equation VI-15 and we have replaced R by T which is 

* 

valid in energy units for relativistic electrons . (F(b) is chosen to make 

the energy dependence of k continuous and is normalized such that f(l)=l. 
For. example, if k is represented by two joined power-law segments such 
that for T < T , b = ~, and for T > T., b = 1, then f (|) =\/F.) If we 
assume the interstellar spectrum j(D,T) a T ^ and if we approximate 
1 C OC R^ we find a relative maximum in the near-Earth spectrum at energy 
T given by: 1 i ‘ 


T 

m 


= (-^-) 
r (b) 


1 /b 



(VI-18) 


In terms of this condition reduces to 


to-V-b w-w 

Thus an estimate of the absolute modulation at the observed peak energy 

T is given simply by the ratio of the spectral index of the inter- 
m 

stellar spectrum to the exponent of the rigidity dependence of the 

diffusion coefficient. From the radio data (Section IV.E.l) we find 

that y ~ 1.8 below ~2 GeV; hence, we have a( 1,T ) ** 1.8 and, 3.6 for 

y m 

values of b of 1 and 1/2, respectively. These values of ^ correspond to 


R-*S 

ze 

(v c, 


where p is the electron momentum, 
z = 1), pc ^ T(MeV) ~ R(MV) . 


For relativistic electrons 
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absolute modulation factors, et } of approximately 6 and 36, respectively. 
From the observed spectra (Figure VI -1) we note that these values apply 
to an energy near 1 GeV. 

We can also apply this procedure to the positron spectra 
since, again, the energy dependence of the calculated interstellar 
positron spectrum is known more accurately than its absolute intensity. 
Near 1 GeV the interstellar spectra presented in Section VI. E. 2 
(Figure VI-11) have a spectral index of ~2.4. Thus, we have 
^ (1,T ) 2.4 and 4.8 for b = 1 and 1/2, respectively. 

It is interesting to note that the peak in the positron 
spectrum at 1 AU is expected to occur at a lower energy than the peak 
in the electron spectrum. From equation VI -18 we have 


T (e + ) 
m 

T m (e) 

m 


Y(e ) 

= 0.75 


1/b 



= 1 


= 0.56 b = 1/2 


(VI-20) 


where we have assumed y(e) = 1.8 and y(e ) = 2.4. Unfortunately, 
positron data of sufficient accuracy to observe this difference do not 
presently exist. Hopefully, future observations will confirm the 
prediction. Along these same lines, we note that if the energy 
dependence of the diffusion coefficient does not change significantly 
from year to year, the peak energy should move to higher energies 

- I ’ „ 

with increasing modulation as is observed. (See Figure VI-1.) 

Larger modulation implies larger ^ implies larger f) implies larger 


These estimates of the modulation parameters and peak-energy 



78 


locations are based on the DC approximation. Figure VI-3 indicates 

that for the choice of parameters considered the peak-energy locations 

are nearly the same for the FF and FN (full numerical) solutions and 

that both occur at a somewhat lower energy than that of the DC 

solution. Thus we might expect an improvement in the estimates of T 

m 

and ^(1,T ) by examining the FF approximation. 

The FF solution is given by equation VI-9. If we assume b = 1 
(i.e. k oc R^), the parameter § in equation VI-9 has the simple value: 

§ = J 71 (VI-21) 

(Gleeson and Axford, 1968) 

Equation VI-9 may be maximized to yield the peak energy T^: 

T, ’ 3 y C FF > b * l ) < VI - 22 > 

Note that this value of is just 2/3 of the DC estimate (equation VI-18) 
The corresponding ^ value is given by: , 

,|r(l,T n ) = |y [FF, b = 1] (VI-23) 

Equations VI-22 and VI-23 represent our best estimates of T^ 

and ^'(1,T ) for the case k a R* 3 ^ . The parameter § in equation VI-9 

has a complicated energy dependence for other rigidity dependences of 

K. Hence, it is difficult to evaluate T m and ^(l,T m ) in the FF 

approximation for b ^ 1. However, we mention that a comparison of the 

numerical and DC solutions of the transport equation (described below 

and shown in Figure VI- 7) indicates that T = 2/3 of the DC estimate 

m 

b = 1 / 2 

(equation VI-18 ) is also appropriate for «! « R 
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In the electron spectra shown in Figure VI-1 we find the best- 
resolved peak occurs in the observed 1970 spectrum. The spectrum 
begins to deviate from a power-law below ~4 GeV and a relative peak is 
observed in the region near 1 GeV. If we assume the rigidity dependence 
of tCj we can determine the interstellar electron intensity at 
~1 GeV as a function of its spectral index, i.e. of its energy de- 
pendence. From the non- thermal -radio-background data we estimate that 
y is approximately in the range 1.7 to 1^9 at energies below ~2 GeV 
(see Section VI.E.l). Assuming «; o: R'*’ we’ find from equation VI-23 
^,(1 GeV) ~ 2.55 - 2.85. Since (~1 GeV) ^ .006 electrons/ (m^ sec 

sr MeV), we estimate from equation VI-17, j (~1 GeV) ~ .077 - .104 

0D 

2 

electrons/ (m sec sr MeV). These values are consistent with our 
calculation of the galactic electron spectrum from the radio data (see 
Figure VI-9). We note that with an accurate determination of the 
rigidity dependence of e.g. from the power-spectra data, we could 
place more stringent limits on the interstellar electron intensity near 
1 GeV than those shown in the figure. (The limits indicated in 
Figure VI-9 result from uncertainties in^ the galactic parameters and do 
not reflect the analysis discussed here.) This particular refinement 
in the analysis will not be pursued in this thesis. 

We are now in a position to interpret the flat portion of the 

observed electron spectra at Earth during the years 1965-1971. From 

the power-spectra data (see Figure VI-19) we can infer roughly the 

rigidity dependence of the diffusion coefficient above a few hundred 

MV. As ah illustration, we idealize k(R) such that below a break-point 

1/2 1 

rigidity R , .^.0£ R , and above , k a R . If the values of 
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<n and R are such that T (equation VI-18) is greater than 
1 c m 

2 n 

, we expect a peak in the near-Earth spectrum at T^ = - -[ 

1 GeV in 1970). If R £ is too small, the diffusion coefficient 

becomes small enough that the modulation is very large. In this case a 

well-defined peak is obtained in the numerical solution, which in 

general is not observed in the data. For larger F the diffusion 

c 

coefficient is larger at low energies and the well-defined peak becomes 
a broad flat region as is observed. In this example a second "peak", 
corresponding to the b = 1/2 segment of the diffusion coefficient, 
would be obtained in the near-Earth solution of the transport equation. 

In terms of R^ this peak, T^^j can be calculated from equation VI-22 
and VI-18 (multiplying by the factor 2/3 on the right-hand side of 
equation VI-18): 

' 3T 2 

T l/2--8iT < VI ' 24 > 

c 

where is the b = 1 peak given approximately by equation VI-22. For 

R ~ T, ~ 1 GeV, we have T, ~ 3/8 T, * 375 MeV. 
cl If £■ 1 

In Figure VI- 7 we demonstrate the general features described 
in the example above. The unmodulated spectrum is derived from the non- 
thermal-radio -background data (Section VI.E.l). Both the full numerical 
solution and the DC approximate solution are shown for the case where 
R^ = 750 MV and f| = 1950 MV. From equation VI-22 we compute ~ 722 MeV 
and from equation VI -24 we obtain = 26 1 MeV. These values, which 

agree with the position of the two "peaks" in the FN solution, are shown 
in the figure. Note that these values are about 2/3 of the DC peak 
energies, which are also indicated in the figure. The excellent agree- 



81 


; » 


ment of the solution with the 1968 data indicates that the particular 
choice of parameters used in the calculation is reasonable. 

We feel that the above interpretation of the flat portion of 
the spectra forms a good, self-consistent picture with our present 
knowledge of the rigidity dependence of the diffusion coefficient and 
of the energy dependence of the interstellar electron spectrum. We 
also note that by knowing the energy dependences of both the diffusion 
coefficient and the interstellar electron spectrum and by observing a 
relative maximum in the near-Earth spectrum, we could, with the 
analysis presented here, determine the absolute interstellar electron 


intensity. 
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E . Quantitative Study of Solar Modulation 

In the previous section we have qualitatively discussed the 
parameters of solar modulation theory. In this section we present a 
quantitative study of these parameters using: 

1) the expected range of the interstellar electron spectrum, 

2) the calculated interstellar positron spectrum, 

3) a derivation of the modulation parameters and, hence, the diffusion 

. ! I 

coefficients, using the interstellar spectra of 1) and 2), 
numerical solutions of the transport equation, and the spectra 

i ■ 

measured pear Earth, 

4) a comparison of these diffusion coefficients with those derived 

j 

from the available power-spectra data, and 

5) numerical solutions of the transport equation for cosmic -ray protons 

and He nuclei using assumed interstellar spectra and the 

diffusion coefficients derived from the electron modulation study. 

1. Determination of the Interstellar Electron Spectrum 

from the Galactic Non-Thermal-Radio-Background Emission 

The galactic non-therma 1-radio-background data are generally 
ascribed to synchrotron emission from relativistic electrons spiraling 
in the galactic magnetic field. Several authors have examined the 
radio data and, under certain assumptions on the galactic parameters, 
have calculated the corresponding galactic electron spectrum 
(Anand et al., 1968a, b; Verma, 1968; Webber, 1968, Goldstein, Ramaty 
and Fisk, 1970, Burger, 1971). However, knowledge of the galactic 
parameters is incomplete and previous studies have not attempted to in- 
dicate the possible range of interstellar intensities which are consistent 

with the acceptable range of the parameters. We shall assume that cosmic- 
ray electrons are distributed uniformly in the galaxy and estimate a 
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reasonable range of interstellar electron spectra consistent with the 
non-thermal -radio-background data and consistent with our knowledge of 
the galactic parameters. With this range of spectra we shall estimate 
the expected limits on the cosmic-ray modulation parameter. 

In Figure VI-8 we show the non- thermal -radio data in the 

* 

galactic anticenter direction . The data above MHz are from the 
compilation by Webber (1968). Above ~40 MHz the data are primarily from 
high resolution surveys (~1° aperture) in which the galactic disk is well 
resolved. From ~ 5-40 MHz Webber has adjusted the available low (~30°) 
and medium (~10°) resolution data. He has used the high -resolution 
measurements of Blythe (1957) at 38 MHz to normalize the emissivities 
found in the lower -resolution studies to the standard anticenter 
direction. These adjustments amount to increasing the observed low- 
resolution intensities by ~10-30«g. Below 5 MHz we have plotted the 
recent data of Alexander et al. (1970) from their instrument on board 
the RAE-1 satellite. The data from this low-resolution instrument are 
representative of emission and absorption from a broad region on the 
order of ~100 degrees in angular extent. Since the disk of the galaxy 
is only ~0.5 kpc thick (compared to a radius of ~15 kpc) it may be 
inappropriate to regard these data as representative of emission and 
absorption in the disk of the galaxy. However, these authors find the 
intensity below ~ 5 MHz to be nearly isotropic. This finding may 
indicate that due to interstellar absorption the radio emission at these 

The radio data in the anticenter direction are used in this analysis 
because they are somewhat easier to interpret than the data available 
in other directions (e.g. towards the galactic center or the direction 
of minimum brightness). 
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frequencies is coming from a local region. Our assumed range of galactic 
parameters includes values for which the absorption arises within both lo- 
cal (^100 pc) and extensive. (~4 kpc) regions. However, even in the latter 
case the size of the region depends on the assumed galactic structure, such 
as the distance to absorbing cold clouds, etc. We shall interpret the low- 
frequency radio data as corresponding to ! emission and absorption in the 
disk of the galaxy. (Stephens (1971) has made an alternative interpretation 
in which emission from a spherical galactic halo is assumed to account for 
most of the low-frequency emission observed in the halo directions. Below 
~200 MeV the interstellar electron intensity he assumes is slightly lower 
than that derived in this analysis.) 1 

The general features of the radio spectrum in Figure VI-8 are: 

- 4 

1) a segment from ~10 to ~150 MHz which is proportional to \> ' with 

_ ~j _ g 

evidence for a steepening to v * above ~150 MHz and 

2) a segment from 0.4 to 1 MHz which is approximately proportional to 

1.6 

v • 

There is a smooth connection between these power-law segments in the 
intermediate frequency range of 1 to 10 MHz. Above about 10 MHz inter- 
stellar absorption becomes negligible (see Appendix B). If we assume 
the interstellar electron spectrum to be a power law in energy with 
spectral index y, the intensity of synchrotron radiation is: 

I ±1 

I(vj.B ,L) a C L B 2 v" a V ~ 10 MHz (VI-25) 

X I 

(Ginzburg and Syrovatskii, 1964), where 
u is the frequency 

a = - is the resulting power-law index of the radio emission, 

L is the line-of-sight emission length. 
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B is the perpendicular component of the magnetic field, 

X 

C is the constant used in defining the spectrum of cosmic -ray 
electrons, i.e., 

N(W) = CW" Y ' 

where W is the total ienergy and N(W) represents the number of electrons 
3 

per cm per unit energy. Thus the power-law segment of the radio data 

-.4 

above 10 MHz proportional to v ’ implies that y 1.8. We estimate, by 
constructing straight-line segments through the data points of 
Figure VI-8, that a range of slopes from -0.35 to -0.45 is consistent 
with the data. Hence] y is roughly in the range 1.7 to 1.9. 

At low frequencies free-free absorption by interstellar 
electrons and the Razin suppression of emission (Razin, 1960) must be 
considered. The Razin effect is due to the ambient electron density 
which causes the index of refraction to be greater than 1. We include 
both the Razin effect and free-free absorption in our calculations. We 
find, as did Ramaty (1971), that the Razin suppression represents in 
general a much smaller effect than free-free absorption. In Appendix B 
we show that in the case of large optical depth, and if we ignore the 
Razin effect, 

±±X 

2 2-a 

I( V ,B,L) a C B V v « 10 MHz (VI-26) 

X 

Thus, at low frequencies the intensity is independent of the total line- 
of-sight emission distance, L. To match the observed ^ dependence in 
the low-frequency range we again have a = 0.4 or y ~ 1.8, consistent 
with the result from the high-frequency data. If we use the 15$ error 
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bars of the data of Alexander et al., as they suggest for the relative 
error between points, the range of observed slopes below 1 MHz is 
roughly 1.3 - 1.65, which implies a range of 1.7 to 2.4 for the spectral 
index of electrons, y* The addition of Razin absorption will alter 
this picture slightly;' however, this effect can be roughly compensated 
for by raising the interstellar temperature (see Appendix B) . 

At higher frequencies C5 150 MHz) the radio spectrum 

appears to be steepening to v ^ ° r which implies a cosmic-ray 

electron spectral index y ** 2.4 to 2.6. Most of the emission at 150 MHz 

comes from a region of the electron spectrum near 2 GeV (Appendix B). 

Thus we see that an interstellar electron spectrum consistent with the 

non-thermal-radio-background has a power-law index of ~1.8 below ~2 GeV 

which steepens to ~2.5 at higher energies. We have chosen a set of 

reasonable galactic parameters (see Table B-l), assumed a simple 

galactic model, and calculated the resulting electron spectrum necessary 

to fit the radio data, assuming the power-law indices of 1.8 and 2.5 

mentioned above. In Figure VI-9 we display this nominal spectrum 

(which is essentially identical to one derived by Goldstein, Ramaty 

and Fisk, 1970) along with two spectra which result from a study of the 

•k 

reasonable range of galactic parameters . This calculation is described 
in detail in Appendix B. We note that at high energies there is 
roughly a factor of 4 between the bracketing lower and upper spectra. 

The electron spectra are plotted between 70 MeV and 5 GeV. This energy 
range has been chosen so that there is less than 25^ contribution to the 
radio emission at the minimum and maximum frequencies, 0.4 and 600 MHz, 
from electrons outside this energy range. This definition is consistent 
with the absolute accuracy of the low-frequency radio data (±25$) • 
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(For convenience, the electron spectra are shown as connected power-law 
segments. The coefficients and spectral indices for each segment are 
given in Table B-2.) 

In Figure VI-10 we show the calculated radio background which 
results from each of these electron spectra. All three curves are 
consistent with the radio data. In each calculation the different 
galactic parameters shown in Table B-3 were used. 

The galactic electron spectra shown in Figure VI-9 will be 
used in Section VI. E. 3 to derive the cosmic-ray modulation parameter 
from ~100 MeV to ~5 GeV. 

2. Interstellar Positron Spectra from Galactic Nuclear Collisions 

Due to the limiting lower frequency of the radio data, the 
interstellar electron i spectra derived in Section VI.E.l are restricted 
to energies above ~ 100 MeV. In order to complement the electron 
modulation studies at lower energies, we use the near-Earth Caltech 
positron data and the calculated interstellar positron spectrum. 

Since thereis no evidence to indicate the existence of anti- 
matter stars or galaxies, it is commonly assumed that there are no 
sources of primary, directly accelerated positrons. Above ~10 MeV 
nuclear interactions in the interstellar medium are believed to be the 
only source of cosmic-ray positrons. (Radioactive decays may 
contribute at lower energies.) Based on this collison-source 
mechanism, several authors have attempted to calculate the local inter- 
stellar positron spectrum. It is usually assumed that particles 
propagate by diffusing through the turbulent galactic magnetic field 
and that a steady-state exists between particle production and loss. 
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This process is described by: 

Q(r,W) = (^7 U) - ^-(k5 U) 


(VI -2 7) 


(Jokipii and Meyer, 1968) 

where Q(r,W) is the rate of positron production as a function of 

position, r, and total energy, W, U(r,W) is the positron density, and 

«;(r,W) is the galactic cosmic-ray diffusion coefficient. The rate of 
dW 

energy loss, — , is determined by losses from ionization, bremmstrahlung, 

synchrotron radiation, and Compton collisions with starlight and the 

universal blackbody radiation. A reasonable approximation for the 

dW 

energy dependence of tt is 


dW 2 

—■ = a + bW + dV/ 
dt 


(Beedle , 1970) 

where a is constant (ionization loss), b is constant (bW = bremmstrahlung 

2 2 
T) 13 

loss) and d a — +w , , where — and w , are the energy densities of 
bit ph bit ph 

the magnetic field and the photons, respectively. 

There have 'been two basic approaches to solving equation Vi-27. 
We briefly discuss each of these. 

1) Leakage- lifetime approximation 

In this method the diffusion terms and boundary conditions 
are replaced by a leakage-loss term U/ T , where x is the "lifetime" of 
a particle before it escapes from the confinement volume. The production 
spectrum, magnetic field, hydrogen density, and photon density are 
considered independent of galactic position within the confinement 
region (disk or disk + halo). In the energy region most affected by 
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solar modulation (~L GeV) the spectra calculated under these 

assumptions are essentially only dependent on the source function and 

x = pc x , the amount of material traversed by the particles during the 

time t • Among the recent calculations based on this method are those 

of Perola et al . (1968), Ramaty and Lingenfelter ("R&L") (1968), and 

— 2 

Arai (1971). The spectra of R&L for x = 4 g/cm and Arai for 
— 2 

x = 3 and 5 g/cm are shown in Figure VI-11. These authors have used 
different positron production spectra which accounts for their 
different calculated intensities. For comparison we also show in 
Figure VI-11 the measured spectra for 1965-66 and 1968. 

2) Diffusion model 

Several authors have pointed out that the leakage-lifetime 

approximation may not be physically justified for electrons. (Shen, 

1967; Jokipii and Meyer, 1968; Beedle, 1970.) For example, the term 

U/t implies that all particles are assumed to have a constant 

probability of escape. However, in reaching the boundary an electron 

may lose much of its energy (particularly high-energy electrons since 

2 

their loss rate is a W ). Hence, setting the loss rate to U(W)/ T 
where W is the particle's initial energy is not correct. In addition, 
the assumption of a production spectrum independent of galactic 
position may not be reasonable and several investigators have assumed 
different distribution functions for the matter density in the galaxy 
in making their calculations (Ginzburg and Syrovatskii, 1964; Shen, 
1967; Beedle, 1970). 

Beedle has solved equation VI-27 using an ellipsoid of 
revolution for the distribution of matter in the galaxy. He also 
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assumed the energy-loss rate (equation VI-28) to be independent of 

position. In Figure VI-11 we show his calculation of the local inter- 

3 

stellar positron spectrum using the parameters p = 1 atom/cm , 

? 

B 3 29 2 

— + w = 1 ev/cm and k = 10 cm /sec. Of these parameters the 

calculated spectrum below ~1 GeV is most sensitive to the diffusion 

30 2 

coefficient, being ~50$ lower at 10 MeV for % = 10 cm /sec. 

The positron intensity calculated by Beedle is roughly a factor 
of 10 larger than that derived by Ramaty and Lingenfelter . Below ~1 GeV 

the difference is primarily due to the different models used in the 

29 2 

calculations. As noted by Beedle, his spectrum with % = 10 cm /sec 
is almost identical to the x = « (Le. T = oo) disk -model spectrum of 
Perola et al. (1968). (He used their positron production spectrum.) 

Thus, in his model the particles we observe locally are not being lost 
through boundary escape. It is also interesting to note that all the 
calculated spectra of Figure VI-11 have roughly the same energy 
dependence . 

We shall not attempt to choose between the various calculated 

spectra of Figure VI-11 on the basis of the validity of the models, the 

production spectra, etc. used in the calculations. Rather, we shall 

determine positron modulation parameters using both the bracketing high 

and low spectra of Figure VI-11. By requiring these results to be 

consistent with those of the electron modulation study, we shall 

approximately determine the galactic positron intensity. 

3. Derivation of Modulation Parameters and Implications for the 
Low-Energy Interstellar Electron Spectrum 

We have presented in Sections VI.E.l and VI. E. 2 our calculation 
of the expected range of interstellar electron spectra and the inter- 
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stellar positron spectra calculated by several authors. In addition, 
we have shown the data from the Caltech magnet spectrometer and from 
other experimenters covering the period 1965-71 (see Chapter V and 
Figures VI-1 and 11). We now derive the modulation parameter at 
Earth, 


ird,T) 


D 

f Vd 

A * (r ’ 


Vdr_ 
T) 


for each of the years for which we have presented data. The discussion 
divides itself conveniently into two sections: a) the electron inter- 

stellar spectra and near-Earth data are used to derive the 
modulation parameters and diffusion coefficients above ~100 MeV and 
b) the positron interstellar spectra and near-Earth data are used 
primarily for the discussion of the modulation below 100 MeV. 
a. Comparisons of Measured and Calculated Electron Spectra 

Following the discussion in Section VI. D, we assume the 
diffusion-convection model is a reasonable first approximation to the 


r KP.T) . 


transport equation and calculate ^(1,T) ,= ^n f f' h *~rf ] using the 

J > 1/ J 


near- 


Earth spectra and the range of galactic spectra of Figure VI-9. In 
Figure VI -12 a-e we show ^ as a function of energy for the periods 
1965-66, 1968, 1969, 1970, and 1971. Hie points are derived from the 
data of Figure VI-1 and are shown by circles which are filled for the 
high and low galactic spectra (connected by dashed lines) and open for 
the nominal spectrum. (Note that in Figure VI-12e (1971) some of the 
points are upper limits and hence the errors in the data points extend 
considerably beyond the dashed lines.) 
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Since the high and low interstellar spectra differ by roughly 
a factor of 4 above a few hundred MeV, the resulting uncertainty in 
At* is 



At lower energies (~300 MeV), where the galactic electron spectrum is 
more uncertain, this error increases. 

These ^-values are derived from the DC approximation . To 
refine these values we adopt the following procedure: we first 

determine the rigidity dependence of the diffusion coefficient by 
drawing power-law segments through the "nominal" ^-points of Figure VI-12 
(open circles). If we assume a constant radial dependence for % out to 
a boundary distance D; we have from equation VI -6 

! i 

K (r - R >-$MO (VI - 29) 

(We arbitrarily choose D = 10 AU in what follows; note from Figure VI-4b 
that the calculated spectra are practically independent of D.) Using 
this diffusion coefficient and the nominal interstellar electron 
spectrum, we determine the numerical solution of the transport equation. 
Since the DC approximate solution, upon which the diffusion coefficient 
is based, differs slightly from the numerical solution, the spectrum 
generated in this way does not represent the best fit to the data. 

Better agreement is achieved by adjusting slightly (in the manner 
described below) the rigidity dependence of the diffusion coefficient. 

To facilitate computation we use the following model for the diffusion 
coefficient (consistent with the power-spectra data - see Section VI. E. 4): 
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r 


P ( r > k 2 ( r > 


K(r,R) = 


C pR 


R > R, 


cp V rr 2 r l < r < r 2 


(VI -30) 


cpV^R 2 R £ r l 


where C is a constant. By adjusting the values of C, R^, and R^, we 
have calculated the electron spectra for the periods 1965-66, 1968, 
1969, 1970, and 1971. i These spectra are 1 shown in Figure VI-13. The 
values of C, R^, R 2 , and 71 (defined by equation VI-15) for each epoch 
are given in Table VI-2 and the corresponding modulation parameters 
are shown as solid lines in Figure VI-12. For the present, below 
~100 MeV, we have used an extrapolation of the nominal interstellar 
electron spectrum and a diffusion coefficient which is arbitrarily 
defined to be independent of energy. Therefore, we do not discuss 
quantitatively the electron modulation below ~100 MeV. In the next 
section we shall use the positron spectra in a discussion of the 
modulation at low energies. 

Because of statistical errors of the data, the modulation 
parameters used in deriving numerical solutions in agreement with 
the data (solid lines in Figure VI-12) are not the only ones possible. 
As examples, in each of Figures VI-12a, b, and c we show limiting 
modulation parameter curves (dotted lines) from which acceptable fits 
to the data were derived using the same nominal interstellar electron 
spectrum. We note that these ^'s we have used in achieving consistency 
with the data differ only slightly from the ones calculated from the 
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TABLE VI -2 

Diffusion Coefficient Parameters -- Electron Modulation Study 



Epoch 

T 

(MV) 

R 1 

(MV) 

R 2 

(MV) 

C*(xl0 18 ) 

2 

(cm sec /MV) 

June -July 

1965-66 

(averaged) 

1350 

64 

900 

4.006 

June-October 

1968 

1950 

160 

750 

2.773 

June -July 

1969 

2400 

182 

1000 

2.253 

June -July 

1 

1970 

3300 

312 ! : 

1100 

1.639 

July 

19 71 

2700 

480 

480 

2.003 

* 

These values 

are based on 

a solar 

wind velocity V = 

400 km/sec and 


on a diffusion coefficient assumed independent of radius with a 
boundary D = 10 AU. 
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simple diffusion-convection approximation (open circles in Figure VI-12), 
confirming the discussion in Section VI.D.l.b. 

The dotted modulation parameter curves of Figures VI-12a, b, 
and c represent an uncertainty in ^ due to the uncertainty of the measured 
spectra. In discussions of the absolute magnitude of k we will need to 
consider the larger limits on y (dashed lines in Figure VI-12) which 
result from the uncertainty in our knowledge of the interstellar electron 
spectrum. 

We have thus determined the rigidity (energy) dependence of 
the modulation parameter for electrons above ~100 MeV for each of the 
periods 1965-66, 1968^ 1969, 1970, and 1971. Under the assumption k 
separable in radius and rigidity, the rigidity dependence of ^ 
determines that of the diffusion coefficient. We now use the positron 
data to discuss the modulation parameters and diffusion coefficients 
at lower energies. 

b. Comparisons of Measured and Calculated Positron Spectra 

In Chapter V and Figure VI-11 we have shown the positron data 
for the same years used in the electron study as well as the calculated 
interstellar positron spectra. As we pointed out in Section VI. E. 2 
the various calculated interstellar positron intensities differ by a factor 
of ~10. We now demonstrate that the calculated spectrum of Ramaty and 
Lingenfelter (lower curve in Figure VI-11) yields modulation parameters 
consistent with the electron results while the spectrum of Beedle 
(upper curve) does not. 

For the period 1965-66 we used the electron data of Fanselow 
et al . (1969) to derive the modulation parameters in the preceding sec- 
tion. Their instrument also yielded the only available positron fluxes 
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above a few hundred MeV in the 1965-66 time period. We have chosen to use 
only the data of Fanselow et al. in that time period, in order to eliminate 
possible intercalibration problems between various instruments and thus 
preserve, as much as possible, consistency between positron and electron 
data. In Figure VI-14a we show a plot of the DC modulation parameter ^ 
based on the positron data and the calculated galactic positron spectra 
of both Beedle (circles) and ML (squares). We have also plotted by lines 
the i|; from the electron study (Figure VI-12a) appropriate for the period 
1965-66 when the positron data were taken. The dashed lines 
correspond to the estimated limits on ^ derived from the uncertainty in 
our knowledge of the ! interstellar electron spectrum. This band is 
consistent with the tyf"data points" of R&L; on the other hand, the points 
of Beedle are in considerable disagreement. In Figures VI -14b, c, d, 
and e we show similar plots covering the years 1968-1971. In the cases 
where the modulation parameters from the electron and the positron 
studies overlap, e.gi ! 1970 and 1971, the ^-points of R&L agree better 
with the electron values than do those of Beedle. Thus we believe 
that the interstellar positron intensity 'is roughly that calculated by 
Ramaty and Lingenfelter . However, we do not imply that the leakage- 
lifetime model on which his calculation is based is necessarily correct. 
Since all the calculated positron spectra have roughly the same energy 
dependence, we only infer from the modulation studies the approximate 
magnitude of the interstellar intensity. In the galactic nuclear 
collisions calculation, this magnitude depends on the assumed positron 
production spectrum as well as on details of the specific model. 

There is other evidence that the positron intensity derived by 
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Beedle is not valid. For example, the intensity of electrons he cal- 
culates from galactic nuclear collisions is so high that it passes through 
the high-energy electron data observed near Earth. Thus, no "primary" 
source of negatrons is needed if we assume his calculation is correct. 
However, several authors have concluded from observed positron fractions 
that there exists a dominant primary source of cosmic-ray negatrons 
(Beuermann et al . , 1969; Fanselow et al., 1969). ■ 

It has been suggested, however, that these observed positron 
fractions may not be representative of the interstellar positron 
fractions since electrons and positrons may be modulated differently, 
e.g. because the energy-loss effect depends on spectral shape (Beedle, 
1970). To investigate this possibility^ we used the numerical solution 
of the transport equation to calculate the positron fraction at the 
boundary and at 1 AU for two different cases. In both cases we used the 
galactic positron spectrum of R&L and, above 100 MeV, the nominal 
galactic electron spectrum (Figure VI-9). In model 1 we extrapolated 
the nominal galactic electron spectrum to low energies and used a diffu- 
sion coefficient appropriate for 1968 (see eq. VI-30 and Table VI-2). The 
resulting positron fractions (labeled MODEL 1) are shown in Figure VI-15. 
We have included the Chicago (1965-66) and Caltech (1969) observations 
for comparison. At low energies these measured points fall above the 
calculated curves. We can achieve better agreement by using an inter- 
stellar electron spectrum which turns over below 100 MeV (solid line in 
Figure VI-16b) and by modifying the diffusion coefficient such that 
below 60 MV, k OL 1/R. (We shall discuss this behavior at low energies 
shortly.) The resulting positron fractions are labeled M0DEL-2 in 
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Figure VI-15. In general, the net effect is a shift of the I AU curve 
to lower energy with respect to the boundary curve. The shape of the 
positron-fraction curve is roughly preserved, however. The shifts are 
not large and the differences are smaller than the statistical errors 
of the data. Thus, we conclude that the 'positron fractions measured 
near Earth are also indicative of conditions in interstellar space and 
hence that a primary negatron source is required. Hence, we believe that 
the galactic secondary electron and positron intensities derived by 
Beedle are too large. 1 ' 

The low-energy values of the positron modulation parameters 
shown in Figures VI-14b, c, d and e indicate that the modulation is 
decreasing at low energies. This decrease is most evident in the plots 
for 1968 and 1969 (Figures VI-14b and c). For example, the solid line 
in Figure VI-14b represents a modulation parameter consistent with both 
the electron and positron studies. The segment below 60 MeV is propor- 
tional to 1/T (or, equivalently, 1/R). In Figure VI-16a we show the 
numerical solution of' the transport equation at 1 AU using the diffusion 
coefficient derived from this modulation parameter and the interstellar 
positron spectrum of Ramaty and Lingenfelter (1968), together with the 
measured spectrum. The calculated spectrum is in excellent agreement with 
the low-energy data. In Figure VI-16b we show the numerical solution 
for the electron flux using the same diffusion coefficient. In 
order to achieve agreement with the data, the interstellar electron 
spectrum was turned over below 100 MeV as shown by the upper solid line 
in the figure. 

The indicated energy dependence of the interstellar electron 
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spectrum in Figure VI-16b is not the only one possible. Hie dotted 
lines in Figure VI-14b show the possible range of the positron modulation 
parameter at low energies in 1968. This range was determined by 
considering both the uncertainty in the observed 1968 positron spectrum 
(Caltech) and the uncertainty in the magnitude of the interstellar 
positron intensity. The latter uncertainty was derived by assuming that 
the range of the modulation parameter derived from the electron 
intensity in 1965-66 (dashed lines in Figures VI-12a and VI-14a) also 
applied to the positron intensity in the region of overlap measured 
during the same period. (Fanselow et al'., 1969). We have used the 
dotted modulation parameter curves of Figure VI-14b to "demodulate" the 
1968 low-energy electron data of L'Heureux et al. (1972) (see 
Figure VI-1), assuming the diffusion-convection approximation, i.e. 
j (D, T) = j(l,T)e'* f ’ 1 . In Figure VI-17 we indicate the resulting range 
of interstellar electron spectra below ~50 MeV as a shaded region 
bounded by dotted lines. Above ~70 MeV we show as a shaded band the 
range of spectra consistent with the analysis of the non -thermal -radio - 
background data (Section VI.E.l and Appendix B) . For comparison we also 
show (solid line) the galactic electron spectrum of Figure VI-16b. We 
conclude from Figure VI-17 that the electron spectrum must flatten 
below ~100 MeV if the positron and electron modulation studies are to be 
consistent. 

These general features, i.e. the turn-over in the galactic 
electron spectrum and \^0t — , which are based on the 1968 positron modula- 
tion study, are also supported by the 1969 results. In fact, the lowest- 
energy positron data point at ~14 MeV is almost the same in 1968 and 
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1969, implying that nearly equal modulation was observed during the two 
years at low energies. The higher-energy positron data for 1969 are 
statistically not as accurate as the 1968 data and therefore a detailed 
analysis is not warranted. 

It is interesting to note that the solar-flare proton studies 
of Lupton (1972) are consistent with a 1/R dependence of at 

low rigidities. In the 1-10 MeV (43-137^ MV) region he finds k is 
roughly independent of kinetic energy . Since k Ct p ^(R); K = constant 
implies ic,(R) a For non-relativistic protons this is equivalent to 

W,(R) 

We note that Lupton' s solar-flare studies determine the 

magnitude of between the Sun and the Earth, whereas the modulation 

studies yield information on the diffusion coefficient beyond 1 AU. 

For the June 7, 1969 event Lupton (1972) derives a radial diffusion 

20 2 

coefficient of ~1 - 3.5 x 10 cm /sec for the 1-10 MeV protons. If we 

assume ^ independent of radius with a boundary at 12 AU (see next 

section) and use the positron modulation parameter from Figure VI-14b 

19 2 

(solid line), we derive k m 7.5 x 10 cm /sec for the low-energy 
protons. Thus the solar-flare result is a factor of ~3 larger than 

this estimate from the positron modulation study. This disagreement may 

imply that a) the boundary is at a larger distance than 12 AU or 
b) k is larger inside 1 AU than beyond. For example, if we assume 
K a ~ outside 1 AU ? then a boundary distance of ~9 AU yields a modulation- 
derived k consistent with the solar-flare result. 

We note that solar-flare studies have generally used boundary 
distances of ~3 - 6 AU in order to fit the observed exponential 
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decay. However, recent evidence (Marshall and Stone, 1972) indicates 
that during the time interval of the flare observations an equilibrium 
condition may not be reached and a larger boundary is also consistent 
with the data. 

In summary, the important conclusions from the positron 
modulation studies are a) the interstellar positron spectrum of R&L 
provides consistency between electron and positron modulation studies 
above ~100 MeV and b), below ~100 MeV the diffusion coefficient must 
increase and the interstellar electron system must flatten considerably 
to maintain consistency between positron and electron modulation. 

4. Relation of Cosmic-Ray Diffusion Coefficient to Power Spectra of 
the Interplanetary Magnetic Field and Implications for the Radial 
Dependence of the Diffusion Coefficient 

The diffusion coefficient derived from the electron and 

positron modulation studies is an "average" k for the entire modulation 

region, i.e. the modulation parameter essentially determines 

On the other hand, measurements of the power spectrum of the inter- 
planetary magnetic field, made near 1 AU, determine the local diffusion 
coefficient. Therefore, a quantitative comparison of these diffusion 
coefficients can provide information on the radial dependence of k, 
including estimates of the size of the modulation region. 

The basic theory relating the cosmic-ray diffusion coefficient 
to the magnetic -field power spectrum is described in the review paper of 
Jokipii (1971). Two methods of calculating the parallel diffusion 

coefficient are given, which yield the same result only if P a v \ 

i 

where P is the spectral density of the perpendicular fluctuations and 
l 
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v is the frequency. The different methods result from the use of 

different approximations to the basic Fokker-Planck equation describing 

the evolution of the particle distribution function. Recently, Earl 

(1972)) has developed an improved approach based on eigenfunctions of the 

operator which describes pitch-angle scattering. The evaluation of the 

lowest-order eigenfunction leads to a precise expression for the 

parallel diffusion coefficient. This method is easily applied when the 

power spectrum can be represented by a power law in frequency with one 

index. For typical values of the index (-.5 to - 2) Earl finds that 

the second method discussed by Jokipii, the perturbation method (1971; 

see also: Jokipii, 1966; Hasselmann and Wibberentz, 1968), yields 

results within 10$ of his calculation. The other method (Jokipii, 1968) 

gives results that differ considerably from those of Earl. In this 

discussion we use the perturbation method of Jokipii rather 

than Earl's method for two reasons: 1) the calculation is in 

terms of an integral of the power spectrum over frequency and hence it 

is more easily applied to various functional forms of P (v) and 2) 

l 

Jokipii' s method allows a determination of the appropriate rigidity 
range of the diffusion coefficient corresponding to the frequency range 
of the observed power spectrum. 

If the magnetic-field fluctuations are approximated as one- 
dimensional waves propagating along the field direction, we can represent 
Jokipii 's integral equation for the parallel diffusion coefficient as: 



(VI-31) 
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where R is particle rigidity, V is the solar-wind velocity, B is the 

magnetic field strength, and p=fJ /p=coscp where p c is the component of 

II . II 

the particle velocity, pc, along the direction of the field and cp is the 

pitch angle- Note that P ( v ) refers to a frequency spectrum defined for 

1 

negative as well as positive frequencies. The published power spectra 

are defined for positive frequencies only and hence must be multiplied 

by a factor of 1/2 before insertion in equation VI-31 . 

Power spectra during the relevant time period have been 

published by Jokipii and Coleman (1968), Sari and Ness (1969), Bercovitch 

(1971), and Quenby and Sear (1971). In Figure VI-18 we show the data 

—6 “ 2 

from these authors in the range 10 - 10 Hz for the indicated time 

periods. 

Sari (19 72b) has noted that, in evaluating « , the power spectrum 

II 

observed by the spacecraft must be converted to the power spectrum 

"observed" by a particle spiraling along the average magnetic field 

line. The power spectra of Figure VI-18 represent power at wavelengths 

which are frozen into the field and convected past the spacecraft in the 

radial direction. Particles, however, are scattered by the power at 

wavelengths along the field direction. It can be shown that if 

P a v q then the corrected power spectrum is: 

„ corrected , _,.q-l „ observed 

P = (cos 0) P 

l l 

where 0 is the angle between the average field direction and the Earth- 
Sun line (~48° at 1 AU). For q = 1.5, typical at high frequencies, the 
correction at 1 AU is about 16$. (Note that the above correction differs 
by a factor of cos 0 from that derived by Sari (1972b).) 
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We approximate the observed power spectra by power-law segments, 
apply the correction mentioned above, and use equation VI-31 to compute 

K for each of the power spectra (except the high-frequency data of 

!! 

Sari and Ness in which q = 2 where equation VI-31 breaks down). In each 

case we have assumed V = 400 km/sec and B = 5y (ly = 10 gauss). The 

results are displayed in Figure VI-19. 

We note that for a given rigidity, R, contributions to ^ (R) 

II 

in equation VI-31 come from the power spectrum at all frequencies 

VB 

greater than v m £ n = 2rtR* ^ or ® = and V = 400 km/sec the numerical 
relation between v . in Hz and particle rigidity in GV is 


= 0-955 x 10 
v min R 


(VI-32) 


Since a given power spectrum only extends up to some maximum frequency, 
-4 -2 

v (*^10 - 10 Hz), it follows from equation VI-32 that for 

max ^ 

R ^ R . = 0-955 x — , — none of the contribution to «; is derived from 

mm -j ! ii 

max II 

the measured power density of the field fluctuations. In calculating 
K from equation VI-31 it is necessary to extrapolate the measured 

H i 

power spectra of Figure VI-18 to higher frequencies. The diffusion 

coefficients shown in Figure VI-19 are plotted to a lower-limit 

rigidity, R such that less than 50^ of the contribution to k is 
low || 

from the extrapolated portion of the power spectrum. 

These diffusion coefficients can be characterized by a para- 


K a R 


meter b where 
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In general, we find that b ^ 1.5 - 2 near 100 GV and is slowly decreasing 

to perhaps .5 in the region near 1 GV. 

We now compare these diffusion coefficients with those 

derived from the electron modulation study. In that study we used the 

radial diffusion coefficient, k , which depends on both k and k and 

rr II x 

is given in terms of these quantities by equation VI-3. By substituting 
6 = 48° (i.e. the 1-AU value) in the equation we have 


K rr = 0.45 k + 0.55k 

it II 


(VI-33 


At present there is no consensus on the value of y . It has 

x . 

been suggested that perpendicular diffusion is dominated by the random 
walk of the field lines (Jokipii, 1966; Jokipii and Parker, 1969), as 

v 

measured by the power at zero frequency: 


r 


k ( r ) 

x 


B x 

(v-0, 

B X 


R <; 1 GV 


» 


1 GV 


(VI -34a) 


(VI -34b) 


Jokipii and Parker (1969) have noted that P (v) must have zero slope at 

x 

low frequencies and have estimated k by making a low-frequency extra- 

x 

polation of the observed power spectrum of Jokipii and Coleman (1968). 

However, the power at zero frequency is difficult to measure and such 

extrapolations represent only estimates of the perpendicular diffusion 

coefficient. Furthermore, these extrapolations may represent a large 

overestimate of k for the following reason. A large part of the observed 
x 

power at low frequencies may be due to tangential discontinuities 


being swept past the spacecraft. Such discontinuities may be visualized 
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as representing an interplanetary medium ' composed of many relatively 

disordered flux tubes 1 of plasma. Within a given flux tube the field 

has roughly the same average direction. It is not clear whether the 

contribution to the magnetic -field power spectrum from such 

discontinuities should be removed before; extrapolating to zero 

frequency in estimating k . 

■ l 

Because of the systematic nature of the uncertainty of k we 

A 

shall consider two limiting estimates: 1) we shall extrapolate the 

power spectra to zero frequency and use equation VI- 34 to estimate k , 

A 

and 2) we shall also consider the case k ~ 0 (i.e. ^ « k )• 

A A II 

As an example, we compare the radial diffusion coefficient 
from the 1968 modulation study with that derived from the power spectrum 
of Quenby and Sear (1971). In Figure VI-20 we show the ic rr 
derived from the power-spectra data for the period 12/68-3/69 using 

equation VI-33. Curves 1 and 2 correspond to the results for k =0 and 

21 2 X 
K = 4 x 10 cm /sec, respectively. The latter estimate is the zero- 

A 

frequency extrapolation result using equation VI-34a. (The high- 

rigidity result for % given by equation VI-34b is only a factor of 2 

A 

larger than the low-rigidity value. Since « k at high rigidities, 

A II 

the use of equation VI-34a for the entire rigidity range is reasonable.) 
The error bars on the two curves indicate the 2 CT uncertainty in the 
observed power spectrum assigned by Quenby and Sear. Since the 
quantity determined from the cosmic -ray modulation is an 

integral, ^ (r, R) = 

•Jr 

tion-derived k at 1 AU depends on its assumed radial dependence 


K~(r < ^b ', t * 16 actua ^ roagnitude of the modula- 
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including the boundary distance, D. Assuming a constant radial 
dependence with a boundary at 12 AU we obtain from the 1968 modulation 
study the diffusion coefficient shown as the solid line in Figure VI-20. 
The 3 representative error bars indicate the approximate uncertainty 
derived from the limiting modulation parameter curves of Figure VI-12b. 
Although there is only a limited region of overlap, the rigidity 
dependences of the diffusion coefficients derived from the power-spectra 
method and from the modulation study are consistent. 

We can place limits on the possible value of the boundary 
distance, D, by requiring that the magnitude of the modulation- 
derived if agree with that from the power spectra study. As an 
illustration, we show in Figure VI-21 a comparison at 1 GV 
of these diffusion coefficients as a function of boundary distance 
assuming if independent of radius. Hie two power-spectra estimates of % 
are shown as horizontal bands, corresponding to the 2 c uncertainty of 
the data. The boundary dependence of the modulation-derived % is: 


|C(D, 1 GV) 


V» (D-l) 
i|r (1 AU, 1 GV) 


and the band in this case results from the uncertainty in our knowledge 

of the interstellar electron spectrum. Hie crosshatched areas represent 

the intersections of the bands. We find that if if is negligible, bound- 

X 

ary distances of 6-15 AU are required for consistency between the two 

diffusion coefficients, under the assumption if independent of radius. 

21 2 

If If is 4 x 10 cm /sec (the value inferred from the zero- frequency 


x 

extrapolation of the Quenby and Sear power spectrum), we obtain the 
boundary range 11-25 AU. 
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If we assume a different radial dependence for k, the integral 

definition of the modulation parameter (equation VI -6) still determines 

the magnitude of the diffusion coefficient at 1 AU as a function of 

boundary distance. By requiring this magnitude to be consistent with 

that derived from the power spectrum we can calculate the limits on D for 

any specified radial dependence of k* As a simple example, we consider 

the case k « k and assume k is a separable function of radius and 
i. II 

rigidity with the radial dependence K^( r ) a r”.. We calculate the 
minimum and maximum D for different values of the index n. At the 
comparison rigidity of 1 GV in 1968 it can be shown from the definition 


of ,i ( (equation VI-6) that the functional forms of D . and D in AU 
™ min max 


for the case considered are given by: 


D ^- l + 5 -° f 1 -") 


D mx ■ 1 + 14 - 2 O-n) (VI-35b) 

In Figure VI-22 we show plots of these limiting boundary distances 

as a function of the index n. The horizontal bar at n = 0 

indicates the 6-15 AU range we obtain for k independent of r. If the 

1-n 

c 

index n has a value n^ such that £ 0, we cannot obtain consistency 

between the diffusion coefficients derived from the magnetic-field power 
spectrum and from the electron modulation study for any value of D. 

From equation VI-35a we obtain 

n > 1.2 
c — 
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21 2 

for the case k « k • For k; = 4 x 10 cm /sec the condition is 

x 11 x 1.1 
n > 1.1. Thus, if ir- is assumed to increase with r faster than ~r ' , 
c — 

there is not enough calculated modulation of electrons beyond 1 AU to 

agree with the observed modulation. 

Recently, Jokipii (1972) has calculated the radial dependence 

of k for two types of fluctuations: 1) Alfven waves and 2) frozen-in 

irregularities. Beyond about 1 AU he finds for Alfven waves, k Qt r°, 

and for frozen-in fluctuations, a 1/r. From equation VI-35 (or 

Figure VI-22) we find that a 1/r dependence would imply a boundary 

range of ~3.3 - 5.5 AU (4.6 - 7.0 AU for «; = 4 x 10 ^ cm^/sec). The 

x 

r behavior gives the 6-15 AU range we derived above. 

5. Comparisons of Measured and Calculated Spectra of Cosmic-Ray 
Protons and He Nuclei 

In Section VI. E. 3 we discussed the numerical solutions of the 
transport equation for electrons and positrons for the periods 
June-July 1965-66 (averaged), June-October 1968, June-July 1969, 
June-July 1970, and July 1971. Electron spectra consistent with the 
data were calculated using the nominal interstellar electron spectrum 
(Section VI.E.l) and the diffusion coefficients described by 
equation VI-30 using the values of the parameters listed in Table VI-2. 

The transport equation for cosmic-ray nuclei is the same as 
that for electrons (equation VI-1). It follows that we should be able 
to use the electron-derived diffusion coefficients in deriving 
numerical solutions appropriate for the nuclei. In this section we 
present such solutions under the following restrictions and assumptions 
1) rigidity dependences of the diffusion coefficient derived 

from the electron modulation parameters based on the nominal 
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interstellar electron spectrum 

2) diffusion coefficient independent of radius with boundary 
at 12 AU. (Calculated spectra for boundaries in the 
range 6-25 AU, derived in Section VI. E. 4, are identical 
above a few hundred MeV/nucleon and differ by less than 
20$ above ^0 MeV/nucleon . ) 

3) constant solar-wind velocity of 400 km/sec 

4) interstellar spectra of nuclei given by 

j = A(W - m/4) 2-65 p/(m^ sec sr MeV/nucleon) 

where W is the total energy per nucleon, m is the nucleon 

9 7 

rest energy and A is 1.07 x 10 and 7.67 x 10 for protons 
and He nuclei, respectively. (See Garrard (1973) for , 
discussion of this particular form of the interstellar spectra.) 
Cosmic-ray proton and He-nuclei data similar to those 
compiled by Garrard (1973) for the periods 1965-1970 are shown in 
Figure VI-23 a-d and Figure VI-24 a-d, respectively. In each figure 
we show one or two calculated spectra using diffusion coefficients 
described by equation VI -30 with the parameters listed in Table VI-3. 
Each calculated spectrum is marked with a number corresponding to an 
entry in the table. (Garrard used parameters similar to those shown 
in Table VI-3. His calculated spectra are slightly steeper at low 
energies, however, since he used boundary distances of 2.7 and 6.1 AU.) 
A complete discussion of the relation between the calculated and 
observed spectra is given by Garrard (1973). We include a short 
description for each of the epochs listed in Table VI-3. 



Ill 



* 

These values are based on a diffusion coefficient assumed independent 
of radius with a boundary at 12 AU. 
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1965-66 (Figures VI-23a and VI-24a): 

The calculated curve is based on the same parameters as used in 
the 1965-66 electron study. The curve is slightly above the low-energy 
proton measurements and slightly below the corresponding He-nuclei data 
points. We regard the fits as adequate for this "two-year" epoch. 

1968 p and 19 6 7- 68a (Figures VI-21b and VI-24b): 

The same parameters are used as in the 1968 electron study. 

The fit is good for both the proton and He-nuclei spectra. 

1969 (Figures VI-23c and VI-24c): 

We include two curves. Curve 3 is derived using parameters 
consistent with the electron study. Both the proton and He-nuclei 
measurements fall below the curve. A better fit is obtained by increasing 
by ~7$ (curve 4). This change is justified since the 1969 nuclei data 
were taken in August-September whereas the electron data were collected 
in June- July. Since may have changed by as much as 40^5 from summer 
1969 to summer 1970 (see Table VI- 2), a 7$ change over a two-month period 
is reasonable. (The turn-up in the observed proton spectrum below ~40 MeV, 
which is not reproduced in the calculated spectra, may be due to sblar 
emission (Garrard, 1973).) 

19 70 (Figures VI -23d and VI-24d) 

The calculated He-nuclei curve is slightly below the data 
points but the curve does fall within the error bars. The observed 
proton spectrum is much flatter than the calculated curve. Since this 
period is near solar maximum, the flattening may result from a combina- 
tion of depressed galactic fluxes and possibly enhanced solar emission. 

In general we regard the fits for the nuclei spectra for 1965- 
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1970 as adequate. We note that we have achieved reasonably good agreement 
despite the many restrictions imposed, i.e. the use of the ^'s derived 
from the electron data using only the nominal interstellar electron 
spectrum, the particular interstellar spectra of nuclei assumed, and 
the restriction k independent of r. We have not, however, ruled out 
the possibility of different interstellar nuclei spectra or more 
complicated radial dependences of k. On the other hand it is not 
necessary to invoke them. In addition, the argument (e.g. Burger and 
Swanenburg, 1971) that a diffusion coefficient which is non-separable 
in its rigidity and radial dependences is necessary to fit the electron 
and nuclei data is not seen to be true. However, we cannot rule out a 
non-separable diffusion coefficient. 

In summary, we find that the good agreement between the cal- 
culated and measured nuclei spectra indicate that 

a) the interstellar proton and He-nuclei spectra used are 
reasonable (although it must be remembered that at low 
energies (< few hundred MeV) the near-Earth spectra are 
relatively insensitive to the interstellar spectra) and 

b) the diffusion coefficients derived from the electron 
modulation studies are appropriate for the nuclei 

as well . 

Thus the nuclei form the final element in our consistent picture of the 
solar modulation of cosmic rays. 
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VII . SUMMARY 

In this thesis we have derived the expected range of the 
interstellar spectra of positrons and electrons (e + + e ) and have 
discussed the mechanism of solar modulation of cosmic rays. 

We have based our studies on cosmic-ray positron and 
electron spectra, measured by us and by other investigators. The 
observations covered an energy range of ~10 MeV to ~10 GeV and the 
time period 1965-1971. The studies presented here have led to the 
following conclusions: 

1) Analytic Approximations to the Cosmic-Ray Transport 

Equation 

We have used numerical solutions of the full transport 
equation describing cosmic-ray propagation in the interplanetary 
medium to discuss the validity of several analytic approximations 
to the equation. We have found that: i> 

a) In order for the convection-adiabatic deceleration 
approximation to be valid the interstellar intensity 
of electrons is required to be a factor of ~8000 
greater than that inferred from the analysis of the 
galactic non-thermal-radio-background data. 

b) The force- field approximation is inadequate at 
low energies. The diffusion coefficients and 
interstellar spectra derived by Meyer et al. (1971) 
and Schmidt (1972) using this approximation are 
inconsistent with our conclusions for energies 
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below '"'100-200 MeV. ; 

c) The diffusion-convection approximation yields a 
reasonable first-order solution of the transport 
equation for both electrons and positrons. Energy 
loss by adiabatic deceleration (& 5C$) leads to 
a shift in energy of' the numerical solution at 1 AU 
from the DC approximate solution, but, on the whole, 
the spectral shape is preserved. In the DC approxi- 
mation the logarithm of the ratio of the inter- 
stellar cosmic-ray intensity to the near-Earth 

r^Vdr 

intensity is the modulation parameter ^,(1 , R) = Jx R) 
Thus, if the diffusion coefficient, 1C, is assumed 
to be a separable function of radius, r, and rigidity, 
R, the DC approximation may be used to estimate 
the rigidity dependence of K from a knowledge of the 
near-Earth and interstellar electron (or positron) 
spectra. 

In addition, we have used the diffusion-convection and 
force-field approximations, together with numerical solutions of the 
full transport equation^ to discuss the expected behavior of the electron 
spectrum at 1 AU at energies above ~100 MeV. Assuming a nominal 
galactic electron spectrum, we found that the flat portion of the 
near-Earth electron spectrum from ~100 MeV to 1 GeV can be attributed 
to a change in the rigidity dependence of the diffusion coefficient 
near 1 GV. We have also shown that a knowledge of only the energy 
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dependences of the interstellar spectrum ^and the interplanetary cosmic- 
ray diffusion coefficient is sufficient to estimate the absolute solar 
modulation of electrons (or positrons) at the energy of a relative 
maximum in the near-Earth spectrum. Further analysis may lead to 
improved estimates of the absolute interstellar positron and electron 
intensities. 

2) Interstellar Electron and Positron Spectra 

We have made a new derivation of the approximate range of 
the interstellar electron spectrum at energies between ~100 MeV and 
~5 GeV from the non- thermal -radio -background data. Uncertainties in 
our knowledge of the galactic parameters used in the analysis lead to an 
uncertainty of about a factor of 4 in the electron intensity above ~300 
MeV and to larger uncertainities at lower energies. 

We have discussed several interstellar positron spectra 
calculated by other investigators for the energy range 10 MeV - 10 GeV. 
Since the calculated absolute intensities differ considerably, we 
determined the appropriate spectrum by requiring consistency between the 
electron and positron modulation studies at energies above ~100 MeV. We 
have used the most consistent interstellar positron spectrum (i.e. that 

proposed by Ramaty and Lingenfelter (1.968)) and the Caltech positron data 

• i. 

I 

to study the modulation of both positrons and electrons at low energies. 
From this study we concluded that the interstellar cosmic-ray electron 
spectrum must flatten considerably below ~100 MeV. Our solar modulation 

i ' 

studies indicate that the ratio of positrons to electrons in interstellar 
space is nearly the same as that at 1 AU. We thus conclude that the low 
values of the positron fraction measured near 1 AU imply that cosmic-ray 
electrons with energies above ~10 MeV have a predominantly primary origin 

i 
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3) Diffusion Coefficient and Size of the Modulation Region 

From a comparison of the interstellar and near-Earth 
electron (and positron) spectra we calculated the modulation parameters 
for the periods 1965-66, 1968, 1969, 1970, and 1971. These parameters 
were used to derive the approximate rigidity dependences of the 
diffusion coefficients for these periods. These rigidity dependences 
were compared with those calculated from measurements of the power 

spectrum of the interplanetary magnetic field. In the limited rigidity 
range where the comparison is possible, these rigidity dependences 
were consistent. For rigidities below ~60 MV we derived diffusion 
coefficients which increased with decreasing rigidity. This increase 
at low rigidities is consistent with the rigidity dependence inferred 
from the solar-flare proton studies of Lupton (1972) . A comparison 
of the magnitudes of the diffusion coefficients derived from the 
solar-flare and the modulation studies indicates that at low energies 
K may be larger inside 1 AU than beyond or that the distance to the 
boundary of the modulation region may be relatively large (~30 AU) . 

We have also derived limits on the possible radial 
dependence of (C by requiring that the magnitude of the modulation- 
derived diffusion coefficient be consistent with that derived from 
the power-spectra study. Assuming K(r) a r°, we found that n < 1.1. 

For K independent of radius, consistency between the magnitudes of the 
diffusion coefficients requires the boundary of the solar modulation 
region to be in the range 6-25 AU. 

We have also applied the diffusion coefficients derived 
from these electron modulation studies to the cosmic-ray nuclei. 
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Assuming a particular form for the interstellar spectra of protons and 
the nuclei, we have calculated spectra of these particles at 1 AU 
which are consistent with the observations. Since our complete 
analysis was done assuming a diffusion coefficient which is a separable 
function of radius and rigidity, and consistency with the measured 
spectra was achieved, we found (as have Gleeson and Urch (1972)) no 
necessity to invoke the non-separable diffusion coefficients proposed 
by some authors. 
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APPENDIX A 

Details of Data Analysis 
1 • Selection Criteria for Data Analysis 

The selection criteria for "acceptable events" are essentially 
the same for both detector configurations MOD-1 and MOD-2. The criteria 
have been explained in detail by Rice (1970). A brief review of the 
procedure will be given as well as a description of the differences 
between MOD-1 and MOD-2 selection criteria, 
a. Spark Chamber Performance 

Initially, ‘the trajectory of a particle in each spark chamber 
is determined by making a least-squares fit of the measured spark 
locations to a straight line. In some events either no spark or a 
spurious spark is registered in one or more modules within a spark 
chamber. In these cases the module is ignored in the least-squares 
fit. If more than one plane in a chamber malfunctions, the event is 
categorized as a "multi-error" event and rejected from analysis. 
"Perfect" events are those in which all 8 planes determine the 
trajectory in the two chambers, and "one-error" events are those in 
which an error is detected in one plane in either or both of the 
chambers . 

The trajectories of the "perfect" and "one-error" events are 
subjected to further tests to determine their acceptability. The 
average deviation of the measured spark locations must be within 
1.25 mm of the best-fit straight line or the event is rejected. In 
addition, the extrapolated trajectories in each chamber must fall within 
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the acceptance cone of the detector. 

In connection with the above criteria, multiple-particle 

events can be recognized by the multiple-spark-indicator (MSI) bit 

(see Chapter II). In analyzing the 1968 data it was necessary to ignore 

this bit because some of the modules developed persistent spurious 

sparks at the edge away from the pick-up coil. As explained by Rice, 

this condition did not significantly affect the data. The frequency 

of occurrence of these spurious edge sparks was subsequently reduced 

and in later years it was necessary to ignore the MSI bit in no more 

than one module per flight. The MSI feature is somewhat more important 

2 

in the analysis of MOD-2 data since the 2 g/cm of material above the 

v 

upper spark chamber (gas Cerenkov counter) is a possible source of 
contaminating particles. The contamination due to this effect is 
discussed in Appendix A.4.c. 

b. Trajectory- Consistency Check 

This simple test, which utilizes the symmetry of the detector 
and the magnetic field, determines whether the calculated trajectories 
in the two spark chambers are consistent with the bending expected in 
the field for the computed deflection angle. In Figure A-l we show a 
projected particle trajectory assuming no scattering and an idealized, 

uniform magnetic field that is completely confined to the gap. Outside 

l 

the gap the trajectories are straight line segments whereas inside the 
field region of the magnet the path is an arc of a circle. From 
simple geometry it can be shown that for the idealized path shown 



(A-l) 
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This equation (which is also valid for a field with symmetrical 
fringing above and below the magnet gap) is a necessary and sufficient 
condition for cd to be joined smoothly to the straight-line segments 
ca and db (see Figure A-l for definition of symbols). The angles 
and ^ are calculated from the trajectories in the spark chambers. 
However, the idealized values of and k ^ are not determined because 
of multiple scattering and the intrinsic angular resolution of the 
detector (see Appendix A. 2). The expected angular distribution for 
A due to these effects has been calculated (Rice, 1970), and the result 
for the standard deviation of a is 

■>*- W +4 ° ffi 2 (a - 2 > 

where cjd and 9 (the deflection angle) are defined in Figure A-l. The 
uncertainty in 0, CT g, depends on the angular resolution of the detector 
and is derived in Appendix A. 2 (equation A-ll). The uncertainty in to, 
CT ^, results from the uncertainty in the spark locations of the modules 
defining the line segment ab . Rice (1970) derived ~ .0014 radians. 
Using this value for and equation A-ll, equation A-2 reads: 


CT r: 

A 


\A .170) 2 + ( .004) 2 


,0680) 2 + (.004) 2 



MOD-1 
MOD- 2 


(A- 3a) 


(A-3b) 


A possible further contribution to due to non-uniformities in the 
magnetic field was found to be negligible. 

In the case of M0D-2 observations we adopt a selection 
criterion such that events with a a more than approximately 2^ away 
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from zero are rejected. Using a Gaussian approximation for the 

distribution with the a of equation A-3b, we would obtain for the 

A , ■ 

selection criterion: 

' A ^ A c 

Ac 2 = (.1360) 2 + ( .008) 2 

However, at low energies (large 0) where multiple scattering is 
important the Gaussian approximation is not very accurate . For example. 
Figure A-3b shows the probability for electrons to scatter through 
angle greater than cp versus momentum x cp. (This distribution function 
is calculated in Appendix A. 2) . At the "2<J" level (ordinate = 0.025) 
there is considerable disagreement in the actual computed distribution 
(curve 1) and the Gaussian approximation (curve 2) used in deriving 
in equation A-2. Based on these considerations we adopt as a 
criterion on A for acceptance of an event 

A £ A c 

Ac 2 = ( . 1550) 2 + ( .008) 2 M0D-2 (A-4a) 

Thus at high energies (0 0), A .008 ** 2Q and we reject ~5 $ of 

C Zj | 

valid MOD-2 events. At low energies (large 0), ^ .1550. From 

Figure A-3b at p A = .155p0 = 1.37 (equation II-2b implies p0 = 8.85) 
c I 

we find that we are rejecting about 7$ of the valid events. 

The selection criterion for MOD-1 events is based on similar 
considerations (Rice, 1970): 
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A s A c 

Ac 2 = ( . 320) 2 + ( .008) 2 MOD-1 (A-4b) 

However, the scattering angle distributions shown in Figures A-3a and b 
are slightly different from the distributions calculated by Rice. At 
low energies p A ^ .32p0 ** 1.15 (equation ll-2b implies p0 = 3.55). 

Thus we are rejecting about 11$ of the valid MOD-1 events at low 
energies. The high-energy rejection is ~5$. These limits and standard 
deviations apply equally well to the cosmic-ray nuclei since the 
scattering term is negligible for these events. Thus, the criteria 
introduce essentially no bias according to particle species or rigidity. 
On the other hand, the criteria are such that there is very little 
probability that an event will be accepted which includes spurious 
sparks in the trajectory determination. 

In Figure A- 2 we show the A distributions for the analyzable 
(perfect + one-error) events as measured at the Caltech Synchrotron for 
positrons of 85 and 790 MeV energy using the MOD-1 detector. The smooth 
curves represent the calculated Gaussian distributions based on the 
standard deviation given by equation A-3a. The shaded areas show the 
events which are rejected because of the /^-criterion (equation A-4b). 
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2 . Rigidity Resolution 

The resolution of the MOD-1 detector system has been described 

in detail (Rice, 1970). Much of that discussion is relevant to the 

present description. However, some refinements in the calculations 

have been made and therefore a general discussion will be presented. 

The ability; of the detector to measure the rigidity of a 

particle is principally affected by 1) multiple scattering within the 

chambers or magnet gap and 2) the intrinsic angular resolution. Multiple 

scattering of the electrons adds a random angular deviation to the 

true deflection angle. Most of the scattering occurs at the wire 

planes and aluminized mylar covers adjacent to the magnet gap. A rough 

calculation of this effect was made previously by Rice. We have made a 

more refined calculation which uses a better approximation to the true 

mass distribution of the wires and which also includes scattering in the 

gas of the chambers and magnet gap. In Figures A-3a and b we show 

(curve 1) the distribution of projected scattering angle cp in the y-z 

plane (see Figure A-l for definition of this plane) calculated for 

electrons of momentum p according to the theory of Moliere (Galbraith 

dN 

and Williams, 1964). Both the differential distribution — , which is 
normalized by dividing by the momentum p, and the integral angular 
distribution N (>p) are shown plotted vs . pep. As plotted, curve 1 can 
be used for all electron momenta above a few MeV/c. We also show in 

Figure A-3a and b (curve 2) a Gaussian distribution with o = .60 MeV/c 

pep 

radians. This distribution will be used to approximate the true 
scattering-angle distribution. In particular, note that the integral 
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distributions for curves 1 and 2 are equal at pep = .cr = .60 MeV/c 

P9 

radians (Figure A-3b). The long large-angle scattering tail, which is 
not well reproduced in the Gaussian approximation, is accounted for in 
the trajectory-consistency checking (see Appendix A.l.b). 

If ° is the standard deviation of the distribution as 

pep 


plotted, then 

a = 

cp. p 


(A-5) 


We thus have: 


.60 




radians 


(A- 6) 


By substituting for p the values from equation II-2 we obtain 



.179 

MOD-1 

(A- 7a) 

.0680 

MOD -2 

(A- 7b) 


where 0 is the deflection angle of the particle. 

The intrinsic angular resolution derives from the approximately 
Gaussian distribution of the measured spark locations about the true 
trajectory position in each module. If we let represent the standard 
deviation of the deflection angle due to the intrinsic angular 
resolution, then the standard deviation of the measured angle is given 

by 

°0 =\Ap 2 + ’A 2 <A - 8) 

In Figure A-4 we show the angular distribution obtained for 
790 MeV positrons (nominal 0 = .0045 radians) at the Caltech Synchrotron 
using the detector configuration MOD-1. The smooth curve is a least- 
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squares fit of the data to a Gaussian distribution; the best-fit 
standard deviation is = .0023 radians, which implies = .0021 
radians from equations A- 7a and A-8 . However, the calibration runs 
were made with the beam aligned with the detector system, whereas 
during a flight particles have incident angles of as much as 30°. The 
distribution of the measured spark locations about the true trajectory 
position is expected to be broader at larger incident angles since the 
ion pairs are distributed over a large transverse distance. This argu- 
ment is substantiated by the fact that the average of the mean deviations 
of the measured spark; locations about the least-squares fit trajectories 
is approximately 60$ larger for flight data and ground-based muon runs 
than for the calibration runs. Therefore, we expect that a appropriate 

A 

for flight data might be slightly larger than that deduced from the 
calibrations. 


It is possible to determine o directly from the flight data. 

A 

More than 90$ of the particles which trigger the detector system during 
a flight are nuclei with energy greater than 400 MeV/nucleon, the 
threshold of the lucite Cerenkov counter. Of these particles approximately 
90$ are protons and 10$ He nuclei. For these particles the effect of 
scattering is small and, hence, their distribution reflects the intrinsic 
angular resolution described by 

2 

The expected rate of protons (p/m sec sr) in the deflection 


interval 0. 

i 


is 


N. 

l 



dR j(R) P(0,0’) 


(A-9) 
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where 

' 2 

j (R) = differential rigidity spectrum of protons (p/m sec sr MV) 

and P(0,0’) is the Gaussian probability distribution function defined 
by equation IV-1. In Figure A-5 we show the angular distribution of 
particles observed during the local nighttime interval of flight 71C2 
with deflection angles in the range -.008 to +.015 radians. The smooth 
solid curve is the predicted angular distribution of cosmic-ray protons 
calculated from equations A-9 and IV-1 using a proton spectrum appropriate 
for solar maximum and = .0025 radians. The dashed curve represents 
a similar calculation except that the proton spectrum appropriate for 
solar minimum is used. It is seen that the variation in the proton 
spectrum over the solar cycle shifts the peak of the distribution but 
does not significantly alter the width. Thus the standard deviation, 

° , which we derive in this manner does not strongly depend on the 

n 

assumed proton spectrum. Curves were calculated with the solar maximum 

proton spectrum for several different values of In each case the 

location of the peak of the calculated distribution was shifted to match 

the peak location of the observed data. For each distribution we then 

- 2 

calculated chi- squared, ^ , defined by 

X 2 = S [(y ± - n^A/n.] 2 (A-10) 

where 

til 

n. = observed number of particles in i channel. 

1 (Each channel is .001 radians wide.) 

y. = calculated number of protons in the i^ channel. 
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A plot of the resolution versus rigidity is shown in Figure A- 7. The 
filled circles represent resolution measurements made with positrons at 
the Caltech Synchrotron in configuration MOD-1 . 
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3 . Raw Flux Parameters 
a. Live Time (t ) 

The detector is insensitive for a fraction of the time 
interval over which the data are summed in the raw flux computation. 

For example, the phase-B one-minute rate counting period occurs every 
16 minutes. In addition, the total live time, t , during the 15 minute 

Li 

phase-A period is given by: 

t = (900 -nt) (1-nt) seconds (A-13) 

L e e cl q . 

where . 

n = number of events recorded 
e 

t g = time required to write a word (.35 sec) 

n •= total guard counter rate (cts/sec) 

3 

t = dead time following anti-coincidence (2 psec) 

3 

provided 


n t « 1 
a a 

which applied throughout all flights. The fractional dead time during 
a typical phase-A period at float altitude ranged from .14 in 1969 to 
.22 in 19 71. 

b. Spark Chamber Efficiency (D^) 

v 

Since every triple coincidence, T1AT2ALC, is a potentially 
valid event, the spark chamber detection efficiency is: 

D - f (A- 14) 


where 



n is the number of analyzable events 
N is the total number of triple coincidences 


In the case of MOD-2 data, the detection efficiency can be defined in 
terms of gas Cerenkov events, i.e. T1AT2ALCAGC, 


d g£ N 


n G? 


gE 


(A-15) 


where 


n v is the number of analyzable G? events 
GG 


N^< is the total number of GC-event coincidences 


Typically, D was .7 - .8 for a flight and D jf was approximately given by 

GG 

D ~ * D + -05 (A-16) 

GG 

This small difference is probably due to the different species of 
particles which make up the G& and NON-G? events. (Less than 10$ of 
the events are GC events; roughly 3/4 of these are electrons above 
15 MeV and the remaining 1/4 are nuclei above ~22 GeV/nucleon. The bulk 
of events triggering the detector are of the NON-GC type and consist 
mostly of protons above ~400 MeV.) Further analysis is in progress to 
determine the exact cause of the difference between D and D^* However, 
since the difference is small compared to the statistical accuracy of 
the data and since D can be obtained over short time intervals with much 
greater statistical accuracy than D v, we adopt the following method for 
estimating the spark chamber detection efficiency: we first determine 

D for the ascent and float intervals and then apply a correction factor 
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based on the difference of D and D v over the total float period (where 

uL 

better statistics prevail). 

V V 

c. Gas Cerenkov Efficiency Factor (C 

The gas Cerenkov counter was fabricated in January, 1970. 
Shortly afterward, it was calibrated at the Caltech Synchrotron as 
described in Section II. B. The efficiency was determined to be 
approximately 98$. However, oxidation of the mirror surfaces that 
reflect the light inside the counter can degrade the efficiency. 

V 

Although ground-based muon runs provide a check of the Cerenkov counter 
operation, changes in the efficiency of less than 10$ are masked by the 
statistical accuracy of the data. A comparison of ground-based muon 
runs in 1970 and 1971 show no significant differences; however, in order 

V 

to correct for possible smaller changes in the gas Cerenkov counter 
efficiency we use the flight data to directly calculate the efficiency 
factor. As an example, we show in Figure A-8 the hourly count rate of 
both gas-^erenkov and non-gas-Cerenkov events in the lowest three 
energy ranges for flight 71C2. We make the following interpretation of 
the particles making up the two classes of events 

V 

1) The GC events in the energy ranges considered consist 
of electrons only. 

V 

2) The NON-GC events consist of misfit nuclei, back- 
ground produced by nuclei, and also the electrons 

V 

which did not trigger the GC counter because its 
efficiency is less than 100$. 

V/ 

Under the above interpretation, a night-day difference in the NON-GC 
count rate is due to its electron component. All other N0N-G(5 events. 


t 
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which always have rigidities above cutoff, should retain a constant 
count rate from night to day. The electron counting rate increases by 
a factor of 3 or 4 from night to day due: to the large flux of re- 
entrant albedo electrons present at these low energies during the day- 
time interval. In the following we use the two-component model for 
the NON-g£ events and compare the night-day ratios of both classes of 
events. From these ratios we shall determine the fraction of electrons 
in the NON-GC data. Thus, both the count rate of g 6 electrons and 
N0N-G?5 electrons can be determined and, hence, the gas Cerenkov counter 
efficiency can be computed. 

We define the following symbols: 

e = nighttime GC rate (electrons) 
n 

e d = daytime GC rate (electrons) 

— V 

e^ = nighttime NON-GC rate of electrons 

— V 

e , - daytime NON-GC rate of electrons 
a 

p = time -independent rate of all other NON-GC events 


y “ = 


e d + p 


e + p 
n 


= ratio of day to night NON-GC rates 


a d 6 d 

x = — - = =- = ratio of day to night electron rates 
n n 

Both quantities x and y are determined from the data. It follows that 

SZ 

the ratio of the two components of the NON-GC class of events is given 
by: : 

-E- = 2L_-..y 

- y - 1 

n 


Since we measure the counting rate at night for NON-GC events, p + e^, 
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we derive the counting rate of the electron component of the NON-GC 


events at night : 

p + e 

_ n 

e„ = Z 
n £- + 1 

e 

n 


V V 

The nighttime GC rate, e^, is measured and thus the gas-Cerenkov 

v 

efficiency at all energies above GC threshold is: 


v 


C 


eff 


e 

n 

e + e 
n n 


Table A-l shows the result of this calculation for the flights of 1970 
and 1971. The error bars result from the statistical errors in the 
average day and night counting rates. The night-day transitions are 
most distinct in flights 70C2 and 71C2, and hence, we take as Cerenkov 

V 

efficiency factors, c e ff> the values .93 and .84 for 1970 and 1971, 
respectively . 
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4. Background Corrections ■ 

a. Upward-Moving Particles 

With the detector in the MOD-1 configuration, a fraction of 
the upward -moving particles (splash albedo and those due to y-ray 

V 

interactions in the lucite Cerenkov counter) contribute a small 
contamination to the low-energy data (due to the ~4$ backward detection 
efficiency in l£) . These corrections have been described in connection 
with the 1968 data (Rice, 1969, 1970). 

As a result of further calibrations at the Caltech Synchrotron 
we have made improved estimates of the y-ray contamination. The new 
corrections are given in Tables V-l and V-2 (Chapter V). The magnitudes 
of these corrections are not significantly different from those used by 
Rice (1970); however, the estimated errors have been reduced because 
of the more extensive machine calibrations. 

These calibrations also showed that the -'•4$ of backward- 
moving electrons which trigger the detector are rejected from analysis 
about twice as often as forward -moving electrons. Therefore, we have 
correspondingly adjusted the earlier splash-albedo corrections of Rice 
(1970). If we assume the typical detection efficiency of 0.75 for for- 
ward-moving particles (see Appendix A.3.b), it can be shown that the 
earlier corrections of Rice should be multiplied by the factor 2/3. The 
new splash-albedo corrections are listed in Tables V-la and b and V-2a 
and b . 

No corrections for upward -moving particles were necessary for 
the MOD-2 configuration (1970 and 1971) since the gas Cerenkov counter 
completely discriminates against these particles. 
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b. Atmospheric Muons and Pions 

The contamination of the MOD-1 data due to atmospheric muons 
and pions was shown to be negligible by Rice (1970). In the case of 
MOD-2 data, the muons and pions must have energy greater than about 
2.5 GeV and 3.4 GeV, respectively, to be above the effective threshold 

v 2 

of the gas Cerenkov counter. At 2.4 g/cm altitude pions of energy 

2 

3.4 Gev decay within about 0.06 g/cm of their point of production. 
Therefore, the flux of pions compared to that of muons may be ignored. 

From the pion production spectrum of Perola and Scarsi (1966) 

'.f 

and the formulas of Verma (1967), we have calculated the muon spectrum 
2 

at 2.4 g/cm (see Rice (1970) for details of the method). Above 2.5 GeV 
this spectrum is approximately: 

j (T) = .853 t" 2,94 u^/m 2 S ec sr GeV (A-17) 

where T is muon kinetic energy in GeV. Most of these high-energy 
particles have smaller bending angles than the .006 radians threshold 
value used in the data analysis. Folding the spectrum given by 
equation A-17 with the resolution function (equation IV-1) we find only 
2$ of these particles fall within the deflection interval 
.006 ^ |0| ^ .009. The percentage contribution to any of the other 

intervals is much smaller. We derive a count rate in the .006 - .009 
radians interval of 5.7 x 10 2 prVsec which is less than 0.2$ of the 
measured count rate of positrons or negatrons in this interval and there 
fore negligible. 

v 

c. Secondaries Produced in the Gas Cerenkov Counter 

2 

The 2 g/cm of material above the upper spark chamber is 
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potentially a source of contamination arising from interactions by 

cosmic-ray nuclei and yrays. We treat first the problem of nuclear 

interactions . '■ 

We consider two possible ways in which products of nuclear 
v 

interactions in the GC counter could simulate electron events 
(T1 AT2ALCAGC): 

1) Two or more particles (pions or protons) above 

V 

LC threshold could be produced with one traveling 

V 

through the detector system triggering LC, while 
another passes through one of the 1/2-inch 
quartz windows (which protect the GC phototubes) 

V 

triggering the GC counter. 

. V 

2) A high-energy particle above GC threshold could 
be produced in the material 1 above the flat 
mirrors of the GC counter and travel through 
the detector system triggering both Cerenkov 

V V 

counters (LC and GC) . 

The first possibility requires that one of the particles be emitted 
at a relatively large angle (order of 90°). Using information 
in the tables of Bertini (1967) on the angular and energy distribution 
of secondary protons and pions from interaction p+O^ (fluorine and 
sulfur tables were not available), we estimate that the upper limit on 
the rate of such events is ~5xl0 ^ particles/sec in any one energy range. 
This rate is less than 1$ of the measured count rate and therefore neg- 
ligible. In the second case we are only concerned with particles above 

v 

the effective gas Cerenkov threshold (Section II. B). We illustrate 



138 

the magnitude of the correction by considering the production of pions . 
The mean lifetime of 3.4 GeV charged pions is 6.1 x 10 ^ seconds. Thus 
a pion of this energy will travel roughly 180 meters before decaying. 
Since the entire spectrometer is just over a meter long, we ignore 

V' 

pion decay in the calculation. In order to trigger the gas Cerenkov 

2 

counter, the particle must be produced in the 1 g/cm of material 

above the mirrors. Considering the rather complex geometry of the 

detector, we shall calculate an upper limit to the contamination. We 

v 2 

first replace the GC counter by a 1 g/cm slab of air. Then we use 

the pion production spectrum of Perola and Scarsi (1966) to calculate 

the flux of pions that emerge from the bottom of such a layer which is 

exposed to the cosmic-ray nuclei flux. Above 3.4 GeV, the 

2 

differential flux of pions at 1 g/cm can be represented by the 
power law: 

j (T) = 16. 5T 2*35 sr Ge v (A-18) 

When we fold this spectrum with the resolution of the detector we derive 
a count rate of 1.85 x 10 ^ jt^/sec in the highest energy interval. 
Considering both charge signs, this represents only 2^ of the measured 
count rate for positrons in this energy interval in 1971 and less than 
.5$ of the negatrons . In addition, the high-energy nuclear inter- 
actions which produce the pions have a high multiplicity and even the 
very small contamination derived above is drastically reduced because 
of the anti-coincidence counters and the multiple-spark-indicator 
feature. Therefore we consider the contamination from nuclear inter- 
actions negligible. 

Cosmic -ray nuclei also produce knock-on electrons as they pass 
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through the gas Cerenkov counter. When expressed in terms of g/cm the 
probability functions, describing knock-on production are proportional 
to Z/A where Z and A are the charge and mass numbers of the material 
traversed. Thus the production spectrum' of knock-ons will be essentially 
material-independent. We use the production spectrum derived by 

Beuermann (1971) for air appropriate to the cosmic-ray nuclei flux 

j 

level of 1968. Since knock-on electrons of 15 MeV are produced by 
protons with energy greater than ~3 GeV, where solar modulation effects 

: I 

are not large, the use of the 1968 proton spectrum for the period 1968- 
1971 introduces negligible error. Beuermann' s production rate for the 
interval 10 ^ T ^ 100 MeV is given approximately by: 

I 

-2 7- 

Q(T) = .035 T * e /g sec sr MeV (A-19) 

2 

This rate corresponds to a flux of ~.10 e /m sec sr MeV in the lowest 

energy interval (14.3 - 29 MeV). This flux represents some 25<£ of the 

2 

measured flux in this energy interval at 2.4 g/cm in 1971. However, 

using the formulas of Rossi (1952), we calculate that a 15 MeV knock-on 

electron emerges at no more than 3 1/2° from the forward direction. 

Hence, we would expect a very large fraction of these events (proton + 

knock-on electron) to' be multiple-particle events which are rejected 

2 

from analysis. In addition, the 1 g/cm of material is an upper limit 

since particles produced near the mirrors will not have sufficient 

v 

pathlength to trigger the gas Cerenkov counter. Taking these effects 
into account, we consider the contamination due to knock-on electrons to 
be negligible. 

A possible source of contamination could arise in the high- 
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v 

energy intervals if a proton above lucite Cerenkov counter threshold 

V 

produces a T1AT2ALC coincidence and also produces a knock-on electron, 

. V 

which triggers GC but which fails to be registered in the spark chambers. 
This could happen, for example, if the electron experienced a large 
single scattering or if the multiple-spark-detection efficiency were 
less than 100$. We estimate from the knock-on production spectrum that 

V 

less than 0.1$ of the nuclei above the lucite Cerenkov counter threshold 
will produce a knock-on electron of sufficient energy to trigger the gas 

V 

Cerenkov counter. However, approximately 4$ of the cosmic-ray protons 

V V 

above the LC threshold are also above the GC threshold. Thus the above 

contamination is only ~l/40 of that due directly to the high-energy 

nuclei (see Appendix A.4.d) and therefore negligible. 

Electrons can also be produced in the gas Cerenkov counter 

by Compton scattering and pair-production from y-rays . In order to 

calculate the fluxes of positrons and negatrons resulting from these 

effects we assume the following: 

2 

1) The 1 g/cm of material above the mirrors has an average 
charge number, Z, of 10 and an average mass number. A, of 
20 . 

2) The probability functions of Rossi (1952) are used. In 

the case of pair production we use the complete-screening 

approximation which gives a larger flux of electrons for 

our conditions than does the no-screening approximation. 

2 

3) We use the atmospheric y-ray spectrum at 2.4 g/cm 
residual atmosphere from the calculations of Beuermann 

! (1971) with the electron and cosmic-ray nuclei spectra 
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appropriate for 1971. 

In Figures A-9 and A-10 we show the assumed Y-ray spectrum and 

the resultant positron and negatron spectra, respectively. For 

negatrons the sum of the Compton scattering and pair-production processes 

is shown. In the 14.3 - 29 MeV interval we derive positron and 

2 2 

negatron fluxes of .010 p/m sec sr MeV and .012 p/m sec sr MeV, 
respectively. These values represent only ~4<^ of the measured fluxes 
at float altitude in 1971 and therefore represent a negligible contribu- 
tion. 

d. High-Energy Cosmic-Ray Nuclei 

Cosmic-ray nuclei above approximately 22 GeV/nucleon also 

V 

trigger the gas Cerenkov counter. Most of these particles are confined 
to bending angles smaller than those used in the calculation of electron 
fluxes (i.e. <; .006 radians). Because of the resolution of the detector, 
however, a small fraction of these particles are observed with larger 
bending angles. As ah example, we show in Figure A-ll the angular 

V 

distribution of GC events for the total float period for 1971. The 
dotted curve represents the calculated distribution of the sum of 
primary and secondary electrons. The remaining events are high-energy 
nuclei. The smooth curve is a Gaussian distribution with o -= .0025 
radians and this curve was used to calculate the contamination in the 
highest energy interval. The results for 1970 and 1971 are given in 
Tables V-3b, V-4b, V-5b, and V-6b. The proton contamination in the 
98 3 - 1475 MeV interval is ~25=£ for positrons and ~6$ for negatrons. 

The contribution is negligible at lower energies. The error in the 
values is estimated at ~50^ on the basis of uncertainties in the 
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instrument resolution and spark chamber alignment. 

v 

e. Accidental Gas Cerenkov Coincidences 

About 90$ of the particles which trigger the detector are 
cosmic-ray protons above ~1000 MV, the threshold of the lucite 

V V 

Cerenkov counter. Accidental GC coincidences tag some of these 

particles as electrons. Since 1000 MV corresponds to a nominal 

9 ~ .009 radians, these events primarily contaminate the high-energy 

positron data. We have used the measured rate of cosmic-ray nuclei, 
v 

the accidental GC -coincidence rate (Section II. b), and the resolution 

of the detector to calculate the expected rate of these events in the 

highest energy intervals. The results are given in Tables V-3b and 

V-5b. The maximum contribution to the data is 16$ in the 

983 - 1475 MeV positron interval in 1971. The 50^ estimated error is 

based on uncertainties in detector resolution and spark-chamber 

alignment. 

f . Spark Chamber Alignment 

Because the fiducial wires cannot be precisely lined up for 
all 8 planes it is possible to have a built-in offset in the deflection 
angle. An initial alignment is made by using the deflection-angle 
distribution of the cosmic-ray nuclei above lucite Cerenkov counter 
threshold. This procedure is similar to the one described in 
Appendix A. 2 in determining the instrument resolution. However, in 
the resolution calculation we were interested in determining the width 
of the deflection-angle distribution of the nuclei. In the alignment 
procedure we are interested in comparing the locations of the peaks 
of the measured and calculated deflection-angle distributions. These 
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distributions are shown in Figure A-5 for flight 71C2. Alignment 
factors have already been introduced in the data in computing the 
observed histogram. The two curves were calculated using solar 
minimum and solar maximum proton spectra, respectively. Since the 
peak in these curves differ by only .0015 radians over a solar cycle, 
we feel that our deflection zero is accurate to approximately .0005 
radians . 

In the case of MOD-2 high-energy data the alignment was 
refined by using the G? proton distribution (see Figure A-ll) . The 
average bending angle of these particles (T > 22 GeV/nucleon) is 
calculated to be ~0. 00018 radians. It was necessary to adjust the 
data by about .0005 radians in 1970 and 1971. 


In Table A- 2 we summarize the corrections to the data 

discussed in this section (A. 4). Only the last two entries (high- 
s' 

energy protons and accidental GC coincidences) are considered non- 
negligible. These two corrections and those due to atmospheric 
secondaries are included in the data tables of Chapter V. 
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TABLE A- 2 


Summary of Corrections Discussed in Section A. 4 


* 


Probable 


Description of Correction 
(MOD -2) 

Section 

Discussed 

Energy Intervals 
Affected 

Contribution 

(tf) 

Upward -moving particles 

A. 4. a 

low energies 

0 

Atmospheric muons and 
pions 

A .4 .b 

983-1475 MeV 

0.2 

Products of nuclear 
interactions in G? 

A.4.c 

medium and 
high energies 

<2 

Knock-ons in G? 

A.4.c 

low-energy e 

~0 

V 

Knock-ons in GC scattered 
out of acceptance cone 

A.4.c 

high energies 

0 

Compton- scattered electrons 
and pair-produced electrons 
in G? 

A.4.c 

14.3-29 MeV 

<4 

Cosmic-ray protons above 
G? threshold 

A.4.d 

98 3-1475 MeV 

25 e + 
6 e" 

V 

Accidental GC coincidences 

A .4 . e 

98 3-1475 MeV 

16 (e + -1971) 


Unless otherwise noted a correction applies to both charge signs. The 
correction for spark chamber alignment was made before computing the 
raw fluxes (Chapter V) and hence no correction for this effect is listed. 
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5 . Correction of Fluxes to The Top of The Atmosphere 

In the final stage of data analysis we correct the locally 
observed fluxes to their values at the top of the atmosphere. This 
correction is complicated by the energy loss (bremmstrahlung and 
ionization) experienced by the particles in the material above the 
spectrometer and by the decreasing resolution of the detector at high 
energies (or, equivalently, small bending angles) . In computing the 
"raw" fluxes at float altitude (data tables of Chapter V), we have not 
taken into account the resolution, i.e. we have assumed that the 
resolution is perfect and have calculated the flux by dividing the 
measured rate in a given energy interval by the width of that interval 
(and by the geometrical factor) . However, in the highest two energy 

°Q 

intervals of the MOD-2 observations the ratio —x~ is relatively large 

A °i 

(o is' the standard deviation of the deflection-angle distribution 

C7 

til 

given by equation A-llb and is the width of the i deflection-angle 

bin given in Table IV-1), and there is considerable probability, 
particularly above 1 GeV, that a particle is assigned to the wrong 
energy interval. Thus, at high energies, dividing by the numerical 
width of the energy interval is not necessarily a correct way of 
relating the observed rate to the true particle flux. We describe a 
procedure which accounts for the effect of resolution as well as 
energy loss in correcting the measured fluxes to the top of the 
atmosphere. The method is based on a similar calculation by 
Fanselow (1968). 

We begin with a trial primary electron spectrum incident at 

■ i '■ ■ 

the top of the atmosphere. (Such a spectrum is also used in the 
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secondary subtraction procedure described in Chapter IV.) We then 
write in matrix form the equation relating the expected rates at float 
depth to the assumed incident spectrum. We invert the resulting 
matrix and thus derive the spectrum at the top of the atmosphere from 
the observed rates at float altitude. Our procedure also allows us 
to calculate the uncertainty in the fluxes at the top of the 
atmosphere. Because of the effect of resolution, these uncertainties 
in the data of the two highest energy intervals are larger than those 
of the raw fluxes at float altitude. 

The assumed incident primary spectrum j^CT) is modified by 

energy-loss effects as the particles pass through the atmosphere and 

2 v • 

the ~2 g/cm of material between the top of the gas Cerenkov counter 

and the upper spark chamber. We denote the residual primary spectrum 

at the top of the upper spark chamber by j (T) . We calculate j (T) 

from j.j.(T) by folding in the bremmstrahlung energy-loss probability 

distribution and by including the average ionization energy loss. 

We first calculate the effect of bremmstrahlung energy loss. 

The probability P(T,T')dT' that a particle with kinetic energy T will 

have energy between T' and T' + dT' after passing through x radiation 

lengths of material is given by Rossi (1952): 

Y-l 

P(T,T') dT' = dT' ^ (A-20) 

where y - x/^n 2 and r(y) is the gamma function (Abramowitz and Stegun, 

1964). The spectrum after correcting for bremmstrahlung loss is 

00 

j(T') = J j ] .(T) P(T, T' ) 

T' 


dT 


(A-21) 
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We then obtain jg(T), the residual spectrum at the top of the upper 
spark chamber, from j(T') by correcting for ionization loss using the 
tables of Berger and Seltzer (1964) . 

The expected rates at float altitude are obtained through the 
response of the magnet spectrometer to j (T) . We have 

b 


where 


M. 

1 


f R^T) G(T) j g (T) dT 
0 



R i <T) G(T) jg (T) dT 


(A-22) 


PL = expected rate in the i th energy interval (p/sec) 

R_ (T)dT = probability that particle with kinetic energy T 

til 

will be observed in the i energy interval 
G(T) = geometrical factor as a function of T 


S S til 

Tj, T^ + ^ = lower- and upper-limit energies defining the j 

energy interval . 

I = number of energy intervals. (Electrons with energy 

greater than the largest measured interval are placed in 

a bin from 1.475 - 40 GeV.) 

S S 

The superscript S on T and denotes the fact that the quantities 

are measured at the upper spark chamber. 

The response function IL (T) in equation A-22 is related to 
the resolution function described earlier in Appendix A. 2. P(0,0')d0', 

the probability that a particle with rigidity corresponding to a bending 
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angle 9 will be observed with bending angle between 6 ' and 9' + d0', 

. ■ . 4 - • 

is given by 


P(0,0')d0' 


d<9' -(9-9') 

exp 

0 


(A- 23) 


where 9 and o are given by equations II-2 and A-ll, respectively. 
Then the response function is 


R.(T) 



P(9,9 ' )d0 ' 


(A- 24) 


x 

(Note 9 is a function of T by equation II-2) 

In Figure A-12 we show plots of versus T for each energy interval 
for the MOD-2 configuration. A curve corresponding to perfect 
resolution would be a rectangular box with amplitude 1 and with 
vertical sides at the two energies corresponding to the energy-interval 
end points, which are indicated by dotted vertical lines in the figure. 
For low energies the resolution is good while for the highest energy 
interval used in the analysis the relatively large ratio of 
causes Ry(T) to significantly overlap the adjacent energy intervals. 

The of equation A-22 represent the expected rates at float 
altitude for the assumed incident primary spectrum, taking into account 
the energy loss of the electrons and the resolution of the detector. 

It proves useful to write an equation for the corresponding rates at 
the top of the atmosphere . N^, assuming that the detector has perfect 
resolution and that no energy losses occur in the G(? counter: 


1 


T 1 

i+1 


N. 

l 


G(T) j J .(T) dT 


(A- 25) 
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The superscript I denotes the fact that the quantities are measured at 
the top of the atmosphere. It is usual in the literature to shift the 
energy intervals by the ionization energy loss which occurs for each 
particle. The breiranstrahlung loss distribution is such that most 

I 

particles lose very little energy and hence an average shift for this 

effect is not applied. Thus, the energy at the top of the atmosphere, 

I S 

T^, is related to the energy at the upper spark chamber, T. , by 

T* = T^ + (A-26) 
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If H. . is an element of the inverse of the F matrix, the rates at 
the top of the atmosphere are given by: 


N. = 
1 


i 

E 

j=l 


H. . M. 
ij J 


Furthermore: 




i 

E 

j=l 




(A- 30) 


(A-31) 


where cr is the standard deviation for the observations and cr i s 
M. N 

J i 

the resulting standard deviation for the N. . 

l ■ 

Equations A- 30 and A-31 represent the desired results. 

Equation A-30 gives the corrected rates (N^) at the top of the atmosphere 
in terms of the rates at the detector at float altitude (M. ) . If we 

l 

replace the M. in equation A-30 by the actual measured rates, . , then 

we derive a set of from which a better approximation for the trial 

input spectrum can be made. From this new trial spectrum the F„ and 

H. . are recalculated and a new set of N. derived. The process can be 
ij i 

repeated until the N^, converge. (Note that since the rates, N^, 
correspond to a detector with perfect resolution, it is appropriate 
to divide by the width of the energy interval in converting the rates 
to differential fluxes.) However, in some cases the measured rates are 
only upper limits and in these cases the matrix- inversion procedure 
for determining the corrected rates cannot be applied. Therefore, we 
modify the procedure for determining the corrected fluxes in the 
following way. For a given trial input spectrum we determine the ratios 
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N. 


til 

— , where N. and M. are the calculated rates in the i energy interval 

M. ’ ' x r 

1 

at the top of the atmosphere and at float altitude, respectively. If 

fljiyi.. are the measured rates at float altitude, then we determine the^H . , 

the estimated corrected rates at the top of the atmosphere, by 
N i 

. = r— . These^fl . are used to define a new trial input spectrum 
* x M . i x 


x 

and the procedure is repeated until the^YV are consistent with the 
input spectrum. 

We do not modify the procedure (equation A-31) for calculating 
the expected uncertainties in the fluxes at the top of the atmosphere. 
This procedure takes full advantage of the matrix- inversion technique 
and yields larger errors in the corrected high-energy data than those 
listed for the raw fluxes (data tables of Chapter V), which are based 
on statistical errors only. In deriving ! the errors we need a knowledge 
of the error in the data for each of the energy intervals, including the 
1.475 - 40 GeV interval. The number of electrons in this interval is 


not measured directly since a significant fraction of the gas Cerenkov 
events with |0| ^ .006 radians are high-energy protons (j>22 GeV). It 
is possible to roughly deduce the number from the observed deflection - 
angle distribution of these events and from a knowledge of the resolu- 
tion of the detector (see Figure A-ll). We estimate the error in the 
data from this interval to be ^25$. In making the corrections we have 
assumed a 25$ error in this interval (>1.475 GeV) which increases the 
relative error in the highest measured energy interval (0.983-1.475 GeV) 
by ~25$ over the result obtained if zero error is assumed for the 
integral flux above 1.475 GeV. The error in the data at lower energies 
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is not affected by the estimate. 

We found that the above matrix- inversion procedure was not 
necessary for the MOD-1 observations, which extend to only 200 MeV, i.e. 
the matrix was essentially diagonal because of the excellent resolution 
at low energies. The fluxes corrected to the top of the atmosphere 
for both MOD-1 (1968 , 1969) and MOD-2 (19 70, 1971) observations are 
given in the tables and figures of Chapter V. 
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APPENDIX B 

Interstellar Electron Spectrum from Non-Thermal- 
Radio-Background Data 


The synchrotron emission (ergs/sec -cm .sr-Hz) by electrons 
spiraling in a uniform magnetic field is given by Ginzburg and 
Syrovatskii (1964) as: 


e(v) 4it 2 


*/3 e\ f W 2 - 1/2 

' [1+ (1-T1 ) (-“) ] 

me J 2 me 

me 


-*7 f K 5/3 <”> / N 

v c Jv/v„ ' ! 


(W)dW 


(B-l) 


where v is the frequency of emission 
e is the charge of the electron 

B is the mean value of the magnetic field perpendicular to the 

x 

line of sight 

m is the mass of the electron 

c is the speed of light 

W is the total energy of the electron 

' I 

N(W) is the number density of relativistic electrons 
K 5 / 3 (T|) is a modified Bessel function of the second kind 

■\IT- 


T| 


ne 


is the refractive index of an ionized gas with 


rtm v 


n electrons per cm and 
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W 2 

V c - v c [1 + (1-H ) 2> ] 

me 


-3/2 


(B-2) 


154 


where v £ is the critical frequency given by 


3eB 


= L 

v c 4jtmc 2 } 
me 


(B-3) 


The terms involving l-^ in equations B-2 and B-3, which result from the 
ambient electron density in the interstellar medium, cause a suppression 
of the emission at low frequencies (Razin, 1960; Lerche, 1971). 

For the case of a vacuum (Tj = 1) equation B-l has been solved 
exactly for an electron spectrum which is a power law in energy over a 
sufficiently large range (Ginzburg and Syrovatskii, 1964). The result 
is 

111 i±2 

3 3 5 ) 2 C B 2 v 2 (B-4) 

2im c -L 


(v) = 


e kv 


£_ 

l+Y 


r<^)r(^> 


12 


4itmc‘ 


where C and y are parameters defining the electron spectrum: 

N(W) = C W _Y 

and r(x) is the gamma function (Abramowitz and Stegun, 1964). 

To obtain the total intensity of radiation over a given line- 
of-sight distance we must include the free-free absorption by the 
medium. The absorption coefficient for the radio frequencies of interest 
to us is given by: 

? n 2 r 3/2 

k (v) = 10 -377—2 c 1 7 - 7 + £ n (B- 5 ) 

T v 

(Ginzburg, 1964) where T is the electron temperature. The total 
intensity is then given by solving the differential equation: 


77 = e (v) - k( v )I 


(B-6) 
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where ds is an element of length along the line of sight . The solution 
of equation B -6 depends on the structure of the interstellar medium. 

Several recent reviews have dealt with the physical properties 
of the interstellar medium (Field, 1970, 1971; Heiles, 1971; Dalgarno 
and McCray, 1972; Wentzel, 1972). The observations indicate a medium 
composed of dense, cold clouds with diameters of a few parsecs separated 
by a hot, rarefied intercloud region with a scale of hundreds of parsecs. 
This structure has also been predicted by theoretical studies (Field 
et al., 1969; Hjellming et al., 1969; Shu et al., 1972). We shall take 
as a model of the interstellar medium a uniform distribution of cold 
clouds with diameter and separation distance (This model with 

£ = 1 pc and j 0 . = 1 kpc is identical to that used by Goldstein, Ramaty 

and Fisk (1970)). The subscript convention is "c" for "cold" clouds and 
"i" for "intercloud" . The first cold cloud is assumed to lie at a 
distance In Figure B-l we show a schematic diagram of the assumed 

galactic structure. Our position in the galaxy is labeled s. We assume 
there are m clouds (and thus m intercloud separations). 

In the hot intercloud region we have both emission and 
absorption. If there were only one such intercloud region, the 
solution to equation B -6 would be 


- T. 



where the optical depth t. is defined by 



0 


(B-7) 


(The optical depth of a cold cloud is similarly defined by changing the 
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subscripts to c.) 

Since, typically, X c « we may assume the emission in cold 

clouds is negligible. Solving equation B-6 for e = 0 implies that the 

radiation penetrating through a cold cloud is partially absorbed with 

" T c 

an absorption factor e . Thus the solution of equation B-6 for the 
assumed galactic structure (see Figure B-l) can be written as the 
following series of equations 


GALACTIC POSITION 
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2 
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2m 


I 2m-1 ° 


^Zm = ° 


where we have assumed that there are an integral number of intercloud 

separations with — = m = integer where L is the total emission distance. 


By substitution we obtain: 
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The factor r^-Q-e ) is common to all terms. We thus have 
k. 
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r s - *> t 1+e 
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The term in braces is just 
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Hence, the solution to equation B-6 for our model is: 


t. -“(t.+t ) 

( 1 -e *) ri-e ] 
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(B-9) 


At high frequencies k becomes small (equation B-5) and thus j is small 
(equation B-7). Equation B-9 becomes 




(B-10) 


Using equation B-4 for e ( v ) we obtain 

1±1 

2 -CL 

Kv) a C L B v (high frequencies) 

X 


V— 1 

where we define OL = -‘y - as the power-law index of the radio emission 
spectrum. 


At low frequencies, t is large and equation B-9 becomes 



and using equations B-4 and B-5 we have 
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i±2 

I(v) OC C B 2 v 2 ~\ 3/2 v{ 2 
4 


I 2 (low frequencies) (B-ll) 


T. 

[17- 7+^n— ] 


v 


Thus at low frequencies interstellar absorption changes the spectral 

2-a 

shape to roughly v • (Note that the logarithmic term varies quite 
slowly with frequency.) We note that the intensity at low frequencies 
is independent of the total line-of -sight distance, L. 

At intermediate frequencies and for the general case including 
the Razin effect and an arbitrary electron spectrum, we must solve 
equations B-l, B-5, and B-9 numerically. (The term in braces^ in 
equation B-l is available in tabular form (Ginzburg and Syrovatskii, 
1964)). If one knew the value of all the parameters involved in the 
equations - principally B , L, 1 ^, T^, T c , n^, and n £ - one could 

vary the energy dependence and magnitude of the electron spectrum until 
the computed radiation intensity matched the observations (Figure VI -8). 
However, there is considerable uncertainty in some of the parameters - 
particularly the intercloud ambient electron temperature T^ - and this 
uncertainty should be reflected in an uncertainty in the interstellar 
cosmic-ray electron spectrum. To illustrate the approximate range of 
interstellar electron spectra possible we: 

1) choose a nominal set of galactic parameters and calculate the 
galactic electron spectrum necessary to account for the radio 
emission in the galactic anticenter direction, 

2) vary each parameter through its range and for each variation 
calculate, using the nominal electron spectrum, the resulting 
radio -emission spectrum, and 
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3) compare these spectra with the observed radio-background emission. 
Since the emitted power is directly proportional to the number of 
electrons, we can derive a multiplicative factor versus frequency that 
can be applied to the electron spectrum once the conversion factor 
between electron energy and emission frequency is established. Using 
these multiplicative factors we can estimate the electron spectrum that 
is required to produce a radio-emission spectrum in agreement with the 
observations for a given set of parameters. 

We now discuss the possible range of the parameters involved 
in the calculation. 

1) B - Theoretical arguments concerning the dynamics of the galaxy 

x 

place the average magnetic field between 3 and 5 p-gauss. (Parker, 
1969b). Recent studies of dispersion and rotation measures 
observed for 18 pulsars indicate an average interstellar field of 
~3 ,5 p-gauss with an estimated error of .5p-gauss. (Manchester, 
1972). The value o ( f 3.5 p-gauss corresponds to the largest fields 
found and it is thought that the lines of sight to the two pulsars 
involved lie along the direction of the magnetic field. Moffet 
(1971) from a survey of the polarization properties of pulsars also 
derives magnetic fields of a few microgauss. In this study we use 
a nominal value of 5 p-gauss and illustrate the variation in the 
calculated radio spectra for the range 3-5 p-gauss. 

2) L - The emission length for a uniform galaxy would just be the 

distance to the edge of the galaxy. The size and structure of the 

I 

galaxy and our location in it have been deduced from optical and 
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radio astronomy observations (Blaauw and Schmidt, 1965; Allen, 
1963). The solar system is approximately 8-10 kpc from the 
galactic center and the diameter of the galaxy is ~25-30 kpc. The 
distance to the edge of the galaxy in the anticenter direction is 
estimated to be ~4-5 kpc. However, the galaxy is quite inhomogenous , 
being composed of spiral arms. Thus the use of an average emission 
length of 4 or 5 kpc might be misleading. We have chosen to use 
4 kpc as our nominal value but include the range L = 2-6 kpc in our 
calculations . 1 

3) and ^ - The average fraction of the line of sight intercepted 
by cold clouds depends on the viewing direction. A model of the 
interstellar medium which uses the data (dispersion measures and 
21 -cm absorption) from two pulsars, NP 0532 and CP 0328, which lie 
roughly in the anticenter direction, yield a range of (the 

"filling factor") of 0.006-0.081 (Hjellming et al . , 1969). 

Other models yield similar results - for example. Field et al . 

(1969) obtain a range 0.02-0.04 and Dalgarno and McCray (1972) 
in their recent review paper use values of 0.007-0.041 for the 
filling factor. We shall take the larger range 0.006-0.08 in the 
calculations with a nominal value of 0.02. The recent observations 
of 21-cm hydrogen absorption by Greisen (1973) indicate that clouds 
having dimensions of i ~ 1 pc are common. In our model we shall 
consider variations in of from 1-10 pc. The possible values of 
^ we shall use are 13-1670 pc corresponding to the assumed filling 
factor range of 0. 006-0 . t 08 . 
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4) T. and T c - The observations of 21-cm absorption yield an estimate 

* 

of the spin temperature of neutral hydrogen. The high resolution 

studies of Hughes et al . (1971) and Radhakrishnan et al . (1972) 

indicate that the spin temperatures in clouds range from ~15° - 250°K. 

This range is somewhat larger than those considered in the 

theoretical models . We shall take T = 250°K as the nominal value 

c 

and consider the effect of varying T c through this total range. 

There is little accurate observational evidence for the spin 
temperature of the intercloud medium but the lower limit is roughly 
500°K (from observations of emission features in which no detectable 
absorption is observed - Hughes et al., 1971; Radhakrishnan et al., 
1972). However, the spin temperature may be much lower than the 
gas kinetic temperature in the low density intercloud region 
(Dalgarno and McCtay, 1972). Upper limits to the kinetic 
temperature may be estimated from velocity dispersion measurements . 
Heiles (1967) has found an emission feature with a velocity 
dispersion corresponding to T^ < 4000°K, although Field (1971) finds 
none with upper limits less than 8000°K. Mention should be made of 
the theoretical models of the two-component system which generally 

t ! 

* 

fix the temperature with a rather high value. For example, Habing and 
Goldsmith (1971) use T^5000-8000°K. and Shu et al. (1972) obtain 
T. ~ 7500°K. The theoretical models, however, are not, as yet, 

* 

The spin temperature relates the number of hydrogen atoms in each of 
the two hyperfine ground state levels. In dense clouds, where 
collisional excitation is important, the spin and kinetic temperatures 
are expected to be equal. (Field, 1958) 
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experimentally well verified (Greisen, 1973). We shall take the 
range T\ ~ 500-10000°K and use the upper limit for our nominal 
value . 

5) n^ and n^. - The frequency dispersion of pulse arrival time from 

pulsars provides a measure of the density of free electrons along 

the line of sight. For example, for the Crab Nebula pulsar 

NP 0532 J n e dj£ = 57 cm 8 pc (listed by Maran, 1969). The distance 

to the pulsar is 2020 pc (Trimble, 1968) and hence we derive 
~ -3 

n g ~ 0.03 cm . Most theoretical models predict little contrast 

in the cloud and intercloud electron 1 densities. Field et al . (1969) 

use n^ = 0.02 in their model and Hjellming et al . (1969) use 

n^ ~ 0^03 and n c x: 0.04-0.05. Bridle and Venugopal (1969) find 

that n^ = 0.02 is consistent with a variety of data. Dalgamo and 

McCray (1972) intheir review use values of n^ ~ 0.03-0.05 and 

n ** 0.06 (for cosmic ray heating). We shall use the nominal 
c 

values n^ = 0.03 with a range 0.02-0.05 and n £ = 0.02 with a range 
0.02-0.06. 

In summary, we present in Table B-l the nominal values of the 
parameters together with their possible range considered here. 

Using the nominal values of the parameters we obtain an inter- 
stellar electron spectrum which yields from equation B-l, B-5, and B-9 

a radio spectrum in agreement with the observations. This spectrum can 

2 

be represented in p/ (m sec sr MeV) by 

C 1.34 x 10 4 W‘ 1,8 W < 2000 MeV 

2.75 x 10 6 W~ 2 ' 5 W > 2000 MeV 


j(W) = 


(B-12) 
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TABLE B-l 

Galactic Parameters Used in the Analysis 


Parameter (Units) 

B (p-gauss) 

1 

L (kpc) 

p ± (pc) 

i c (PC) 

T. (°K) 

T (°K) 
c 


Nominal 

5 

4 

50 

1 

10000 

250 


Range 

3-5 

2-6 

13-1670* 

1-10 

500-10000 

15-250 


n. (cm ~ > ) -03 .02-. 05 

n (cm ^) *02 .02-. 06 


*In our model the value of t. is determined by the assumed range of the 

l 

filling factor (ff = j .006 - .08. 
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where W is the total energy of the electron. 

Using this electron spectrum we first calculate numerically a 

* | 

rough correspondence between electron energy and emission frequency. 

An approximate analytical expression for the frequency of maximum 
intensity from an electron of energy W is: 

v (MHz) ~ 0.29 v = 4.6B W 2 (B-13) 

m c l 

where B is in p-gauss and W is in GeV (Ginzburg and Syrovatskii, 1964) . 
However, this correspondence does not take into account the spectrum 
of electrons, and since there is contribution to a given frequency from 
a considerable range of electron energies we might expect some 
deviation from equation B-13 for a steep interstellar spectrum. In 
Figure B-2 we show a calculation of the relative contribution to the 
synchrotron emission at v = 10 MHz from electrons with energies between 
100 MeV and 10 GeV using the nominal interstellar electron spectrum given 
by equation B-12. The peak contribution comes from W ~ 330 MeV. The 
arrow at 660 MeV corresponds to the prediction of equation B-13, which 
is based on a flat electron spectrum. The difference by a factor of 2 
in the estimates roughly applies throughout the range of observations, 
0.4-600 MHz. We note that 0.4 MHz, which is the lower limit of the 
radio observations, corresponds to ~60-80 MeV electrons. Thus below 
this energy the radio data provide little information on the inter- 
stellar electron spectrum. 

In Figure B-3 we show a plot (solid line) of the peak energy 
versus frequency which we obtained from plots similar to Figure B-2 
covering the entire range of radio data for the nominal 5 p-gauss field . 
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The dashed line refers to a calculation using B =3 g-gauss. We shall 

x 

use Figure B-3 as a guide in estimating the electron energies at which 
changes in the electron intensity are necessary in order to produce cal- 
culations of the synchrotron spectrum in agreement with the observations. 

As mentioned earlier , the change of a single parameter in 
equations B-l, B-5, and B-9 will produce* a change in the computed radio 
spectrum. If at each frequency we compute the ratio of the observed 
radio spectrum to the computed spectrum, we obtain a set of multi- 
plicative factors F(W) which can be applied, using Figure B-3, to the 
nominal interstellar spectrum to derive the adjusted galactic electron 
spectrum necessary to produce the observed synchrotron emission. We 
show in Figure B-4 the result of such a calculation. The nominal 
electron spectrum of equation B-12 has been used; each line corresponds 
to a change of the labeled parameter. 

The effects of the interstellar medium at the low frequencies 
(low energies) are clearly shown. The range of parameters considered 
implies that below ~300 MeV (corresponding to ~10 MHz) the interstellar 
electron spectrum becomes increasingly uncertain due to uncertainties in 
the properties of the medium. 

At high energies (> 300 MeV) only B and L cause uncertainties 

X 

in the interstellar electron spectrum. From equation B-10 we have: 

I( v ) a C B 1,4 L 
. x 

where we have set the interstellar electron spectral index y = 1.8 
(W < 2000 MeV) and C determines the magnitude of the electron spectrum. 


Thus 
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c a. 


1.4 


We can estimate the error in C from the range of values of B and L. 

x 

If we assume that these parameters are Gaussian distributed about the 
mid-points of their ranges with the end-points of the range taken as 

the l-o limits, we derive 

2 2 2 
a 


C _ .. , .2 ^B_ , 

C B L 


0.375 


or 


” = 0.61 (B-14) 

Below ~300 MeV the functional dependence of the calculated 
radio spectrum on the many interstellar medium parameters cannot be 
easily determined. At these energies it is not obvious how to perform 
an appropriate statistical analysis, and, instead, we characterize the 
variation by an envelope which encloses the maximum variation for a 
single parameter. The total envelope of variation considered, which 
at energies above ~300 MeV is determined from equation B-14, is shown 
as the dashed curve in Figure B-4. Our estimate of a reasonable range 
of interstellar electron spectra is obtained by multiplying the 
nominal electron spectrum (equation B-12) by the factors corresponding 
to this envelope. These high and low spectra are shown in 



Figure Vl-9 . The explicit energy dependence, represented by power- law 
segments, is given in Table B-2 for all three spectra. For certain 
choices of galactic parameters each of the electron spectra can produce 
a radio spectrum in agreement with the data. The resulting radio spectra 
for the three electron spectra are shown in Figure VI-10; the particular 
galactic parameters used in each calculation are given in Table B-3. 

It must be remembered that the galactic parameters used are 
not all independent quantities, e.g. Hjellming et al. (1969) show that 
if the clouds are in pressure equilibrium then the temperatures, 
neutral hydrogen densities, and electron densities for both the clouds 
and intercloud medium are uniquely related. However, not all the 
theoretical models predict the same values for the galactic parameters. 

In particular, we note that T^, the intercloud temperature, is highly 
uncertain and it accounts for a very large variation of the spectrum at 
low energies. We feel that our procedure for calculating a reasonable 
range of the interstellar electron spectrum is the best that presently 

i 

can be done . 


^ 5 * 

The electron spectra are plotted for energies between 70 MeV and 5 GeV. 
We have calculated the relative contribution to the synchrotron emission 
at v - 0-4 MHz from electrons of different energies. (The relative 
contribution to the intensity at v = 10 MHz is shown in Figure B-2.) 
Electrons with energies below 70 MeV contribute ~ 25$ of the emission at 
frequencies above 0.4 MHz. This percentage contribution is the same as 
the quoted absolute accuracy of the low-frequency radio data (Alexander 
et al., 1970). Similarly, we have chosen 5 GeV as the upper-limit 
energy since electrons with higher energy contribute ~ 25$ to the radio 
emission at frequencies below 600 MHz, the upper-limit frequency of the 
observations used in the analysis (Figure VI-8). 
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Power 


SPECTRUM 


LOW 


NOMINAL 


HIGH 


TABLE B-2 


■law Approximations of the Calculated Galactic Electron 

-Y 2 -1 

Spectra: j = AT (m -sec.sr.MeV) 


ENERGY RANGE 

A 

Y 

1 (MeV) 

/-s 

o 

I— 1 

X 



70-2000 

.0134 

1.80 

2000-5000 

2.75 

2.50 


70-2000 

.0254 

1.80 

2000-5000 

5.19 

2.50 


70-100 

3.45 

2.16 

100-150 

• 43.4 

2.71 

150-300 

1.38 x 10 5 

4.32 

300-500 

3.61 

2.47 

500-2000 

.0559 

1.8 

000-5000 

11.42 

2.5 
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TABLE B-3 

Galactic Parameters Corresponding to the Galactic 



Electron Spectra 

of Figure VI-9 



LOW SPECTRUM 

NOMINAL SPECTRUM 

HIGH SPECTRUM 

PARAMETER 

MODEL 

MODEL 

MODEL 

B (p-gauss) 

5 

5 

3 

l 




L(pc) 

7 

4 

4 


^i (pc) 

833 

50 

50 

* C (PC) 

5 

i 

1 

1 

T.(°K) 

10000 

10000 

3000 

' 1 

T (°K) 

; l 

250 

250 

70 

n^(cm ) 

.03 

.03 


n (cm "*) 

.02 

.02 
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11“1: Cross-section of the e^ 1 detector system 
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Figure II-3: Exploded view of the MOD-2 magnet and magnet guard 

counter. For clarity the origin of the coordinate 
system is shown displaced from its location at the 
center of the magnet gap. 
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Figure II-5: Gas Cerenkov counter light output as a function of 

energy for singly-charged particles. 
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Figure II -6: Magnetic flux density in. the gap of the analyzing magnet 

versus position. The curves represent the field components 
along three paths parallel to the z axis in the magnet gap. 
See Figure II-3 for definition of the coordinate system 
employed. 
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Figure II- 7: Geometrical factor of the MOD-2 detector versus 

particle rigidity. 



Figure II-8: Range of trajectory deflection angles in the MOD-2 

detector as a function of particle rigidity. The mean 
value, r.m.s. deviation (solid bar), and extreme values 
(dashed bar) of rigidity x deflection angle are shown for 
a random distribution of 1000 incident trajectories at 
each of 7 rigidities. 
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Figure III-l: Trajectories of three typical balloon flights. Invariant geomagnetic 

latitude contours are derived from the internal geomagnetic field 
only (Cain et al., 1967). 
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Figure III-2: Two representative time-altitude 

profiles. The curves are labeled 
with the flight number. 
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Figure III-3: Deep River Neutron Monitor counting rate versus time. 

The vertical bars indicate the times of balloon flights. 
The approximate periods of solar minimum and solar maximum 
are indicated by the horizontal bars. 




Figure III- 
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Figure IV-1: Typical event rate versus local time (Flight 71C2) . 

Negatrons (solid histogram) and positrons, (dotted 
histogram) are shown separately. Typical l-o error 
limits are indicated. The nighttime period used in 
the analysis is indicated. 
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Figure IV-2: Representative examples of the measured event rate versus 

atmospheric depth. The energy intervals indicated are 
those measured at the magnet. Also shown is the separation 
into primary and atmospheric secondary components as 
determined by the least-squares fitting technique 
described in the text. 

Dashed curve: best-fit primary contribution 

Dotted curve: best-fit secondary contribution 

Solid-curve: best-fit total positrons or 

negatrons. 

2 

The y probability, P, is indicated for each fit. 

Figure lV-2a: Energy range with a relatively large contribution of 

residual primaries at float altitude. 


Figure lV-2b: Energy range with essentially zero primary flux at 

float altitude. 
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Figure V-l: Primary positron and electron spectra (at 0 g/cm ) derived 

from the 1968 observations. 
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Figure V-2: Primary positron and electron spectra (at 0 g/cm ) 

derived from the 1969 observations. 
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Figure V-3: Primary positron and electron spectra (at 0 g/cm ) derived 

from the 1970 observations. 
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Figure VI-1: Selected near-Earth electron spectra for the period 1965- 

1971. The Caltech data are shown as filled squares (1968) 
and filled circles (1971). The Chicago data are represented 
by open diamonds (Fanselow et al., 1969), open squares 
(L'Heureux et al., 1972; Schmidt, 1972), and open circles 
and triangles (Schmidt, 1972). Data from the Goddard 
Space Flight Center (GSFC) experiment on the IMP-IV 
satellite are shown as crosses (Simnett and McDonald, 1969). 


For clarity the data have been separated into two graphs. 
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Figure VI-2: Illustration of the dependence of the modulated spectrum 

on the magnitude of the diffusion coefficient. Numerical 
solutions are calculated using an interstellar electron 
spectrum of the form 

j (T) = 2.70x10 T p/m^ • sec • sr *MeV 

03 r 

and a diffusion coefficient of the form 

R > R 

c 

R < R 

c 

Figure VI-2a: Calculated electron intensity at 1 AU versus kinetic 

energy for R^ = 440 MV . 

Figure VI-2b: Contours of constant phase - space density, F, for 

R =440 MV. 
c 

Figure VI-2c: Calculated electron intensity at 1 AU versus kinetic 

energy for R^ = 100 MV. A j a curve (dashed) is 

shown for comparison. 

Figure VI-2d: Contours of constant phase-space density for R c = 100 MV. 

Contours of constant \|f (r,T) (defined by equation VI-6) 
are shown as dashed lines. 

Figure VI-2e: Phase-space density contours derived from the diffusion - 

convection approximation for R^ = 100 MV. 
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Figure VI -2a and b 
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Figure VI-2c and d 
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Figure VI-3: Comparison of the 1 AU spectra derived from the force-field 

(FF) approximation, diffusion-convection (DC) approximation, 
and the numerical (FN) solution of the full transport 
equation for two different galactic spectra. The diffusion 
coefficient used in deriving each spectrum is described 
in the text . 

Figure VI-3a: Electron spectra at 1 AU derived from a galactic electron 

spectrum proposed by Meyer et al. (1971). 

Figure VI-3b: Positron spectra at 1 AU derived from a galactic positron 

spectrum calculated by Ramaty and Lingenfelter (1968). 
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Figure VI -3a 
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Figure VI -4: Radial and rigidity dependences of the diffusion coefficient used in calculating the 

numerical solutions of the transport equation shown in Figure VI-5 a . A boundary 
distance of 3 AU has been assumed for convenience. 
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Figure VI-5: 


Calculated electron spectra at 1 AU for different forms 
of K^( r ) and different values of the boundary distance, D. 

In each calculation the nominal galactic electron spectrum 
from the analysis of the non-thermal-radio-background data 
has been used with a power-law extrapolation below ~100 MeV 
(equation B-12) . The rigidity dependence of the diffusion 
coefficient used is given in Figure VI-4b. The magnitudes 
of the different radial functions K^(r) at r = 1 AU have been 

adjusted so that each calculated spectrum is derived using 
the same value of 



Vdr 

K(r,R) 


The near-Earth electron spectrum observed in 1968 (references 
in Figure VI-1) is shown for comparison. 


Figure VI-5a: The calculated spectra using the 6 different radial 

dependences of k shown in Figure VI-4a. 


Figure VI-5b: The calculated spectra assuming ^ independent of radius 

with assumed boundary distances of 3, 5, 10, 15, and 30 AU, 
respectively . 
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Figure VI-5b 
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Figure VI-6: Illustration of a calculation of the energy loss of 

positrons in diffusing from the boundary to 1 AU. The 
assumed unmodulated spectra are shown as solid lines. 

The corresponding near-Earth spectra, representing 
numerical solutions of the transport equation, are shown 
as dotted lines. The dashed curve indicates the galactic 
positron spectrum derived by Ramaty and Lingenfelter 
(1968). The diffusion coefficient used is independent 
of radius within a boundary of 3 AU; the rigidity dependence 
is given by equation VI-13. 
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Figure VI -6 
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Figure VI- 7: Calculated electron spectra at 1 AU for a diffusion 

coefficient of the form: 

R > R = 750 MV 
c 


(The complete description of the diffusion coefficient used 
is given by equation VI-30 with parameters of Table VI-2 for 
1968.) Both the numerical solution (FN) of the full transport 
equation and the diffusion convection approximate solution 
(DC) are shown. The positions of the relative "peaks" 
discussed in the text are shown. The values of these 
peak locations are: 

(Eq. VI-18, b = 1) 

(Eq. VI-18, b = 1/2) 

(Eq. VI-22) 

(Eq. VI -24) 

The near-Earth electron spectrum observed in 1968 
(references in Figure VI -1) is shown for comparison. 
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Figure VI -7 
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Figure VI-8: The frequency spectrum of the non-thermal- 

radio background in the galactic anticenter 
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Figure VI-9: Interstellar electron spectra derived from 

the non-thermal-radio-background data. The 
high and low spectra reflect the discussed 
range of galactic parameters. (Appendix B) 
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Figure VI-11: Interstellar positron spectra calculated 

from galactic nuclear collisions by different 
investigators. The near-Earth spectra measured 
in 1965-66 (Fanselow et al., 1969) and 1968 
(Caltech) are shown for comparison. 
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Figure VI-l2a-e: Electron modulation parameters, ty(l,T), for the period 

1965-1971. The data points are calculated from the 
DC approximation using the near -Earth electron data 
shown in Figure VI-1 and the calculated galactic 
electron spectra of Figure VI-9 (open circles - 
nominal galactic spectrum; upper filled circles - 
high spectrum; lower filled circles - low spectrum). 
The dashed lines indicate the approximate range of 
from the possible range of galactic spectra. In each 
figure the solid line represents the modulation 
parameter used (together with the nominal galactic 
electron spectrum) in deriving a numerical solution 
of the transport equation in agreement with the 
observed near-Earth spectra. The dotted lines in 
Figures VI-12a,b, and c correspond to limiting 
modulation parameter curves from which acceptable 
numerical solutions (using the same nomina 1 galactic 
electron spectrum) can also be derived. 
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Figure VI-14: Positron modulation parameters, ^(1,T), for the period 

1965-1971. The data points refer to the ^(l.T) calculated 
from the DC approximation using the spectra observed near 
Earth (1965-66, Fanselow et al., 1969; 1968-1971, Caltech) 
and the interstellar positron spectra of Ramaty and 
Lingenfelter (1968) and Beedle (1970). Above ~100 MeV 
the solid and dashed lines are the electron modulation 
parameters from Figure VI-12. In Figure VI -14b the 
modulation parameter represented by the solid line was 
used to derive the numerical solutions of the transport 
equation shown in Figure VI-16a and b. The dotted lines 
in Figure VI-14b correspond to the possible range of 
positron modulation parameters which was used to derive 
the range of the interstellar electron spectrum at low 
energies (%0 MeV) shown in Figure VI-17. 
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Figure VI -14a 
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Figure VI -14b 
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Figure VI-14d 
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Figure VI-14e 
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Figure VI-15: Calculated and measured positron fractions as a function 

of energy for two different models of the galactic 
electron spectrum and diffusion coefficient. The 
calculated fractions at the boundary (solid) and at 1 AU 
(dashed) are shown. Above 100 MeV the diffusion 
coefficient for both model 1 and model 2 is given by 
equation VI-30 with the parameters listed in Table VI-2 
for 1968. The galactic positron spectrum of Ramaty and 
Lingenfelter (1968) is assumed in the calculations. 

MODEL 1: The nominal galactic electron spectrum from 

analysis of the non-thermal-radio-background 
data (Section VI.E.l) is used with a power- 
law extrapolation below 100 MeV (equation 
B-12). The diffusion coefficient for all 
energies is described by equation VI-30. 

MODEL 2: The interstellar electron spectrum shown in 

Figure VI-16b (solid line) is used. The 
rigidity dependence of the diffusion 
coefficient is described by equation VI- 30 
except that below 60 MV k(R) is assumed 
proportional to 1/R. 
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Figure VI-16: Positron and electron spectra derived from the modulation 

parameter of Figure VI-14b (solid line). Below 60 MV, 

K (R) OL 1/R. Above 60 MV, ^ (R) is given by equation VI-30 
using the parameters listed in Table VI-2 for 1968. 


Figure VI-16a: The calculated positron spectrum at 1 AU derived from a 

galactic spectrum calculated by Ramaty and Lingenfelter 
(1968) . Hie near-Earth spectrum observed in 1968 
(Caltech) is shown for comparison. 


Figure VI-16b: The calculated electron spectrum at 1 AU derived from 

the indicated galactic electron spectrum. In order that 
the calculated and observed spectra at 1 AU agree 
(references fur observations in Figure VI-1), the 
nominal galactic electron spectrum cannot be 
extrapolated by a power -law below ~100 MeV but must be 
modified as shown. 
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Figure VI-17: Approximate range of the interstellar electron 

spectrum. The shaded region (bounded by dotted 
lines) below ~50 MeV indicates the range of de- 
modulated electron data using the limiting positron 
modulation parameters shown as dotted lines in 
Figure VI-14b. The shaded region above ~70 MeV 
is the approximate range from the analysis of the 
non-thermal-radio-background data (Section VI.E.l). 
For comparison the assumed galactic electron 
spectrum of Figure VI-16b is shown as the solid 
line . 
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Figure VI-18: Observed power spectra of the inter- 

planetary magnetic field. 
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Figure VI-19: Parallel diffusion coefficients calculated from the 

power-spectra data (Figure VI-18) of Jokipii and 
Coleman (1968) (solid line), Bercovitch (1971) 
(dotted line), and Quenby and Sear (19 71) (dashed 
line) . 
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RIGIDITY (GV) 


Figure VI-20: Radial diffusion coefficients derived from the 1968 

electron modulation study and from the magnetic power 
spectra of Quenby and Sear (1971) (12/68-3/69). Two 
power-spectra-derived curves are shown. Curve 1 is 
derived under the assumption « K 1( > and curve 2 

corresponds to the case k x = 4 x 10 ^ cm2/sec. The 

error bars on curves 1 and 2 correspond to a 2<? un- 
certainty in the power-spectra data. The electron 
modulation result is derived from the modulation param- 
eter of Figure VI-12b (solid line) assuming k in- 
dependent of radius with a boundary at 12 AU. The three 
error bars correspond to the limiting ^(1,T) (dashed 
lines in Figure VI-12b) based on the possible range of 
galactic electron spectra (Figure VI-9). 
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BOUNDARY DISTANCE, D (AU) 

Figure VI-21: Comparison (at 1 GV) of the radial diffusion coefficient 

from Figure VI-20 as a function of boundary distance, D. 
The bands reflect the uncertainties indicated in 
Figure VI-20. The power-spectra results (curve 1 was 
derived assuming « k and curve 2 corresponds to the 

case = 4 x 1021 cmVsec) are independent of the 

assumed boundary distance. The modulation result is 
derived assuming k independent of radius inside D. The 
crosshatched areas indicate the range of boundary 
distances required for consistency between the two 
diffusion coefficients for the case «;^(r) = constant. 
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Figure VI-22: Limits on boundary distance, D, for various 

radial dependences of the radial diffusion 
coefficient, k • The minimum and maximum D 
are plotted for different values of the power - 
law index n for the case k « K • The 

i II 

horizontal bar at n = 0 indicates the range 
6-15 AU obtained for independent of radius 
(see Figure VI-21). 
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Figure VI-23: Comparison of the measured and calculated proton spectra 

at 1 AU for the time periods shown. The same interstellar 
spectrum has been used in. deriving each calculated 
spectrum. The numbers associated with the calculated 
curves refer to entry numbers in Table VI-3. The data 
collected with Caltech instruments are shown as filled 
circles (Garrard, 1973). Data from other references are 

Figure VI-23a: Open circles - Fan et al. (1966) 

Crosses - Ormes and Webber (1968) 

Triangles - Fan et al. (1968) 

Open diamond (for solar minimum in 1954)- 
McDonald (1958) 

Figure VI- 23b: Open squares - Lezniak and Webber (1971) 

Figure VI-23c: Open squares - Hsieh et al. (1971) 


Note that the low-energy portion of the interstellar 
spectrum is shown as a dashed line. Due to adiabatic 
deceleration in the interplanetary medium the calculated 
spectrum at 1 AU is insensitive to the interstellar 
intensity below ~100 MeV. 
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Figure VI- 24: Comparison of measured and calculated He nuclei spectra 

at 1 AU for the time periods shown. The same inter- 
stellar spectrum has been used in deriving each 
calculated spectrum. Hie numbers associated with the 
calculated curves refer to entry numbers in Table VI-3. 
The data collected with Caltech instruments are shown as 
filled circles (Garrard, 1973). Data from other 
references are: 

i 

Figure VI-24a: Open circles - Fan et al. (1966) 

Crosses - Ormes and Webber (1968) 
Triangles - Fan et al. (1968) 

Figure VI -24b: Open squares - Lezniak and Webber (1971) 

Figure VI-24c: Triangles - Mason (1972) 

Note that the low-energy portion of the interstellar 
spectrum is shown as a dashed line. Due to adiabatic 
deceleration in the interplanetary medium the calculated 
spectrum at 1 AU is insensitive to the interstellar 
intensity below ~100 MeV/nucleon. 
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Figure A-l: Schematic view of a particle trajectory seen 

in projection. Parameters used in trajectory 
self-consistency checking are shown. 
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Figure A-2: Measured distributions of the trajectory parameter ^ using 

the MOD-1 detector. Also shown are the theoretical 
Gaussian distributions using the standard deviation of 
equation A-3a. The crosshatched areas are the rejection 
zones according to equation A-4b. 


Figure A- 2a: Mono-energetic beam of 790 MeV positrons. 

Figure A-2b: Mono-energetic beam of 85 MeV positrons. 
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Figure A-2a 



Figure A- 2b 
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(a) , (b) 

MOMENTUM (MeV/c) x SCATTERING ANGLE (RAD'ANS) 

Figure A-3: Calculated electron scattering-angle distribution. 

Curve 1 is the angular distribution of the projected 
scattering angle. As plotted, the curve is valid 
for all electron momenta above a few MeV/c. Curve 
2 is a Gaussian distribution with r, = .60 MeV/c 
radians. Only one sign of the scattering angle is 
shown since the distribution is symmetric about 
zero. 
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DEFLECTION (MILLIRADIANS) 


Figure A-5: Distribution of small deflections measured 

during the nighttime period of flight 71C2. 
The predicted angular distributions of 
cosmic-ray protons (solid curve-solar 
maximum spectrum; dashed curve-solar minimum 
spectrum) are shown for the instrument 
angular resolution function with 
a . = .0025 radians. 




(Ml LLIRAD I ANSI 

2 

Figure A-6: ^ versus a for fitting the calculated 

angular distribution to the measured 
distribution of cosmic-ray protons 
shown in Figure A-5. The proton spectrum 
appropriate for solar maximum has been 
assumed . 
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Figure A- 7: Deflection resolution P, FWHM, of the detector versus 

rigidity. The data points represent measurements in 
a mono-energetic positron beam using the MOD-1 
detector. Calculated curves for the detector in 
both the MOD-1 configuration (solid) and MOD-2 
configuration (dashed) are shown. 
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used in calculating the flux of negatrons and 
positrons produced in the gas Cerenkov counter. 
The gamma-ray spectra from the calculations 
of Beuermann (1971) resulting from cosmic-ray 
nuclei alone (dashed) and nuclei + electrons 
(solid) are shown. 
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Figure A-12: Response functions, Rj_(T), for the MOD-2 detector system. Each curve represents 

the probability that a particle with kinetic energy T will be observed in the 
ith energy interval. The dotted vertical lines indicate the end-points of the 
energy intervals described in Table IV-1. 
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Figure B-2: Relative contribution to synchrotron intensity 

I at v = 10 MHz from cosmic-ray electrons of 
different energies. The nominal galactic 
electron spectrum (equation B-12) has been 
used . 
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Figure B-3: Correspondence between radio frequency and electron 

energy. The energy at which electrons make the maximum 
contribution to the synchrotron intensity at the 
frequency v is plotted for two different values of the 
magnetic field strength. In each calculation the nominal 
galactic electron spectrum (equation B-12) has been 
assumed . 
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Figure B-4: Relative variation of interstellar electron spectrum for 

the range of galactic parameters discussed in the text. 
F(W) is the ratio of the calculated interstellar electron 
intensity to the nominal interstellar intensity 
(equation B-12) at electron energy W. The nominal set of 
parameters is indicated by the bracket. Each labeled 
solid curve is calculated by changing the value of only 
the indicated parameter from the nominal set. The 
dashed lines correspond to the assumed range of variation 
used in computing the high and low galactic electron 
spectra shown in Figure VI- 9. 
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TOTAL ENERGY, W (MeV) 

Figure B-4 




